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Bei Bearbeitung der nachfolgenden astronomischen Probleme 
ist der Verfasser vielfach von neuen Gesichtspunkten ausgegangen. 
So wurde bereits die einfache Planetenbewegung, abweichend von 
der traditionellen Behandlungsweise, durch das Verfahren der 
Variation der Parameter begrü;ndet; in der mechanischen Qua- 
dratur der Störungen, teilweise auch in der Bahnbestimmung 
fand zur Berücksichtigung der Glieder höherer Ordnung eine 
aus der Eule.r sehen Reihe sich sehr einfach ergebende Reduk- 
tionsmethode Anwendung. Überall wurde jedoch aus didaktischen 
Gründen — da die Schrift auch dem Anfangsstudium dienen 
soll — ein möglichst enger Anschlufs an die gebräuchlichen, 
insbesondere die Encke -Lagrang eschen Rechnungsformen zu 
erreichen gesüpht. Auch glaubte der Verfasser aus denselben 
Gründen mehr Gewicht auf klare und bündige Darstellung des 
Gedankengangs legen zu müssen als auf eine erschöpfende Ent- 
wickelung des analytischen und numerischen Details, die zweck- 
mäfsiger dem Spielraum der eigenen Thätigkeit des Lesers 
überlassen bleibt. 

Frankfurt a, M., November 1889. 

Der Verfasser. 
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A. 
Theorie der elliptischen Bewegung. 



Einleitung. 

Das unter dem Namen „Methode der variierenden Kon- 
stanten" bekannte Verfahren der Integralrechnung leistet auf 
astromechanischem Gebiete in allen denjenigen Fällen unschätz- 
bare Dienste, in welchen die Veränderung (Variation) der Kon- 
stanten — wie beispielsweise die Änderung der elliptischen Bahn- 
elemente eines Planeten — nur langsam mit der Zeit fort- 
schreitet. Denn unter dieser Voraussetzung werden die auf- 
tretenden Differentialausdrücke stets einer, wenn auch nur empi- 
rischen, Integration fähig. 

Zur Herstellung allgemeiner endlicher Integrale scheint das 
Verfahren der Variation in der Mechanik der Himmelskörper 
bis jetzt keine Verwendung gefunden zu haben. Die eingehende 
Erörterung eines hierher gehörigen Falls — nämlich der Ent- 
wickelung der Integrale der elliptischen Bewegung aus denen 
der oskulierenden geradlinigen Bewegung — dürfte deshalb wohl 
von einem allgemeinern Interesse sein. Zwar ist für die Er- 
forschung der elliptischen Bewegung schon in anderer Weise 
hinreichend gesorgt. Da indessen das Problem der drei Körper, 
also der gestörten elliptischen Bewegung, stets im Anschlüsse an 
das Problem der zwei Körper und dann meist, namentlich in 
Fällen der Praxis, durch die Variation der Konstanten behandelt 
zu werden pflegt, so liegt schon in der immer zu erstrebenden 
Einheitlichkeit der üntersuchungsmethode eine Aufforderung, 
auch dies letztere Problem wieder an ein einfacheres anzu- 
schliefsen, also das Problem der zwei Körper auf das eines ein- 
zigen Körpers zurückzuführen. Es bietet so die Lösung des ein- 
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2 A Theorie der elliptischen Bewegung. 

fachern Falls bereits den Schlüssel zur Lösung des zusammen- 
gesetzten Problems, welches den zweiten Gegenstand der vor- 
liegenden Abhandlung bilden wird. 

Überdies werden wir uns in der Folge überzeugen, dafs 
das den Kepler sehen Gesetzen entsprechende Problem der zwei 
Körper bei Anwendung der Methode der variierenden Konstanten 
sich nicht blofs unter einem neuen, sondern auch unter einem 
erheblich einfachem Gesichtspunkte darstellt. 

Was die Variationsmeihode selbst betrifft, so kann man 
zwei Formen von wesentlich verschiedener Beschaffenheit, eine rein 
mathematische und eine ausgesprochen dynamische, unter- 
scheiden. Die erstere wird zunächst stets bei mathematischen 
Problemen angewendet, aber auch bei astronomischen Entwicke- 
lungen fast bis zum völligen Ausschlüsse der zweiten Form 
bevorzugt. In der folgenden Untersuchung werden wir hingegen 
lediglich von der dynamischen als der bei weitem geschmeidigeren 
und energischeren Form Gebrauch machen. 

Nur der Vollständigkeit wegen wird auch die mathematische 
Form eine kurze, aber zu ihrer Charakteristik hinreichende Be- 
sprechung erfahren. 

IL Allgemeine Betrachtungen 
über die Methode der yarilerenden Eonstanten. 

Jede Bewegung kann und mufs in ihren unendlich kleinen 
Teilen als einfach , d. h. als geradlinig und gleichförmig gedacht 
werden — mag im übrigen der Bewegungsvorgang, als Ganzes 
betrachtet, ein noch so zusammengesetztes Phänomen bilden. 
Nach dieser Auffassung erscheint deshalb auch die verschlungene, 
in ihrem Gesamtverlaufe unfafsbare Planetenbewegung als eine 
Summe geradliniger Elementarbewegungen. Wäre man imstande 
die Einwirkung aller einen Planeten angreifenden äufseren Kräfte 
zu irgend einer Zeit aufzuheben, so würde derselbe von da ab, 
zufolge d9s Trägheitsgesetzes, in der Richtung der Bahntangente 
mit der ihm gerade innewohnenden Geschwindigkeit seine Be- 
wegung fortsetzen. 

Nehmen wir an, die Bewegung vollziehe sich in einer Ebene 
— ein Fall, der unserer Betrachtung zunächst unterliegt — so 



AUgem. Betrjujhtungen üb. d. Methode d. variierenden Konstanten. 3 



können wir den Ort des bewegten Punktes zu einer bestimmten 
Zeit durch die rechtwinkligen Parallelkoordinaten 

und die momentanen Geschwindigkeiten in der Richtung der 
Koordinatenachsen durch 

dt ' dt 
ausdrücken. Nach Verlauf 
der Zeit t wird dann der 
dem Einflüsse äufserer 
Kräfte entrückte Punkt den 
Ort 



X ^ x^-\- 



dxp 
dt 



t 



einnehmen. 



dpi 



?-'t 



dt 
Durch Elimi- 
nation von t ergiebt sich 
hieraus die Gleichung der Bahn: 




y - 



dx^ 



x-\^ 



ypdx^ — x^dy^ 
dx^ 



demnach die Gleichung einer Geraden , sowie die G eschwindig- 
keit in der Bahn 

dt y [dt )'^\dt ) ' 

Setzt man zur Abkürzung 



dyp 
dxa 



m und 



y^dx^-'X^ dy^ 
dxo 



- 2» 



so hat man m und ^ als Konstante oder Elemente der 
Bahn zu betrachten. Es verdient hier besondere Beachtung, 

dafs diese Elemente lediglich durch ^o^yo^-^y -1^, also durch 

die Koordinaten und ihre Geschwindigkeiten ausgedrückt sind. 
Werden vermöge der Gleichungen 

a; «= r cos 1// ; y = r sin 1// 
Polarkoordinaten eingeführt, so erhält man als Gleichu ng der Bahn : 



r = 



sini// — mcosxp 



(I)- 
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Machen wir nun aber die Annahme, dafs bei Beginii der Zeit t 
gewisse äufsere Kräfte anfangen oder fortfahren auf den Punkt 
einzuwirken, so wird derselbe zwar auch im folgenden Zeit- 
differentiale eine geradlinige Bahn beschreiben, Indessen trerden 
die Elemente oder Konstanten dieser neuen Bahn — infolge des 
Einflusses der äufeeren Kräfte — von den früheren etwas ver- 
schiedene Werte annehmen. Im dritten Zeitelemente wird sich 
dieser Vorgang wiederholen und sofort durch alle folgende Zeit- 
teilchen. 

Besäfsen wir nun ein Mittel, das Gesetz der Veränder- 
lichkeit der Konstanten ausfindig zu machen, mit anderen 
Worten, letztere als bekannte Funktionen der Zeit t darzustellen, 
so wäre der Ort des Punktes für alle Zukunft als gegeben 
anzusehen. Man könnte dann sagen: 

Die Bahn des Punktes sei durch Bestimmung 
der Variation der Konstanten gefunden worden. 

Die Gerade, welche dabei den Ausgangspunkt unserer Be- 
trachtung bildete und bei Beginn der Zeit t mit der Bahn selbst 
einen Augenblick zusammenfiel, liefse sich, nach Analogie eines 
in einem verwandten Falle eingeführten Sprachgebrauchs, als 

oskulierende Gerade 

bezeichnen, d. h. als diejenige Gerade, welche sich im betrachteten 
Augenblicke der Bahn am engsten anschliefst. 

Wir haben im Vorhergehenden angenommen, die Bewegung 
gehe in einer Ebene von statten. Es leuchtet indessen ohne 
Weiteres ein, dafs dieselbe Auffassung auch auf jede, beliebig im 
Räume verlaufende Bewegung übertragbar ist. Nur haben wir 
alsdann den ganzen Vorgang auf ein dreiachsiges Koordinaten- 
system zu beziehen und der oskulierenden Geraden die Gleichungen 

beizulegen, sowie die oskulierende Geschwindigkeit 
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ZU setzen. 

Wählen wir ein anderes Beispiel. In jedem Augenblicke 
hat ein Planet einen bestimmten Ort, eine bestimmte Bevegungs- 
richtung und Geschwindigkeit. Wenn nun von eiuem beliebigen 
Zeitpunkte an — mit Ausschlufs der Attraktion anderer Planeten — 
lediglich die zwischen ihm und der Sonne bestehende Anziehungs- 
kraft thätig wäre, so läfst sich zeigen, dafs die Bahn die Gestalt 
eilies Kegelschnittes, sagen wir einer Ellipse, annimmt und dals 
diese Bahn, so lange keine neuen Kräfte hinzutreten, nach Gröfse, 
Lage und allen anderen Beziehungen unverändert fortbesteht. 

In ähnlicher Weise, wie früher aus den für den Oskulations- 
Z€|itpunkt gegebenen Koordinaten und Geschwindigkeiten die 
E}^mente der oskulierenden Geraden sich ergeben, lassen sich 
jetzt aus 4enselben Gröisen 

^ .. ^ ^0 dyo äz^ 

die sechs Elemeifte der momentanen Ellipse, deren einer Brenn- 
punkt in der Sonne liegt, entwickeln. Dabei können wir, der 
Einf^hheit wegen und ohne Hintansetzung der Allgemeinheit, 
st^ts d^ Ursprung der Koordinaten mit dem Sonnenmittelpunkt 
zusammenfallen lassen. Der bemerkte Zusammenhang zwischen 
den 6 Bahnelementen und den 6 Konstanten der momentanen 
Bewegung ist ein für die klare Auffassung des Problems höchst 
wichtiger Umstand, dessen Begründung wir uns für eine spätere 
Gelegenheit vorbehalten. 

Diese sechs Elemente — aus denen der rein elliptische Ort 
fjir jede gegebene Zeit bestimmbar ist — werden wir in Zukunft 
durch folgende Buchstaben bezeichnen (s. Fig. 2, in welcher S 
die Sonne und P den Planeten bedeutet) : 

1. die gjrolse Halbachse (mittlere Entfernung) durch a, 

2. die Excentricität durch e, 

3. die wahre Anomalie zur Zeit der Oskulatiqn durch <p^, 

4. den Winkel vom Perihel bis zur Schnittlinie (Knotenlinie) 
der Bahnebene mit der Ebene X-Y (etwa der zur Zeit 
der Oskulation stattfindenden Ekliptikebene), gezählt in 
(Je;* Richtung der Bewegung, durch w, 
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5. den Winkelabstand dieser Knotenlinie selbst von ein|^ 
festen Geraden in der X-F- Ebene (etwa der momen- 
tanen Äquinoktiallinie) durch ft, 

6. die Neigung der Bahnebene gegen dieselbe Fundamental- 
ebene durch % 

wozu dann noch 

7. die mittlere tägliche Bewegung n tritt, welche indessen 
zufolge des dritten Kepl ersehen Gesetzes aus der Gaufs- 
schen Konstanten und der Halbachse a gefunden, deshalb 
nicht als selbständiges Element betrachtet werden kann. 

Durch die Länge Sl des aufsteigenden Knotens, der übrigens 
wie gelegentlich bemerkt sei, auch selbst durch Sl bezeichnet 

zu werden pflegt, und 
Fig. 2. die Neigung i gegen 

die Ekliptik ist die 
Lage der Bahn- 
ebene, durch den 
Winkel a> die Lage 
des Perihels gegen 
den aufsteigenden Kno- 
ten, durch die wahre 
Anomalie <f^ der Ort 
des Planeten in der 
Epoche, endlich durch 
a und e sind die Di- 
mensionen der Ellipse 
gegeben. 

Lassen wir nun 
unsere frühere ein- 
schränkende Voraussetzung wieder fallen und nehmen — wie es 
der Wirklichkeit und dem Gesetze der allgemeinen Gravitation 
entspricht — an, dafs nicht blofs zwischen Sonne und Planet 
eine wechselseitige Anziehung bestehe, sondern dafs diese beiden 
Himmelskörper auch mit allen anderen Planeten durch gleiche 
Kräfte verbunden seien, so kann man nur behaupten, dafs der 
Planet ein unendlich kleines Teilchen jener Ellipse wirklich 
beschreibe, im zweiten Zeitelement aber — infolge der hinzu- 
tretenden Einwirkung der anderen Planeten — zwar auch wieder 
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ein unendlich kleines Teilchen einer Ellipse durchlaufe, deren 
Elemente jedoch ein wenig verschieden sind von denen der 
ersten Ellipse und so in allen folgenden Zeitteilchen. Man könnte 
also dem wahren Laufe des Planeten folgen, wenn man für den- 
selben — unter Zugrundelegung jener ersten Elemente — zunächst 
von einer Ellipse ausginge, dann aber diesen Elementen an- 
gemessene Änderungen zuerteilte. Vermöchte man nun wiederum 
diese Änderungen als Funktionen der Zeit zu bestimmen, so wäre 
damit der Ort des Planeten für alle Folgezeit durch einfache 
Rechnung zu ermitteln. Auch in diesem Falle hätte man die 
veränderliche Planetbahn aus der ursprünglichen, s. g. osku- 
lierenden Ellipse durch Bestimmung der Variationen ihrer 
Konstanten abgeleitet. 

Obwohl die Erzeugung der wahren Planetenbahn aus den 
oskulierenden Geraden einfacher und naturgemäfser erscheint 
als die Entwickelung aus der oskulierenden Ellipse, so bietet 
doch die letztere Auffassung in der Astronomie entschiedene Vor- 
teile. Denn die anzubringenden Änderungen der oskulierenden 
Ellipse — welche ja nur der im Vergleiche zur Sonnen- 
anziehung sehr geringfügigen Planetenanziehung ihr Dasein ver- 
danken — gehen äufserst langsam vor sich und schwanken meist 
zwischen sehr kleinen Werten hin und her, während die zur 
Darstellung eines Planetenlaufs erforderlichen Korrektionen der 
oskulierenden geradlinigen Bewegung alsbald und mit sehr 
beträchtlichem Werte hervortreten. Dieser Umstand würde zwar 
an sich kaum für die eine oder die andere Methode ausschlag- 
gebend sein, wenn es in unserer Macht stände, die Variationen 
der Konstanten in beiden Fällen allgemein und streng als 
Funktton der Zeit zu entwickeln. Da dies indessen nicht oder 
doch nur in sehr beschränkter Weise der Fall, man vielmehr 
der Regel nach darauf angewiesen ist, mit Hilfe der s. g. mecha- 
nischen Quadratur für kleinere Zeitabschnitte die Abweichungen 
von der oskulierenden Bahn zu bestimmen, so wäre es gewifs 
recht irrationell bei Anwendung dieses empirischen Verfahrens 
von der oskulierenden Geraden auszugehen, obwohl dieselbe 
uns gewissermaafeen von der Natur selbst dargeboten wird. Es 
gebührt vielmehr demjenigen Wege der Vorrang, welcher die 
geringsten Korrektionen verlangt und die grofse Kluft zwischen 
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der gleichförmigen geradlinigen und der wahren Bewegung der 
Planeten zunächst durch ein allgemeines Verfahren möglichst 
aufhebt 

Wenn nun aber auch die wahre, gestörte Planetenbewegung 
nicht generell aus der gleichförmigen Tangentialbewegung dedu- 
zierbar ist, so kann doch andererseits die Frage aufgeworfen 
werden, ob nicht gerade jener erste und bedeutendste Schritt 
zur wahren Bahn gleichfalls mit Hilfe der variierenden Kon- 
stanten ausgeführt, d.h. die rein elliptische Planetenbewegung 
aus der oskulierenden geradlinigen Bewegung allgemein durch 
Variation der Eonstanten entwickelt werden kann. Im bejahen- 
den Falle wäre dann doch durch blofse Veränderung der Para- 
meter die definitive Planetenbewegung aus der momentanen 
Tangentialbewegung hergeleitet — wenn auch in zwei Annähe- 
rungsstufen, insofern zunächst aus der geradlinigen Bewegung 
die elliptische und sodann aus dieser die eigentliche gestörte 
Bewegung lediglich durch Anwendung unseres Prinzips erschlossen 
worden wäre. 

Es ist bereits in der Einleitung darauf hingewiesen worden» 
dals, abgesehen von den theoretischen Vorzügen dieses Ver- 
fahrens, die in der That mögliche Ableitung der elliptischen aus 
der tangentialen Bewegung sich auch ungleich einfacher gestaltet 
als die allgemein übliche Herleitung der Ellipse aus den Diffe- 
rentialgleichungen der Bewegung. Im letzteren Falle werden 
bekanntlich aus den Differentialgleichungen zweiter Ordnung 
durch geeignetes Kombinieren und einmaliges Integrieren zwei 
Differentialgleichungen erster Ordnung (das s. g. Flächenprinzip 
und das Gesetz der linearen ^Bahngeschwindigkeit) gefunden und 
aus diesen dann erst durch nochmalige Integration die Gleichung 
der Bahn zunächst in Form einer cyklometrischen Funktion — 
erhalten. Dabei veranlafst die Bestimmung der Konstanten noch 
eine recht umständliche Nebenuntersuchung. — Bedient man 
sich hingegen des Variationsprinzips, sp liefert bereits eine ein- 
malige überaus leichte Integration sofort die Polargleichung der 
Bahn, während sämtliche Konstanten aus den Bedingungen des 
Problenas ebenfalls in ungezwungener Weise hervorgehen. Über- 
haupt verdient als eine Eigentümlichkeit der Methode der 
v^,riierendeft Koft^t^nten gleich m dieser St^eUe betont zu werdep, 



Analytischer Ausdruck der dyuamischea Vaiiatlonsmetliode. 9 

dafs durp^ m 4er I^qgel naoh ^ur emfi^che, keine ^OBpelte« 
Integrationen beding wepdei^. 

in. Analytisclier Ausdruck der dynamischen Yariations- 

methode. 

Ehe wir ^ur analytischen Darstellung der Methode schreiten, 
dürfte es zweckmäfeig sein, die wichtigsten, der folgenden Unter- 
suchung zu Grunde liegenden Relationen zunächst in Erinnerung 
zu bringen. 

Dafs — bei Betrachtung einer ebenen Bewegung und 
Annahme eines zweiachsigen Koordinatensystems — für 
den Radiusrektor r die Gleichung besteht (s. Fig. 3). 

sei hier nur deshalb erwähnt, weü hieraus sofort für das Diffe- 
rential des Vektors folgt: 

dr « --dx-{--dy. 
r r ^ 

Nennt man ferner ds das Bahnelement, 
so hat man 

(fc2 _ dx^-^dy\ 
Wird das Winkelelement mitdi/^beaseich- 
net, so ist 

rdip 
der Ausdruck des djöm Radius r ent- 
sprechenden Kreisbogenelements. Es 
milst demnach r^dip 

die doppelte Fläche des unendlich kleinen Bahnsektors 8MN. 
Dieselbe Fläche kann aber (wenn man erwägt, dafs 

{rdipy =f fi^^-dr^ 
mit Hilfe der obigen Beziehungen auch leicht in die Form ge- 
bracht werden: 

xdy — y dx^ 
also nicht, zur Vermeidung* negativer Flächen in die sonst gleich- 
wertige Form 

ydx—xdy. 
EndUch hat man nach dem Gravitationsgesetze für die im Radius- 
vektor zwischen Sonne und Planet wirkende Anziehungskraft — 
yoirausgesetzt^ dafs mit f die Kraftei^heit (also mit yf'^k die 
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s. g. Gaufssche Konstante), mit 1 die Sonnen- und mit m die 
Planetenmasse angedeutet wird, den Ausdruck 

/(l-f m) 

welcher mit dem negativen Vorzeichen zu versehen ist, weil beide 

Kräfte, nämlich die Sonnenkraft -^ und die Planetenkraft ^-, 

den Radiusvektor zu vermeiden streben. 

Zerlegt man dieselbe nach der Richtung der x und y, so 
ergeben sich für die beiden Komponenten die Gleichungen: 

^--^cosv-"-^ (1) 

^^^-f^^rp^-fjl^ (2) 

Es sind dies die Differentialgleichungen der Bewegung, wenn 
lediglich die Attraktion von Sonne und Planet in Betracht ge- 
zogen wird. 

Im Anschlüsse an dieselben können wir unsere Aufgabe 
nun folgendermafeen formulieren: Man kennt die Integrale der 
Gleichungen: 

^ 0^ ^"^' 

also die Integrale für den Fall, dafe überhaupt keine beschleu- 
nigenden Kräfte den bewegten Punkt angreifen, nämlich: 

woraus — wenn zur Abkürzung 

gesetzt wird — sich ergiebt: 



- 2 



* Nicht zu verwechseln mit der Planetenmasse w, welche nur in der 
Verbindung fiy^^) auftritt. 
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ff -^ mx-{-q 
oder in Polarkoordinaten 

^ __ 2 



Biniff — mcosiff (I) 

und ds //dyöy. /q?ä?öy" 

dt^y [dt )'^\dt I' 

Es entsteht nun die Frage: Wie gestaltet sich das Integral 
der Gleichungen (er), wenn an Stelle der rechtseitigen Null die 
relativen Anziehungskräfte treten, demnach die Gleichungen (a) 
übergehen in: 

d^x _/a±m)£ \ 

dt^ "" ^* l (ß) 

d^ fn^m^y f "^^^^ 

dt^ ^ r» J 

Torausgesetzt, dafs die allgemeine Form des Integrals (I) 
beibehalten werden soll. 

In der Analysis nennt man die Gleichungen (a) reduzierte, 
die Gleichungen (ß) TollstSndige lineare Gleichungen der zweiten 
Ordnung, so dafs man unserer Aufgabe auch die Fassung geben 
könnte: Man kennt das Integral zweier simultanen reduzierten 
Differentialgleichungen und sucht das Integral zweier vollstän- 
digen Gleichungen — unter der Voraussetzung, dafe beide Inte- 
grale eine gemeinsame Form haben. 

Es ist ohne Weiteres einleuchtend, dafs die erste Form mit 
den Bedingungen des zweiten Falles nur dann verträglich ist, 
wenn die ursprünglichen Konstanten m und q als variabel be- 
trachtet werden, so dafs es also darauf ankommt, das Gesetz 
der Variation von m und q festzustellen. Von diesem Gesichts- 
punkte aus läfst sich das Verfahren sehr wohl mit der bekannten 
Methode der unbestimmten Koeffizienten vergleichen, wo gleich- 
falls die allgemeine Form der Entwickelung bekannt ist und nur 
die unbestimmten Koeffizienten — in unserem Falle m und q — 
den Bedingungen gemäfs zu ermitteln sind. — Die höhere Ana- 
lysis giebt auch bestimmte Vorschriften über den bei Lösung 
unserer Aufgabe im allgemeinen einzuschlagenden Weg, von denen 
beispielsweise Encke in seiner berühmten Abhandlung über die 
Berechnung specieller Störungen Gebrauch gemacht hat (vergl. 
Berliner Astronom. Jahrbuch für 1837 u. 38). Dies, der reinen 
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Mathematik entlehnte Verfahren erwe^t sich indessen in der 
Praxis als schwer durchführbar, wenigs|;.eii9 bedingt fjsfsßlbe — 
wie die Entwickelungen von Ejicke u. A. beweisen — ganz 
übermäfeige Nebenrechnungen. 

Eine kurze Andeutunc, up aiff einem anderen Wege z[\ 
den Variationen der Konstanten zu gelangen, findet man in den 
Elßwentei^ der Mechanik des Elimmels vpp Möbius, p. 310. 
Die überraschende Einfachheit, mit d^,r es mir gelang ai;f Grund 
dieses Prinzips die Encke sehen Formeln für die Vapationei^ 
der elliptischen Bahnelemente herzuleiten, veranlaf^te mic^i einen 
allgemeinen analytischen Ausdruck für dies Verfahren aufzustellen, 
und zwar durch folgende Bet^achtiijng: 

Da beschleunigende Erä^e stets in der Fqrn^ von Differential- 

quotienten der zweiten Ordnung gegeben sind, z. B. -=^, so 

werden durch die Integration je zv^i Kpnstanten eingefiHirt. 
Nehmen wir ^^n der Einfach^xeit lyegen an, es sei die Wirkung 
der beiden Ji^omponenten 

d^x .d^y 

W^'^'dP 
zu untersuchen — ein Fall, der uns hier vorzugsweise angeht — 
so ergebest sich vißr Integrationsko^stantep. Dieselben sii^4 durch 
vier Bedingungen bestimmbar, also jedenfalls dann, wenn für 
irgend ein.en Zeitpunkt die Koordinaten 

^, y 

und deren Geschwindigkeiten -^ und -^ gegeben sind. Es mufs 

demnach beispielsweise zwischen der Konst^ntßp a und den ge- 
nannten Gröfsen eine Relation b^tehen, die wir allgemein durch 

K"'' ^' $' -f ' ") " ^ (^> 

aussprechen. 

Wirken nun auf den bewegteil Punkt keine andern Kräfte 
als diejenigen, unter 4p^en Einflüsse a als konstant betrachtet 
werden kann, so würde im folgenden Zeitdifferentiale die Glei- 
chung übergehen in: 

, dx j. , dy j, dx d'^x j. dy d^y j. 1 ^ . . 
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SO dafs nach der Regel von den partiellen Differentialen: 

ist. 

Die beiden ersten Glieder dieses Aggregats bedeuten die- 
jenigen Inkremente, welche nach dem Gesetze der Trägheit ein- 
treten; sie würden mit derselben Stärke auch dann stattfinden, 
wenn die auf den Punkt einwirkende Kraft ihre Thätigkeit plötz- 
lich einstellte oder wenn umgekehrt zu dieser Kraft noch eine 

neue hinzuträte. Im ersten Falle würden 37 und -=? konstant, 

dt dt 

da ja dann der bewegte Punkt eine gerade Linie (die Tangente) 

mit konstanter Geschwindigkeit beschreiben würde. — Die beiden 

anderen Glieder unserer Gleichung stellen denjenigen Zuwachs 

dar, welcher aus der fortgesetzten Wirkung der oben defihierteii 

Kraft entspringt. 

Gesellt sich nun aber in jenem ersten Zeitdifferentialte zu 

dieser Kraft noch eine neue, so erfahren die Gröfsen 

dx ., dy 

auch noch ein weiteres Inkrem'ent, das wir durch 

dt^' öezw. ^ 
bezeichnen wollen. 

Mit diesto Inkramewten wird aber die linke Seite der Glei- 
chuüg («) notwendig von Null verschieden, es sei denn, dafs 
diefilelben — wie etwa bleim Übergange von der linearen zur 
elliptischen Flächengeschwindigkeit, den wir bald betrachten 
werden — sich gegenseitig vernichten. Soll die Gleichung in 
der vorliegenden Form dessenungeachtet bestehen bleiben, 
so müssen wir der Konstanten a notwendig eine Veränderung 

beilegeii, wodurch jene neuen Zunahmien ausgeglichen werden, 
wenn wit nämHch mit a den in einer endlichen Zeit erreichten 
Zuwachs von a bezeichnen. Wir haben deshalb die unter diesen 
veränderten Bedinguhgen stattfindende Gleichung folgendermafsen 
anzuschreiben : 
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Die beiden Glieder der ersten geschleiften Klammer be- 
deuten auch hier das Inkrement der Trägheit, die der zweiten 
Klammer den Zuwachs infolge der Einwirkung der ursprünglich 
angenommenen Kraft, die beiden ersten Glieder der dritten 
Klammer stellen die Wirkung der neu hinzugekommenen Kraft 
dar und endlich das letzte Glied der dritten Klammer ist not- 
wendig hinzuzufügen, damit auch unter den neuen Verhältnissen 
die linke Seite der Gleichung noch der Null gleichgesetzt 
werden kann. 

Durch Subtraktion der Gleichungen (n) und'(w,) aber 
ergiebt sich: 

(^-|)-^-(*-|)-^-(l)S-o 0) 

und diese Gleichung bildet den analytischen Ausdruck des 
Princips. Sie lehrt uns die Geschwindigkeit 

da 

dt 
kennen, mit welcher infolge des Auftretens der neuen Kräfte 

dH . d^7] 

-d¥ ^^^ -W 

die Änderung a der Konstanten a sich vollzieht. Durch eine 
einmalige Integration würden wir diese Änderung a selbst er- 
halten — wie schon früher angedeutet worden ist Es bleibt 
nur noch übrig die Gleichung (1) zu integrieren. Dies wird 
erleichtert, wenn wir statt der Zeichen 

wieder die Formen -=7«- und -j^ einführen, uns dann aber stets 
dt^ dt^ ' 

erinnern, da& wir unter diesen Differentialquotienten lediglich 

die neu hinzugekommenen Kräfte zu verstehen haben. Die 

Gleichung (1) verwandelt sich dann in: 
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Die eingeklammerten Ausdrücke stellen offenbar die par- 
tiellen Differentialquotienten der Fundamentalgleichung (m) dar, 
bezw. nach den ersten Derivierten 

dbo dy 
dt' dt 
sowie der Eonstanten a. 

Um sie zu erhalten, hat man also nach und nach nur diese 
drei Grölsen als veränderlich, alle übrigen als konstant zu be- 
trachten. So würde z.B. aus der Gleichung: 



dy dx ^ 



folgen: 




so dafs die Gleichung (2) die Form annimmt: 
d^y d^x dx 
^ dt^ ^ dt^ " dt' 
In Worten wird demnach das Prinzip so lauten: 

Man suche zunächst für die Konstanten, deren Variation 
bestimmt werden soll, einen Ausdruck durch die Koordinaten 
und ihre Geschwindigkeiten 

dx dy 

^' »' ¥' i' 

Wir erhalten damit die Gleichung (m). Hierauf differenziere 
man letztere Gleichung lediglich in Bezug auf die Geschwindig- 
keiten und die Eonstanten. Die resultierende Differentialgleichung 
ergiebt dann unmittelbar die Änderungsgeschwindigkeit der 
Eonstanten. 

Hat man zwei Gleichungen, entsprechend zwei Eonstanten 
a und b: 



jy dx dy ^. 
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— ein Fall, der darin eintreten kann, wenn die Isolierung der 
Konstanten mit zu grofsen Schwierigkeiten verknüpft ist — so 
bleibt das Verfahren im wesentlichen bestehen. Nur erhält man 
für die Änderungsgeschwindigkeiten gleichfalls zwei Gleichungen, 
so dafs schliefslich noch eine Elimination erforderlich wird. — 
Um dieser auszuweichen, wird man deshalb wohl thun, jede 
Konstante womöglich abgesondert zu behandeln. 

Zuweilen bestehen zwischen den Konstanten einfache Be- 
ziehungen, wie zwischen dem Semi - Parameter j), der mittleren 
Entfernung a und der Excentricität € die Gleichung 

p «= a(l — £2). 

Wenn man in diesem Falle fetwa die Variationen von p und a 

bereits gefunden hat, so läßt sich selbstverständlich die Variation 

von € unmittelbar durch die ton p und a ausdrücken — wodurch 

dann eine besondere Behandlung von e entbehrlich wird. 

An dem Verfahren ändert sich nichts, wenn die Geschwindig- 

die du 
keiten 37' t? ^^^ ^^^r ursprünglichen Bahn selbst konstante Werte 

haben. Dieser Umstand würde vielmeäii* nur in den Gleichungen 
(w) und (»1) den Wegfall von Gliedern bedingen, welche in 
der Schlufsgleichutig (1) ohnehii verschtmnden sitid. Bd dei* 
demnächstigen Ableitung der elliptischen aus der geradlinigen 
Bahn werden wir einem solchen Falle näher treten. 



UI. Anwendung auf dasFlftchenprinzip derCentralbewegung. 

Wenden wir uns nun zu den Anwendungen des Variations- 
princips, so werde zunächst ein sehr einfacher Fall betrachtet, 
nämlich die Änderung der Flächengesbhwindigkeit beim Über- 
gange der oskulierenden geradHnigeti ifa die elKptische Bewegung. 
Bekanntlich versteht man unter Flächengeschwindigkeit die von 
dem Radiusvektor in der Zeiteinheit durchlaufene Fläche, und 
nach den Kepler sehen Gesetzen ist dieselbe konstant. Das 
wäre also nun nachzuweisen. 

Wir nehmen zu dem Ende an, der Planet stehe bei Beginn 
der Zeit t im Punkte P (s. Fig. 1) und die Kraft der im Ur- 
sprünge befindlichen Sonne höre in diesem Augenblicke plötz- 
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lieh auf zu wirken. Es würde alsdann der Planet in geradliniger 
Bahn mit konstanter Geschwindigkeit 

PPi ^ P,P2 - P2P3 U.S.f. 

fortschreiten. Aber auch seine, durch die Dreiecksflächen 

OPP^, OP^P^ u s.w. 
dargestellten Flächengeschwindigkeiten wären gleich, d.h. 
konstant, da ja die Dreiecke sämtlich gleiche Grundlinie und 
Höhe haben — eine Eigenschaft der geradlinigen Bewegung, die 
iibrigens auch sofort auf analytischem Wege erkannt wird. 
Denn da bei dieser Bewegung die Koordinateninkremente dx und 
dy konstant sind, so ist offenbar das Differential des vom Radius- 
vektor im Zeitelemente beschriebenen unendlich kleinen Doppel- 
dreiecks xdy—ydx, 
nämlich dxdy-'dydx «= 0. 
Es sind mithin die in jedem einzelnen Zeitelemente, demnach 
auch in gleichen endlichen Zeiten vom Vektor überstrichenen 
Flächen einander gleich. 

Für die geradlinige Bewegung ist folglich die doppelte 
Flächengeschwindigkeit c eine Konstante. Drücken wir dieselbe 
— den Forderungen unseres Prinzips gemäfs — durch die Koor- 
dinaten und Geschwindigkeiten aus, so haben wir hier un- 
mittelbar: 

dy dx 

Diese Gleichung entspricht also dem Falle, wo im Ursprünge 
keine Anziehungskraft thätig ist. 

Nehmen wir nun aber an, im Ursprünge beginne die 
Attraktion der Sonne auf den Planeten und umgekehrt wiederum 
zu wirken, so müssen wir — infolge des Auftretens dieser neuen 
Kraft — die bisherige Konstante c zunächst als variabel be- 
trachten. Indem wir den Zuwachs derselben während einer 
unbestimmten endlichen Zeit durch y bezeichnen, erhalten wir 
für die Geschwindigkeit dieser Zunahme nach dem Variations- 
prinzip: 

dy _ d^y d^x 

Es sind aber nach dem Obigen die neu auftretenden Kräfte 

Israel-Holtzwart. 2 
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dP r^ * 



mithin 



dy _ f{\^m)yx . fi\+m)yx 



dt 



+-'^"S- - -^ (0 



Die Geschwindigkeit -~ und also auch das Inkrement 

dt 
ist demnach zunächst indem der Oskulation folgenden Zeitelemente, 
dann aber auch — weil die Gleichung (l) ganz allgemein in 
allen späteren Zeitelementen besteht — für alle Zukunft = Null, 
d. h. die Flächengeschwindigkeit ist auch bei gegenseitige!: An- 
ziehung von Planet und Sonne konstant, nämlich glöiöh Aßc 
Flächengeschwindigkeit des der Einwirkung der Sonne ent2oge]ieü 
in der Tangente der Bahn sich fortbewegenden Planeten. 

Übrigens ergiebt sich aus der vorstehenden Betrachtung 
noch weiter, dafs die Beständigkeit der Flächengeschwindigkeit 
keineswegs an das Gravitationsgesetz gebund^i ist. Vieltnehr 
wird dieselbe fortdauern, wenn auch die Stärke der Anziehung 
eine beliebig andere Funktion des Radiusvektors wäre, wenn man 
z. B. hätte: 

dt''- ^> r 

dt^ ^' r 
wo durch F eine willkürliche Fanktion angedeatet wird. Dean 
auch in diesem Falle würde offenbar für alle Zeiten 

%-o 

sein. Wir haben demnach mit Hilfe unseres Prinzips (also ohne 
Anwendung einer eigentlichen Integration) detx ersten Funda- 
mentalsatz nachgewiesen: Wenn Planet und Sonne sich 
nach irgend einem von ihremwechelseitigen Abstände, 
also nach dem Radiusrektor r, abhängigen Gesetze 
anziehen, ist die Flächengeschwindigkeit der Pia- 
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netenbewegung konstant, Oder in Form einer oft«an* 
gewendeten Gleichung: 

xdy — xdx -« r*d(// = ccU (s), 

wenn c die doppelte konstante Flächengeschwindigkeit, also cdt 
den doppelten im Zeitelement beschriebenen Elementarsektor 
bedeutet. 

y. Entwickelung der elliptischen Bahn aus der osku- 
lierenden geradlinigen Bewegung, 

Denken wir uns das Koordinatencentrum (Fig. 1) zuerst 
wieder frei von anziehender Masse, so ist, wie schon zu Anfang 
unserer Untersuchung erwähnt, die Polargleichung der tangen- 
tialen Bahn: 

^ g 

sintp — mco8\fß' 
während die Konstanten m und q von den oskulierenden Koordinaten 

und Geschwindigkeiten (^o>yo> "-jf'-Jr) <i^rch die Relationen 
abhängen: 

da wir bereits wissen, daDs der Klammerausdruck eine Konstante 
und zwar die Flächengeschwindigkeit bedeutet, gleichviel ob jias 
Koordinatencentrum Sitz einer anziehenden Kraft ist oder nicht. 
Am elegantesten gestaltet sich nun — wie die Versuche 
lehren « die Konstantenvariation, wenn wir der Polargleichung 
die Form geben: 






1 



r «— 
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sin 
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Alsdann wird: 



wo -=7 und -^ die Geschwindigkeiten der Koordinaten in einem 

beliebig gewählten Oskulationspunkte der Bahn bezeichne» und 

2* 
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Konstante darstellen, wenn der Eoordinatenmittelpunkt ohne 
Anziehung, die Bahn also eine Gerade ist. 

Es beginne in jenem Punkte nun vneder die Anziehung 
der Sonne zu vnrken. Nennen wir 

A die alsdann eintretende Änderung von — 

Ai die Änderung von — > 

so hat man nach dem Variationsprinzipe: '\a 

dA ^ l d-x ^ f^x+m) X 

dt c dt- c r^ 

dAj 1 d^ _ , /(i-hm) y 

dt e dt' c r^ 

Mithin ist: 

Da aber, wie bereits erwiesen: 

cdt =-" r- dtp, also dt => —r- dip 

X ^ rcos ip 

y •« r sin ip, so folgt 

A -= ^/cos i//d^ « ^^J'-^sini/z + iT 

A, « ^^'J^/sin xpdip--- ^f-^ cos i/; 4-Z,. 

Damit ergiebt sich aber für die gesuchte Gleichung der Bahn: 
1^_ 

|- + AJ sin V^ - (^- -^ A J cos ip 

/i -^K-i-^l'P sin V; jsin ip -{^ ^K,-^-^^^^ cos ip\ cos xp 

Dieselbe unterscheidet sich — dem Wesen unserer Methode ent- 
sprechend ^- nur dadurch von der Gleichung der Tangente, dafs 
die Konstanten 
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1 , m 
— und — 

noch einen veränderlichen Zuwachs erhalten haben. 

Ordnet man nach Konstanten und Veränderlichen, so er- 
hält man: 

1 



/(!+' 



^ + (- -f-H:) sin rp-h-^K^] cos ip 



Fafst man endlich statt der Integrationskonstanten j? und K^ 
zwei neue Arbiträre G und k ■ vermöge der Gleichungen ein : 

^-^K = —CsinX 

a 

- + £,= -CcosA 

und dividiert aufserdem Zähler und Nenner mit - -"*"^"-, so geht 
die Bahngleichung aber in: 

^ - -- ^— («)• 

1 -f-TT co8((/;-f-A) 

In dieser Gestalt aber erkennt man sofort die Gleichung eines 
Kegelschnitts, dessen 

Semiparameter p = ^ 

/(l + m) 

Exzentrizität f. ^ G.p 
wahre Anomalie q> =^ \p-{-ly 
so dafs also die gewählte Achse der Abscissen mit dem Radius 
des Perihels (der grofsen Achse) den Winkel X einschliefst. 

VI. Bestimmung der Konstanten. 

Es liegt uns nun noch ob, die Konstanten G und X zu be- 
stimmen, und da dieselben von den ursprünglichen Konstanten 
K und Kl abhängen, so müssen wir zunächst die Werte der 
letzteren ermitteln. 

Betrachten wir die Gleichung (a) des Kegelschnitts, so 
unterliegt dieselbe — wie aus der gesamten Entwickelung her- 
vorgeht — der Bedingung, dafs sie in die Gleichung der Tangente 
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1 



r « 



1 . . w 

— sm V' cos?// 



für den beliebig gewählten Osknlationszeitpunkt übergehen mufs. 
Nennen wir also \p^ 

den Wert des Winkels yp (s. Fig. 1) im Momente der Oskulation, 
80 ist — zur Erfüllung dieser Bedingung — 

C C Y^ 

zu setzen. Hieraus folgt ferner: 

wenn man bei der Beduktion die osknUerende Bahngeschvindi^eit 

du f \ df / \ df ; ' 
ferner /"(i+m) 11 dx„ m dp, 



- v„ = -"• 







c^ jP ö' cd^ g' cdt 

und ciy. öte« . . 

Da die Exzentrizität der Bahn nach Obigem = C'p, 
so hat man: . 

r /(l+m) r^ 

Es ist demnach die Bahn: 

V ^ 2 

eine Parabel (€=1), wenn y = 0, 

/(l+m) r„ 

V ^ 2 

eine Hyperbel (« > 1), wenn >. -" > 0, 

V ^ 2 

eine Ellipse (e < 1), wenn ^ < 0. 

/ll+m) »"o 

Was ferner den Winkel l betriift, von welchem die Lage 
des Perihels (des Anfangspunkts der wahren Anomalie ?// + A = y) 
abhängt, so ergiebt aus dem Vorhergehenden sofort 

tg^-£ — (t). 

2 
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1 m 

Weil —, —^ K und K^ bekannt sind, so ist auch die Lage des 

Perihels, oder vielmehr die wahre Anomalie (i// + A) aus den 
oskulierenden Koordinaten und ihren Geschwindigkeiten, nämlich 

. ^'' ^«' dt dt 
ohne Weiteres bestimmbar. 

YII. IMe anderen Bahnclementc. 

Obgleich die vorstehenden Relationen sich zunächst auf die 
Koordinaten und Geschwindigkeiten im Zeitpunkte der Oskulation 
berieben, so haben dieselben doch allgemeine Giltigkeit, 
da ja dieser Punkt in der zu bestimmenden Bahn ganz beliebig 
gewählt werden kann. Die Gleichungen 

/(i+i«)l> •= c' (1) 

.2^1 + /^-?? (2) 

drückian desba,lb ganss allgemeine Eigenschaften der in Kegel- 
schnitten erfolgenden Planetenbewegung aus. 

Da ferner in der Elllipse der Semiparameter 
p « a(l — «'), so hat man auch: 

l-i-^ (3) 

jene merkwürdige Gleichung, welche lehrt, dafs die Hauptachse 
2a der Ellipse nur von r und t?, nicht aber von der Richtung 
der Geschwindigkeit v abhängig ist. 

Weiter ergiebt sich aus Gleichung (1) 

/(i+m) ap « /(i+m) &2 = ac* « a • —jf-r-' 

wenn U die Umlaufszeit des Planeten und b die kleine Achse 
bezeichnet. Denn da c die doppelte in der Zeiteinheit durch- 
laufene Fläche, abn aber die Ellipsenfläche darstellt, so ist 

2ahn 

'--TT' 

Man hat demnach: 
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Nennt man also den in der Zeiteinheit durchlaufenen Winkel 
(die 8. g. mittlereBewegung) n, so folgt die bemerkenswerte 
Gleichung: 

fii+m) « k\i^m) •= n^a^ (4). 

Dieselbe Relation in der ursprünglichen Form lautet 
Für eine zweite Planetenbahn ist mithin: 
so dafs ? 7^(i-{-m) g^ 

eine Beziehung, welche als der Ausdruck des verbesserten 
dritten Keplerschen Gesetzes zu betrachten ist. In der 
ursprünglichen (weniger korrekten) Form lautet dieses Gesetz 
bekanntlich: 

Die Eeplersche Gleichung. 
Es bleibt nun insbesondere noch übrig, die wahre Anomalie 

als Funktion der seit dem Periheldurchgange (wo 9>«0) ver- 
flossenen Zeit t darzustellen. Zu dem Zwecke müssen wir auf 
die Gleichung 

r^dg) ^ cdt ^ 

zurückgreifen und r^ durch die Polar gleichung eliminieren, wo- 
durch sich ergiebt: 

(l+ficosyP 
oder, wenn man vermöge einer Verbindung der Gleichungen (1) 
und (4), c durch n ersetzt und die Funktionen der halben Ano- 
malie einfuhrt: 

1 d<p 



ndt = -* 



(l-.)l(l+e)* cos^ 






Wird, zur Erleichterung der Integration, ein durch die Gleichung 
. ^u 1 — € . 2y ,. . 

*g-2=m*8 2 (5) 

definierter, exzentrische Anomalie genannter Hilfswinkel 

u 
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ZU Hilfe genommen, so erhält man nach einigen leichten 
Reduktionen 

ndt -= (1 — € cos u)du, 
woraus durch Integration die s. g. Keplersche Gleichung folgt: 

nt «« u — € sin w (6). 

Man hat so zunächst u, durch Vermittelung der Gleichung (5) 
aber auch g> als Funktion der Perihelzeit t — und umgekehrt 
t, wenn ü, bezw. q> gegeben ist. 
Wir kennen nun 

die mittlere Entfernung -« a 
die Exzentrizität =« e 
den durch die Gleichung (i) des Abschnitts VI berechneten 

Winkel X, 
drei Gröfsen, durch welche die Dimension und die Lage der 
Ellipse gegen die Koordinatenachsen bestimmt sind, wenn man 
bedenkt, dafs der eine Brennpunkt im Ursprünge derselben liegt. 
Aus dem gegebenen, zur Zeit der Oskulation stattfindenden 
Winkel xp^ (s. Fig. 1) 

in Verbindung mit l wird erhalten die wahre Anomalie im Zeit- 
punkte der Oskulation 

sodann aus ^^ die excentrische Anomalie «„ und hieraus die 
mittlere Anomalie 



-*>/^ 



nL'^hi/^^-^'t^ 



— alles zur Zeit der Oskulation. 
Man nennt die Gröfsen 

a, e, A, nt^ oder auch 9)^ 
die Elemente der elliptischen Bewegung. Sobald sie bekannt 
sind, lassen sich die Polarkoordinaten 

r und (f 
als reine Funktionen der seit dem Periheldurchgange verflossenen 
Zeit angeben. 

Eliminiert man nämlich mit Hilfe der Gleichung (5) die 
wahre Anomalie aus der Polargleichung der Ellipse, so erhält man: 

r == a(l — ccosw), 
oder a— r 



u = arc cos • 



ae 



26 A. Theorie der eUiptiachen Bewegung. 

Durch Sul^stitution dieses W^rt^s ia die Keplerache Gleichung (6) 
folgt: 

a — r ^/2a/r — r^ — ü 
arc cos ^ 



«(<+<o-=*/^^#+0-(7). 



€b8 a 

Entnimmt man andererseits der Gleichung (5): 

w«2arctg[^j^.. tg|- 

und substituiert in die Keplersche Gleichung, so wird: 

2arctg[^i^tg| -«8iB{?arotg[^[^-^g| }«n(Wo)-{8). 

In den Gleichungen (7) und (8) erscheinen die Polarkoordinaten 

r und (p 
als implicite Funktionen der ßeit der Oskulation verflossenen 
Zeit t oder auch der PerihekQit (^+^o)« 

YIII. Konstruktion der BahncUipse. 

Die relative Bahn eines Planeten um die Sonne, auch wenn 
ihre elliptische Natur und die Eonstanz der Flächengeschwindig- 
keit erkannt ist, kann doch nicht ohxie Weiteres aus den Koor- 
dinaten und Geschwindigkeiten 

dx dy 

eines Punktes konstruiert werden. Denn dieselbe ist, unter 
übrigens gleichen Umständen, auch eine Funktion der Planeten- 
masse m und des Erfahrungsfaktors /". Die Konstruktion, wenn 
man es überhaupt so nennen will, soll hier auch nur deshalb 
kurz berührt werden, um durch eine einfache geometrische Be- 
trachtung zu zeigen, inwiefern die Exzentrizität und Lage 
des Perihels von der Bichtung der momentanen Planeten- 
bewegung beeinflufst ist, während die mittlere Entfernung, 
wie wir gesehen haben, von dieser Richtung unabhäogi^ ist. 

Nimmt man den mittleren Sonnentag zur Einheit der Zeit> 
die mittlere Entfernung der JErde von der Sonne zur Einheit 
der Längen und die So^nenmasse zur Einheit der Massen, so 
hat man nach der an früherer Stelle begründeten Gleichung 
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da die dderischa Umlaufszeit der Erde 

' U « 365,256384 mittl Tagen, 
f(i4-«) - (0,0172021)2. 
Weil die Masse m der Erde in Teilen der Sonnenmasse «0,00000282^ 
so ist auch der Wert der Krafteinheit 

/•« t2« (0,0172021)^ 
also die G aufs sehe Konstante 

Je -= 0,0172021 
zu setzen, wenigstens wenn man sich mit einer Genaunigkeit von 
7 Stellen begnügt. Wählen wir diese Krafteinheit und nehmen 
ferner an, die grofee Achse sei vermöge der Gleichung 

12 v^ 

bestimmt worden, in welcher r und v unmittelbar durch 



x^ 



dx. 



dt 



dyo : dx^ 



gegeben sind, so ist die Ellipse leicht konstruierbar. 

Denn aus x^, y^ 
wird der Ort P(8. Fig. 4) 
des Planeten und dessen 
Eadiusvektor r geftmden, 
während aus 

dt • dt 
die Richtung NM der 
tangentialen Bewegung 
sich ergiebt Da nun der 
zweite Vektor q mit der 
Tangente bekanntlich 
denselben Winkel bildet und die Summe r4-^ der Leitstrahleu 
= 2a, so ist auch der Vektor q nach Lage und Gröfse, mithin 
der zweite Brennpunkt F gegeben. Wir kennen nun die Lage 
und Gröfse der grofsen Achse a sowie die doppelte Excentri- 
cität 2e, a.so die zur Konstruktion erforderlichen Elemente. Das 
Perihel liege in diesem Falle in JT. 
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Ändert sich die Bewegungsrichtung des Planeten — 
während Ort und Geschwindigkeit fortbestehen — geht beispiels- 
weise NM in NiM^ über, so werden auch die Winkel der Leit- 
strahlen von den vorigen verschieden ausfallen. Der zweite 
Brennpunkt rückt etwa nach JP^; die der Gröfse nach unver- 
änderte Apsidenlinie (grofse Achse) nimmt die Lage von A^S 
an. Es verwandelt sich somit die Exzentrizität FS in F^S und 
das Perihel geht von Fl nach /Ij. Die Bewegungsrichtung hat 
demnach zwar keinen Einflufe auf die Gröfse der mittleren Ent- 
fernung und folglich auch nicht auf die mittlere ümlaufszeit, 
wohl aber auf die Exzentrizität und den Ort der Apsiden, des 
Perihels und Aphels. 

IX. Zahlenbeispiel: Bdlm der Geres. 

Um die vorhergehenden Entwickelungen , insbesondere die 
Bestimmung der Konstanten, durch ein Beispiel zu erläutern, 
wollen wir 

die Bahn der Ceres 

behandeln und zwar auf Grund ihrer, am 30. März 1810, abends 
S** Qöttinger Zeit, oskulierenden Koordinaten und Geschwindig- 
keiten, nämlich: 

dx^ 



x,^-^ 1,732505 
y, « 2,064719 



J^ « —0,0086133 
^«- « +0,0062012. 



Da 



1. Charakter der Bahn, 



logr.= 0,43061 
log Ä = 8,23558 — 10, 
so nimmt der Ausdruck 

^— - den Wert -0,36137 

an. Er ist also negativ und die Bahn eine Ellipse. 
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2. FlächengeschwindigJceit 
Die doppelte Fläcliengesch windigkeit ist: 

3. Parameter. 
Für den Semiparameter folgt hiernach: 
l> = |y= 2,75023. 

4. Die Integrationskonstanten, 
Es ergeben sich die Werte: 

jj:« —^-.^ = —0,27854 
c^ r^ 



und da: 



■Ki=+S--~ = +0,23372, 



1 fLr 

- = — ^ - +0,30193 
q cdt ' 

^=.-%« -0,21738, 

q cdt ' 



so wird: 



log C « log r// - + k\\ h + kX^ 8,45534— 10 

logtgA = log| 0,15578 

i 

. f 55» 3' 45" 

'^ "°l 235« 3' 45". 
Ob die Apsidenlinie den einen oder den anderen Winkel mit 
der Achse der x bildet, kann erst bei Bestimmung der wahren 
Anomalie entschieden werden. 

5. Haupttiehse. 
Aus: 12 y/ 

folgt log a = 0,44205. 
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6. Exzentrizität, 
Die Exzentrizität hat den Wert 

£ « C.p ^ 0,07847. 

7. Wahre Anomalie, 
Man erhält: 

Vo- V^« + A == 285« 3' 45". 

Aus sin ipQ «= — folgt nämlich tp^ = 50^,. Da nun ferner nach 

der Polargleichung 

recos(p ^ P^r ^ 2,7502—2,6953 = 0,0549, 
so ist im vorliegenden Falle re cos y, demnach auch cos (p positiv. 

Es kann mithin 

y^ nicht« 500 +55« 3' 45", 

sondern mufe = 50« + 235« 3' 45" « 285« 3' 45" 

gesetzt werden. Damit sind nun aher die vier Elemente der 
Bahn bestimmt, nämUch: 

a, fi, X und y^. 
Durch a und « kennt man die Dimensionen der Bahn, durch 
den Winkel X ihre Lage gegen die gewählte Achse der x, (also 
die Lage des Perihels), durch (p^ den Ort der Ceres im Zeit- 
punkte der Oskulation. 

Wünscht man auch noch die mittlere Bewegung n und die 
Zeit t des Periheldurchgangs zu wissen, so erhält man die erstere 

durch die Gleichung: 

n^a^ «= k^ 

und, nach Berechnung der excentrischen Anomalie u aus der 
Definitionsgleichung 

die Perihelzeit t mit Hilfe der Keplerschen Gleichung: 
nt « u — e sinw. 

X. Erster IN achtrag: Die GescJwnndigkeiten der Koordinaten. 

Die Geschwindigkeiten der Koordinaten finden im Gebiete 
der theoretischen Astronomie und ganz besonders in den vorher- 
gehenden Entwickelungen eine so ausgedehnte Verwendung, dafs 
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es wohl angezeigt erscheint^ die wichtigsten hierher gehörigen 
Beziehungen noch kurz zu erwähnen. 

1. Geschmndigkeiten der Polarkoordinaten r und tp. 
Nach dem Früheren ist die doppelte Flächengeschwindigkeit 



.'■'i-i^/F 



dt 

(wenn man die Masse des Planeten neben der Sonnenmasse ver- 
nachlässigt, hingegen « Tc^l-^m-^/p^ wenn dieselbe berücksich- 
tigt wird), so dafB man sofort für die €resdiirlndigkeit der 
wahren Anomalie in einem Punkte der Bahn, dessen Vektor 
= r, erhält: ^ 

dtp k^p .. 

dt r^ ^^^• 

Ferner folgt aus der Polar^eichung r(l-|-€cos y)«p die (Je- 
schwindigkelt des Vektors " 

dr er^ . dw 

dt p ^ dt 

Da dcp dip kj/p 

dt^ dt ^ ~r^' 
so verwandelt sich die Gleichung in: 

dr __ kesitKp 

dt y/^ 



(11), 



wodurch also die Geschwindigkeit des Vektors als Funktion der 
wahren Anomalie ausgedrückt ist. 



2. Geschwindigkeiten der Parallelkoordinaten 

X '^ rcoQip 
y — rsin ^ 



Aus a; — r cos T^ 



ergiebt sich: 

dx , dip ^ dr kl/p . , ke , . ,„,, 

.. =- — rsmip-^ -i-cos ^^^"^ '^^ sm t//4---Tr- cos i//smy .(III) 

-| « _^cosi// -^4- sm «// ^ -» H — *-^cos xp +7= sm ^ sm y. .(IV) 

Die Gleichungen (I) bis (IV), deren Erweiterung auf ein drei- 
achsiges System keinen Schwierigkeiten unterliegt, bilden die 



32 A. Theorie der eUipÜBchen Bewegung. 

Grundlage für die Rechnung der Geschwindigkeiten und wurden 
auch bei dem obigen Zahlenbeispiele gebraucht, um aus den 
Elementen der Ceresbahn die Data der Rechnung abzuleiten, 
nachdem der Winkel ip willkürlich « 50® angenommen worden 
war. — Dieselben lassen sich mannigfach kombinieren, z.B.: 
,dx ^ , dy he . 

XL Zweiter Kachtrs^. 

Statt die Gleichung der elliptischen Bahn an die Polar- 
gleichung der oskulierenden Geraden anzuschlie&en — wie dies 
im Vorhergehenden geschehen ist — kann dieselbe auch aus 
der Gleichung in Parallelkoordinaten entwickelt werden, was in 
mancher Beziehung vorzuziehen ist. 

Für die Gleichung der Bahntangente hat man nämlich: 

dt^ dt^'^^' 

wo--lr-> -^ die oskulierenden, zunächst als konstant zu be- 
dt dt 

trachtenden Geschwindigkeiten bedeuten. Es ist dies nichts 
anderes als die Gleichung für die doppelte (erwiesenermalsen 
auch in der elliptischen Bewegung konstante) Flächengeschwindig- 
keit c und einerlei mit der durch Elimination von t sich ergeben- 
den Relation (vergl. Blatt 1): 

^ "" dxo dxo 

Läfst man nun die konstanten Geschwindigkeiten variieren, so 
folgt als Gleichung der neuen Bahn: 

Die Integrale können aber, wenn 

d^y _fa±myy 
dt^ ^ r» 

d^X ^ /(l-t-m) ^ 

dt^ r^ 

gesetzt wird, ohne Weiteres angegeben werdien. Man erhält 
dann zunächst: 







Zweiter und dritter Nachtrag. 


(dy„ , 
\dt ^ 


C 




läx^ 
\ dt 

/(l+mj 
C 


C 






K^ + 


/C1+m; 
C 


sini/^„, 


woraus 


durch die Substitution 










cosif/ 


/~^ 


X 



33 



sin i/* -h Z^i y ■= 



sin (/> 






. %.+ /il±2l8inVo=Cr 
gefunden wird: 

also ebenfalls die Gleichung einer Kurve zweiter Ordnung, 
deren Diskussion genau nach der früher gegebenen Anleitung 
gefiihrt wird. 

XII. Dritter Nachtrag: Ätisdruck der beschleunigenden Kräfte 
in Polarkoordinaten; Gravitationsgesetz. 

In der vorhergehenden Untersuchung wurden die beschleu- 
nigenden Kräfte durchgängig auf ein Parallelkoordinatensystem 
bezogen und demgemälis durch 

d^x d^y d^z 
dt^ ' dt^ ' df" 
ausgedrückt, von denen allerdings in dem bisher allein betrach- 
teten Falle einer ebenen Bewegung die dritte Komponente ver- 
mieden werden konnte. 

Häufig aber — sowohl in der Theorie der elliptischen wie 
der gestörten Planetenbewegung — wird verlangt, die äufseren 
den Planeten angreifenden Kräfte nach drei zu einander senk- 
rechten Richtungen zu zerlegen, deren eine mit dem jedesmaligen 
Radiusvektor zusammenfällt, deren zweite senkrecht zum Vektor 

Isiadl-Uoltzwari. ^ 
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in der (ev. momentanen) Ebene der Bewegung und deren dritte 
senkreoht zur Bewegungsebene steht, mit anderen Worten: es 
wird verlangt, die äufseren Kräfte allgemein durch Polarkoordi- 
naten und ihre Differentialquotienten auszudrücken. Die dritte 
Komponente, senkrecht zur Ebene der Bewegung, können wir 
auch hier vorläufig aufser Betracht lassen; im Falle des Bedarfs 
unterliegt die Herstellung ihrer allgemeinen Form kein^rld 
Schwierigkeiten. 

Bei diesen Entwickelungen ziehen wir es vor, den synthe- 
tischen Weg einzuschlagen, da derselbe mehr als der analytische 
geeignet ist, die Bedeutung der allgemeinen Ausdrücke zur An- 
schauung zu bringen; von denen wir schon an dieser Stelle eine 
sehr wichtige Anwendung machen werden. Im Vorhergehenden 
haben wir nämlich auf Grund des Gravitationsprinzips die Natur 
der Bahn entwickelt. Die in Rede stehenden Ausdrücke werden 
uns aber in den Stand setzen, auch die umgekehrte Aufg(ibe 
zu lösen, nämlich aus der Natur der Bewegung (also aus den 
Kepler sehen Gesetzen) die Eigenschaften der Kraft (das Gravi- 
tationsgesetz) und zwar auf dem denkbar einfachsten Wege 
herzuleiten. 

EntwicJcehmg der Differentialgleichungen der Kräfte, 

P1P2 (Fig. 5.) sei das Bahnelement des Planeten im ersten 
Zeitdifferentiale. Macht man i^^a"^ -Pi^2> so würde sich der- 
selbe nach dem Gesetze der Trägheit im Punkte IZj befinden, 
wenn er im zweiten Zeitteüohen lediglich sich selbst überlassen 
wäre, unter der Einwirkung irgend welcher Kräfte gelange er 
jedoch nach P3. Schlägt man mit den Radien SP^ «= r und 
SP2=r-hdr die Kreisbogen P^N und P2C, zieht man ferner 
PjD parallel Pt-M, so ergeben sich — vorausgesetzt, dafs Bogen 
MN= MPi gemacht wurde — leicht die folgenden Beziehungen: 
P^C-P^M^dh- 
^P^SN^^P^SP^-P^SP.^d^^p. 
Zum Unterschiede hiervon setze man 

n^J-P^M^A'^r 

^n.sp^-p^SP.^A^tp, 

indem mit r und (f die Polarkoordinaten bezeichaet werden. 



Entwickelung der Differentialgleichungen der Kräfte. 
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Die äufsere Kraft Kr, welche im Radiusvektor den 
Planeten Pj angreift, ist nun offenbar gleich der im Radiusvektor 
wirkenden Totalkraft ^ 

dV 

vermindert um die aus der Trägheit hervorgehende Centrifugal- 
kraft A'-V 

demnach: ^ d'^r AV 

In ähnlicher Weise ergiebt 
sich für die äufsere Kraft K^, 
welche senkrecht zum Radius* 
Vektor in der Ebene der Be- 
wegung wirkt, die Gleichung: 

^*^^d^-^"d^ • •^^>- 
Es kommt nun darauf an, die 
GröJaen 

AV , AV 

der Gleichungen (1) und (2) 
durch die eigentlichen Diffe- 
rentialquotienten . der Polar- 
koordinaten zu ersetzen. 
Offenbar ist 

=- JF-^dr 
(bis auf Gröfsen höherer Ord- 
nung; denn -^ilaPjF ist nach 
Konstruktion « PjPtJf und 
JC.DP2J- d(p, also 

n,± - '^»^2- "sin IJa^" 




P^^i cosdy 



- dr). 



sinPzPiüf 

Ferner hat man gleichfalls bis auf GrÖfsen höherer Ordnung: 

JF -= P'iJ . d(p «= rd(f . dcp = rdtp'^, 
so dafs schliefslich : 
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Da aber auch 

so folgt: 

AV ^ rd(^^. 

Damit ergiebt sich für die (äufsere) Kraft im Vektor: 

^^"w-'w (^«>- 

Andererseits ist augenscheinlich: 

APtSPj - P2Ä/I3, oder: 
r{r + dr)d(p « (r + dr) (r -i- 2dr + A V) (ö!<)p + A^y), 
woraus mit Vernachlässigung der Gröfsen höherer Ordnung folgt: 

Für die zweite Komponente erhält man also die Gleichung: 

Entwickdung des Gravitationsgesetzes, 

Für die Geschwindigkeiten 

dw , dr 
-jT und -,, 
dt dt 

haben wir an einer firüheren Stelle die Gleichungen gefunden: 
d(p _ ]^J^, ^'^ _ A;e sin y 

Dabei wurde die Richtigkeit des Gravitationsgesetzes voraus- 
gesetzt. 

Jetzt soll gerade umgekehrt letzteres aas den erfahrungs- 
mäfsig festgestellten Kepler sehen Gesetzen deduziert werden. 

Wir haben deshalb jene für die Gleichungen (la) und (2a) 

unentbehrlichen Geschwindigkeitsausdrücke neu zn begründen. 

Jedenfalls ist auch hier (vergl. oben): 

"iahn 271 , — « ,— 

c = — jj— =» 'jj'(^Y^^ ^ na^yp* 

Nach dem dritten Gesetze Keplers ist aber 
nai «= Konstans -* k, also 
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wiB oben. Nur darf man nicht übersehen, dafs das dritte Ge- 
setz in der vorliegenden Form nicht völlig genau ist. Mit dieser 
Einschränkung haben wir demnach wieder: 



Die zweite Gleichung 



dt 

dr ke sin q) 



dt ^f 

ist aber nur eine Folge dieser ersten in Verbindung mit der 
Polargleichung; auch sie besteht demnach lediglich auf Grund 
der Kepler sehen Gesetze. 

Werden nun diese beiden Relationen nochmals differenziert, 
so ergiebt sich: 

d^r ke dw 

dt^ y/p ^ dt 

d'^if 9j,Vpdr 

-dF^~'^^^di' 

Durch Substitution dieser vier Ausdrücke in die Gleichungen 
(la) und (2a) erhält man sodann mit Zuziehung der Polar gleichung : 

K, - - *■ 

Unter der Voraussetzung, dafs der Planet sich in einer Kepler- 
schen Ellipse um die Sonne bewegt, reduziert sich also die 
äulsere ihn angreifende Kraft auf eine einzige im Radius- 
vektor, umgekehrt dem Quadrate der Sonnenentfer- 
nung wirkende Kraft. 

Die zweite Thesis des Gravitationsgesetzes, nach welcher 
die Wirkung der Schwere eine wechselseitige ist und proportio- 
nal den Massen vor sich geht, ist strenggenommen nur eine, von 
den Erfahrungen bestätigte Forderung der Logik. 

Wir glauben hiermit die Theorie der elliptischen Bewegung 
beschliefsen und auf der gewonnenen Grundlage dem zweiten 
Teile unserer Untersuchung, der Variation der Elemente 
einer elliptischen Planetenbahn, zuwenden zu können. 



B. 
Theorie der speziellen Störungen. 



XIII. Entwiekelnng der Eneke- Lagrangeschen Störungs- 
formeln. 

Nach unserer bisherigen Untersuchung erfolgt die relative 
Bewegung eines Planeten um die Sonne in einem Kegelschnitte, 
wenn lediglich die wechselseitige Anziehung dieser. zwei Körper 
in Betracht gezogen wird. Da indessen die Schwere eine all- 
gemeine Eigenschaft der Materie bildet, so werden jene zwei 
Körper auch mit allen übrigen Planeten in derselben Wechsel- 
beziehung stehen. Infolge dessen weicht jede Planetenbahn 
von der Idealgestalt einer Ellipse mehr oder minder ab, und 
diese Abweichungen sind es, welche man Störungen nennt. 
Je nachdem nun diese letzteren entweder als allgemeine 
Funktionen der Zeit t oder nur für begrenzte Zeiträume 
durch die mechanische Quadratur bestimmt werden, heifsen die- 
selben allgemeine, bezw. spezielle Störungen. So sehr auch 
die Überlegenheit der ersten Darstellungsart in theoretischer 
Beziehung anerkannt werden mufs — insofern allgemeine Theo- 
reme nur auf diesem Wege begründet werden können — ebenso 
sehr ist ihre praktische Verwendbarkeit als eine beschränkte zu 
bezeichnen. Dieselbe versagt oder wird wenigstens unsicher, 
sobald die Exzentrizitäten und Neigungen der Planetenbahnen 
eine nur mäfsige Gröfse überschreiten, da in diesem Falle die 
stets gebrauchten Reihen den für praktische Zwecke unentbehr- 
lichen Grad der Konvergenz verHeren. Man ist deshalb ge- 
zwungen, auf Methoden bedacht zu sein, welche dieser Ein- 
SQhränkung nicht anheimfallen, wenn sie auch der Allgemeinheit 
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ermangeln. Alle diese Veffahrungsarten, die, trotz ihres gemein- 
samen Grundgedankens, einer sehr mannigfaltigen Gestaltung 
fähig sindj faist man unter dem Namen 

Methoden der speeiellen Störungen 
2üsanmien. Nur auf diese wird die folgende Untersuchung sich 
einlassen, innerhalb dieser Grenzen aber den Gegenstand mög- 
lichst vielseitig behandeln > deshalb sowohl die Störungen der 
elliptischen Koordinaten als diejenige d«r Elemente in den 
Kreis der Betrachtungen ziehen. Auch hiei* werden wir uns 
vorwiegend auf das uns bereits geläufige Prinzip der Variation 
stützen. 

Wir beginnen mit einer vollständigen Ableitung der Encke- 
scheti Störnngsformeln, nicht blofs wegen ihrer immer noch 
fundamentalen Bedeutung für die gesamte Störungsrechnung^ 
sondern auch um dem Leser durch parallele Behandlung eines 
bestimmten Falls ein sicheres Urteil über die beiden Formen 
de« Variatiansprinzips — der von Encke benutzten mathe- 
matischen und der hier gebrauchten dynamischen Form — zu 
ermöglichen. 

Diese Störungsformeln erfordern eine nachträgliche Ent- 
wickelung der nach Richtung der Polarkoordinaten zerlegten 
störenden Kräfte, zu welchem Zwecke von Encke ein veränder- 
liches Parallelkoordinatensystem eingeführt wird. Da uns zu dieser, 
mit beträchtlichen Nebenrechnungen verknüpften Transformationen 
ein zwingender Grund nicht vorzuliegen schien, so trugen wir 
kein Bedenken, ein, wie wir glauben, einfacheres oder wenigstens 
unmittelbareres Verfahren an die Stelle zu setzen. 

Erst nach Erörterung der Enckeschen Gleichungen werden 
wir zu einer in mancher Beziehung neuen Behandlung unseres 
Gegenstandes übergehen. 

XIV. Allgemeine Bezeichnungen. 

Den auf Seite 5 zusammengestellten und durch Fig. 2 er- 
klärten Bezeichnungen sind noch folgende hinzuzufügen, welche 
teilweise auch schon früher gebraucht sind (s. Fig. 6): 

V = Geschwindigkeit des gestörten Planeten in der Bahn, 
u = exzentrische Anomalie desselben (definiert durch die 
Gleichung 5 auf Seite 24) 



bb 

•i-i 
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l = Länge des gestörten Planeten in der Bahn « dem 

Doppelwinkel FA-H SIP (s. Fig. 6), 
k = SIM '-=' Länge des gestörten Planeten in der Ekliptik, 
ß = Breite in Bezug auf die Ekliptik, 
Sl •= Länge des aufsteigenden Knotens Sl der Planetenbahn 
« F^ der Figur 6, 
PR « Verlängerung des gestörten Vektors r, 
PS «= Richtung in der gestörten Bahnebene, senkrecht 

zum Vektor r, 
P'^ =^ ^Senkrechte auf der gest. Bahnebene im Orte 

des Planeten, 
PQ *= Richtung der störenden Kraft, 
PWi =*- Richtung J_ r in der auf der Ekliptik senkrechten 
Ebene PMS. 
Fig. 6. Femer bedeuten: 

Ekltp^kpol QS. QR. QW und QW, die 

Winkel, welche die Richtung PQ 
bezw. mitP/S, PR, PTFundPTFt 
einschliefst. 
Nennt man also 

P die Totalkraft, welche die 
momentane elliptische Bahn des 
Planeten stört, 
so ist: 

1. die Komponente derselben in der Richtung 

des gestörten Vektors = P-cos^JB, 

2. die Komponente J_r in der gest. Bahnebene « P- cos QS, 

3. die Komponente _L der gest. Bahnebene -= P:COs QW, 

4. die Komponente _L r in der auf der Eklip- 
tik senkrechten Ebene, in welcher also 

die Breite ß gezählt wird = P-cos ^TF|. 

Aufserdem wird von Encke noch eine weitere Komponente 
P • cos QT eingeführt, welche in die Richtung der Tangente der 
gestörten Bahn fällt. Wir werden dieselbe stets — ebenso wie 
die vierte Komponente P - cos QWf — durch die drei zuerst ge- 
nannten Komponenten ersetzen, mit Hülfe der unten a. b. 0. 
entwickelten Gleichungen. 




Serhs Hilfsatze. 41 

Sechs Eilfsätze, 

dl dß dl 
I. Beziehung zwischen den Geschwindigkeiten -^» -^* -^^. 

Aus (8. Fig. 6): 
coscß -= 8in/?sin(/9 + d/^)4-co8/?co8(/y + d/9)co8dA 
ergiebt sich, wenn 

co%dl = 1 5 — !-••• 

cos a A — 1 — H 

gesetzt wird, nach gehöriger Reduktion: 

eine Relation, welche der Figur auch unmittelbar entnommen 
werden kann. 



II. Beziehung zwischen den Geschwindigkeiten ^ und ^• 

Die Gleichung (s. Fig. 6 und die Bezeichnungen): 
tg(i— a) «= eositg(Z— Q) 
führt, wenn i und Q als konstant betrachtet werden, durch 
Differenzieren auf die Relation: 

1 dl cosi dl 

cÖ82(X^^*5f QO%\l—^iydt' 
woraus, wegen 

co8(ü — Sl) «=• cos/yco8(^ — Q), 
erhalten wird: 

dl cos^/J d^ 



dt cosi d^ 



dt ■ 



III. Beziehung zwischen den Geschwindigkeiten tt und 
Durch Kombination von I und 11 ergiebt sich: 

= tg/?cotg(Z— a).;=-. 



Tt-'dty^ 



cm^ß öf 6v / j^ 
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IV. Beziehung zwischen den Komponenten P cos QW^ Pcos C? TT und P cos QS. 

Fig. 7. 

Ist © der Winkel zweier Geraden, 

welche mit drei aufeinander senk- 
rechten Achsen die Winkel «, ß, y^ 
^19 ß\y Yu bilden, so hat man nach 
der analjrtischen Geometrie: 
cos © « cos a cos «1 + cos /5? cos ft 

-l-cosycosyi. 
Hiemach ist (s. Fig. 7): 
cos QWt - cos ÖTrcos(90— iO 
+C0S ö/g8in(90— i,)+cosÖBcosTrtE 
— cos Q TFsin i| + cos QScos if , 
so dafs: 

V. Beziehung zwischen den Geschwindigkeiten — ^— und -^. 

Betrachtet man die in der elliptischen Bewegung konstante 
Länge Q des aufsteigenden Knotens N (s. Fig 7) als veränder- 
lich, indem JV nach N^ rücke, so ist : 

N^P-NP ^ N^P~TP ^ N,T cosiJVJVi, 

oder: 

d(l — Q) =• — cosi'dA, 
Fig, g wo das negative Vorzeichen ein- 

P tritt, weil eine Zunahme des s g. 

Arguments der Breite Z—ft einer 
Abnahme der Knotenlänge ft ent- 
spricht. 

Übrigens ergiebt sich diese 
_^^..^.. Gleichung auch leicht durch Diffe- 

EJcIip renzierung der Relation 

C0S(? — A) = C0S/?C08(A — ft) 

wenn man Si allein als veränderlich betrachtet, wie dies bei 
Entwickelungen nach der Variationsmethode geschehen mufs. 
Man erhält so: 

sin(A - A) d(X — A) = cos ß m^{X —Si)d{k-S^). 



P^ 
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Da offenbar d(A- ß) -= — du, 

so wird: 

jj/7 n\ c<>8 ß sin (^— ft) jn • jn 

d(? — ö?) = . — /y ^ o X ' dft = — oos t dft. 



VI. Beziehung zwischen der exzentrischen Anomalie u und der waliren Anomaliu rp. 
Aus den Polargleicbungen: 

r — a -= «ficosM • 

r — a -• — <iB^-^rBCfy^q> 

findet man darch einige leichte Umformungen: 

a yi> sin w ** r y/a sin y , 
wenn, wie bisher, mit p der Semiparameter =^ a(l — «*) be- 
zeichnet wird. 

XV. StSruug des Parameters und der mittleren Entfernung, 
Reduktion der Enekeschen Tangentialkraft. Störung der 
Exzentrizität, der Enotenlänge, der Neigung, des Perihels, 
des Winkels zwischen Perlliel und Knoten, der mittleren 

Anomalie« 

Durch die vorstehenden Hilfsätze sind wir nun hinlänglich 
vorbereitet, um ohne Unterbrechung die Störungsformeln der 
elliptischen Bahnelemente auf Grund der Variationsmethode ent- 
wickeln zu können. 

1. Störung des Param^ers. 
In der bekannten Grundformel der elliptischen Bewegung 



dl 

r* — 

dt 



«2 



-T.-^kyp 



erscheint, wie es das Variationsprinzip verlangt, die ellip- 
tische Konstante p bereits unmittelbar ausgedrückt durch Polar- 
koordinaten und ihre Geschwindigkeiten. Zwar haben wir im 
ersten Teile unserer Untersuchung stets die Parallelkoordi- 
naten und deren Geschwindigkeiten benutzt, allein vermöge der 
zwischen den beiden Arten von Koordinaten stattfindenden Be- 
ziehungen iet selbstverständlich eine Vertauschung beider jeder- 
zeit gestattet. 
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Diflferenziert man nun, nach Mafegabe der Väriations- 
methode, jene Gleichung in Bezug auf die Konstante p und die 

Geschwindigkeit -^, so folgt: 

, dH h dp 



P cos Q8, 



oder: ^^ 2ryfp 

dH 
da offenbar r-^7T"=der in die Bahnebene fallenden, zum Vektor r 

senkrechten Komponenten der störenden Kraft P ist. 

Es mag bei dieser Gelegenheit nochmals in Erinnerung ge- 
bracht werden, dafs hier, wie im folgenden, die Differential- 
quotienten der Geschwindigkeiten nicht deren Gesamt-Inkremente, 
sondern nur diejenigen Teile umfassen, welche den störenden 
Kräften entspringen — oder, wie wir uns früher ausgedrückt 
haben, den neu auftretenden Kräften (vergl. Teil A, Analyt. 
Ausdruck des Prinzips). 

2. Störung der mittlerenr Enifemung. 
Die an einer früheren Stelle gleichfalls begründete und 
schon wiederholt besprochene Gleichung: 

in die vorgeschriebene Form über: 



a « 1 



r k^\dt) Jc^\dtj 



woraus durch Differenzierung erhalten wird: 

dt k^ ' dt' dt' ■*■ h^' dt' dt^' 

Setzt man nach dem Obigen (Teil I): 

cU k\/p - dr ke 

jx = —^ und -^ = — ^=7 sm cp , 

dt r^ dt Jp ^ 

so ergiebt sich: ^ 

3T = 2 /^ • Pcos QS H 7^ sm w • Pcos QB. 

dt kr ^ Jc^p 



. Beduktion der Enckeschen Tangentialkraft. 4o 

Anmerkungr* BedvMion der Enckeschen Tanyentialkraß. 

Wie bereits erwähnt, hat Encke aufser drei aufeinander senkrechten 
Komponenten der Kraft P noch eine weitere, in die Bahntangente ^fallende 
Seitenkraft zu Hilfe genommen. Auf folgendem Wege kann dieselbe auf die 
anderen Komponenten zurückgeführt werden. 

Differenziert man die vorstehende Gleichung 



<-M:) 



unmittelbar nach der Konstanten a und der Geschwindigkeit t>, so kommt 

da 2ä* dv 

di — 'W^'^^ 

dv 
Offenbar stellt ^ die störende Kraft in der Tangente dar, weshalb 

da 2a« 

da 
Eine Vergleichung dieses Ausdrucks von ^ mit dem oben gefundenen liefert 

die gesuchte Beziehung: 

k V« kf^ sin a> _ 

Pcos QT « — ^- Pcose'S'+ —^ • Pcos QB, 



fr a 



wo V durch k\/— — — ersetzt werden kann. 

3. Störung der Exzentrizität 

Die analytisch -geometrische Belation 

p « a(l — «2) 
ergiebt : 

de ^ p da 1__ dp 

~dt '"2iäFBdi~ 2a€dt' 

so daßs nach Einführung der Werte von , und J\ folgt. 

j^ = ^^-^ /-^- Pcos QS 4- ^^^—^ — - • Pcos QB. 

dt arek ^ k 

4. Störung der Knotenlänge. 
Aus (vergl. Hilfsatz III): 

f = tg^cotg(i-a)|; 

findet man: 
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d'ß igl dl d{l-Sl) . , ^ ^ ,, ^dH 

oder (vergl. IV und V, sowie Fig. 7): 
sini|-Pco8 ^TT-l-cosii'Pcos QS « 

tg/?cos^Y £.^+eo8vPco8(g& 
Schliefslich : 

^^ Ä yi> sin i 

5. Störung der Neigung, 
Nach Hilfsatz III ist: 



Mithin 



oder: 



cos* ß. 

cos* ß, 



..dl dH . . ./dlVdi \dl dH dß d*ß 

'^''dt-w^-''^'''^'\dt)-dr[jä-dr^-dt-w 

der: 
cosUj^^PcosQS^rcostsimi^j^^ . ^^ 

« ^- J'cos QS -^LinH PcoHQW^cosi, Pcos Qs\cos^ß. 

Die verschiedenen Komponenten P oos QS heben sich bei näherer. 
Betrachtung ihrer Koeffizienten aJif, und man erhält nach einigen 
Reduktionen als Resultat: 

^^ "l^cosCi-a). Pcos OTT. 



dt 



k}/p 



6. Störung der wahren Anomalie g> ais Folge der Perihelstorung* 

Die Polargleichung: 

r(l -H« cos (p) ^ p 
liefert : 

d(p cosy de_ 1 dp 

dt"^ e^mfp dt re sin (p dt 
woraus durch Substitution: 



* Vergl. die nächstfolgende Anmerkung. 
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dif QOsg> 



dt € sin (f 



aehr ky^p J . 

r« Sin y k 
und durch Reduktion: 

7. Störung des Winkels lo ewischen Perihel und Knoten. 

Nach Fig. 2 ist: 

o) ^ 4P— (i — ft), 
so dafs 

dw d^ d(? - fl;) 

'di'^ dt dt ' 

oder, da hier von dem letzten Quotienten ausschliefslich die 
Änderung des Sl in Betracht kommt: 

d(ü Ir-hp) sin q> -n r^a 1>C08^ ^^ ^^ . .dSi 

-jT -= — - \^ - — 'PcosQS^^ — T^-Pcosö^ + cos*-^ 

WO fiir -^ der oben gefundene Wert zu setzen ist. 

8. Störung der miMleren Anomalie. 
Es sei 

Wo«" «4o- *offlö«<o 

die mittlere Anomalie und 

k 
Wo— — j 

Oo^ 

die mittlere tägliche Bewegung, beides zur Zeit der Oskulation 
(in der ,,Epoche"). 

Wären diese beiden Elemente — wie in der rein elliptischen 
Bewegung — unveränderlich, so bestände fiir die mittlere Ano- 
malie M zur Zeit t nach der Epoche die Gleichung: 

M ■« mQ-hn^t 
Nun ändert sich aber sowohl m^ als no, ersteres weil €q und Uq 
gestört werden, letzteres weil a^ ^ner Störung unterliegt. 
Nennen wir dn 

dt 

/^ V ' r^y THE 

/ ' . . .-. 1 T V • 

\ ^' . •^^ 
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die Geschwindigkeit der Störung von tig, so ist das Inkrement 
von fio ün ersten Zeitteilchen 

dn 



dt 



'dt 



Demnach beträgt die mittlere Bewegung zn irgend einer Zeit ti 



fio 4- 



j dt 



dt 



Vermöge derselben legt also der Planet während der Zeit t den 
Weg zurück 



^-'^-^f^fti'^' 



Wird ferner mit 



dm 
'dt 



die Geschwindigkeit der Störung von m^ bezeichnet, so hat man 
zur Zeit t nach der Epoche 



mo 



-/w* 



an Stelle von mo zu setzen. Die gestörte Anomalie zur Zeit t 
ist mithin: 



Jf - 



m. 



■r-. 



dm 
'"-^l-dt 



dt 



dn ,.1* 



^'+M^dt^' 







(x). 



Was zunächst 



dn 

dt 



betrifft, so hat man aus der Gleichung: 



n 



li 



durch Differenzierung: 

dn 'i k da 

di~ 2 a'i dt 

Vp 



3 



,L P cos Ö-S + ^^ P cos «ü 
'''ya yap 



(tf). 



* Yergl. die nächstfolgende Anmerkung. 
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um -^ zu bestimmen, differenzieren wir zunächst die Gleichung 

m = w — € sin w. 

Man erhält 

dm ,. .du . de 

-TT == (1 — € cos u)^r — smw^r 
dt ^ ^ dt dt 

r du . de ^ , . 

-ä-5^-"°"äJ ^«> 

Es ist aber nach der Definition 

i/rntg|-i/r:^.tg|, 

also du __ sinte dff a sin u de . . 

dt ^ sin ^ dt p dt ^^^' 

Wird -^ aus (a) und (/?), sowie sin u mit Anwendung des Hilfs- 
satzes VI eliminiert, so folgt: 

dm r^ d(p r(p-hr) . de 

-T- = ;=- -^ -.- - sm y^yT- 

»' ayop »r pyap ^^ 

Endlich ergiebt die Substitution der früher gefundenen Werte 

von -^ und ^ nach gehöriger Reduktion: 

dm 



dt k 



i IL ^P^p,o.QE^l±rsin<f.PoosQ8 



l/. 



•(^). 



Ersetzt man die Quotienten -, und ^ der Gleichung (x) durch 

ihre Werte in (y) und (g), so erhält man die Bestimmungsgleichung 
fiir die gestörte mittlere Anomalie M. 

Anmerl^ung !• Die vorstehenden Entwickelangen würden anschaulicher, 
wenn man die während einer unbestimmten Zeit t entstehenden Störungen 
durch besondere Zeichen hervorgehoben, also etwa den Stönmgsbetrag der 
Exzentrizität a durch 

mithin die Geschwindigkeit der Störung durch 

dt 

angedeutet hätte. Es ist dafür, um auch äufserlich möglichste Übereinstimmung 

ds 
mit den Encke schon Formeln zu erreichen, einfach ^r geschrieben worden. 

Israel-IIoUzwart. 4 
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Aufserdem erforderte der Gang der vorstehenden Entwickelungen , dafa 
auch die Störungen 

dp , dtp 

dt ""d s 

mit aufgenommen wurden , was indessen für die Form der anderen ßtötutags- 
gleichungen ohne Einflufs ist. 

Eine etwas merklichere Abweichung tritt dadurch ein, dafs, wie schon 
erwähnt, die von Encke giebrauchte Tangwitialkomponente überall vermieden 
wurde. Wird dieselbe init Hilfe der unter 2 mitgeteilten Gleichung wieder ein- 
geführt, so gehen unsere Formeln sämtlich in die von Encke über. Denn das 

umgekehrte Vorzeichen in der Störungsgleichung von -^ rührt nur daher, weil 
« bei Encke in umgekehrter Kichtung gewählt wird. 
Wichtiger erscheint ein andrer Differeilzpunkt. 

Nach Encke ist nämlich (vergl. Berl. Jahrb. 1837, S. 315) die gestörte 
Anomalie 

M ^ [m^-j- A»] + [«0 + A«] ^, ^öt 
( i 

hingegen nach unserer Gleichung: 

t t t 

M^\m,J^f^dt + \n,t+fdtf^dt 

L 0» J L 

anzunehmen. Wenn aber die mittlere Bewegung zu irgend einer Zeit t nach 

der Oskulation 

i 

««o + Aw««o+y ^äi 



ist, so kann der aus dieser mittleren Bewegung hervorgehende Teil der mitt- 
leren Anomalie M — wegen der Veränderlichkeit von /\n — nicht wohl 

-=(^o + A«X 
genommen werden. Man mufs vielmehr schliefsen: Während irgend eines Zeitr 
teilcbens dt wird durch die veränderliche mittlere Bewegung das inkrement 

I 
^dn 



»0+/; 



^^dt\.di 



erzeugt; im Verlaufe der Zeit £ entsteht also der Zuwachs: 



^o^+ß^ft 



Differenziert man, so ergiebt sich aus unserer Gleichung (wenn man nur auf 
die störenden Glieder Rücksicht nimmt): 
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dM dm Pdn 



während nach Encke (vergl. S. 316 a. a. 0.): 

dM dm ndn dn 

It-^+J-M^^-^^di 

dm 
folgt. Da aber andererseits nach Encke -^ — dem oben in {z) angegebenen 

dn . dM 

Werte minus ^^, so gelangen schliefslich die beiden Werte von —n- dennoch 

zur Deckung. 

XJnswer unmafsgeblichen Ansicht nach ist indessen das Erscheinen des 

dn dm 

Gliedes t ^ in der Gleichung für —^ ebenso wenig berechtigt, wie in der 

dM 
Formel für -rr. Dies giebt uns Veranlassung, zunächst nochmals auf die 

Störung der wahren Anomalie (mit der in No. 6 angenommenen Beschränkung) 
zurückzugehen. Denn man kann fragen, warum bei der Ableitung der Formel 
aus der Polargleichung 

Wl-f-f C08<^) «= |> 

allein *, p und c/- als veränderlich betrachtet wurden, währ^d r 2X& konstant 
angenommen worden ist. Es hat damit folgende Bewandnis: 

Die Gresamtstörung von (p besteht aus zwei Teilen, einer wahren und 
einer scheinbaren Störung. Die erstere beruht auf einer wirklichen Änderung 
der mittleren Bewegung des Planeten. Zar Einführung der scheinbaren Störung 
sieht man sich nur deshalb gezwungen, um den elliptischen Charakter der 
Bahn aufrecht zu erhalten ; sie läuft denn auch thatsächlich auf eine bloüse 
Verschiebung des Perihels hinaus. Nur diese letztere kommt aber z. B. in der 
Störung dm^ d^ ä{l — g l) 

dt ^ dt~~ dt 
in Betracht, wovon maa sich leicht überzeugt, wenn man sich die Natur des 
Winkels w vergegenwärtigt, dessen Änderung'allein durch die Verrückung des 
Perihels und die Änderung der Knotenlänge cQ, bedingt ist. Übrigens würde, 

dif) 
wenn man auch die wahre Störung von -rr mit berücksitjhtigen wollte, die- 
selbe doch aus der Gleichung wieder herausfallen, da man in diesem Falle auch 
die Störung von l in Kechnung ziehen müfste, welche der eben betrachteten 
Störung von qp gleich ist, aber in der obigen Gleichung das entgegengesetzte 
Vorzeichen hat. — Es kommt also nur noch darauf an, die scheinbare Störung 
von (f mathematisch festzustellen. 

Jeder Vektor r liegt an der Grenze zweier Bahnelemente, eines, welches 
ihm vorangeht, und eines, welches ihm folgt. Von diesem Gesichtspunkte ge- 
hört r zwei benachbarten momentanen Ellipsen an, deren eine den K(»i- 
stanten 4 und p, deren andere den Konstanten * 4"^* "^^^ f-^^P entspricht. 
Derselbe müfste demnach gleichzeitig den Eelationen 

4* 
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r(l -|-€ cos y) «=« j) 
r[l-|-(t-j-(ff)co8y] -= V-^^V 
genügen. Dies ist aber nur dann ohne Widerspruch möglich, wenn auch dem 
<y> eine angemessene Änderung dip b3igelegt, d. h. das Perihel um diese Gro£9e 
verschoben wird — was inuuer zulässig erscheint, da </> eine blofse Bechnungs- 
grofse, nicht aber, wie r, eine von der Natur gegebene, unserer Willkür ent- 
zogene Gröfse darstellt. Man hat demnach als zweite Gleichung des Vektors : 

und ihre Verbindung mit der Gleichung ergiebt unsere Formel für 

lü' 

Dieselbe liefert offenbar dasjenige Inkrement von </•, welches zur fortwährenden 
Wiederherstellung der Ellipse erforderlich ist. 
Wenden wir dies auf die Gleichung 

Jlf s= m-\-ni 
an, so dürfte es wohl keinem Zweifel unterliegen, dafs der Störungsanteil 

dfn 

also die Störung der mittleren Anomalie der Epoche ebenfalls nur zur Korrektion 
der Ellipse dient, so dafs es begründet erscheint, bei ihrer Berechnung aus 
der obigen Gleichung 

dm r« dtf r( j>+r) . de 

~di ^ T~F^ 'dl "" "fT^ ^^^ di 

dw 
nur die scheinbare Störung -^ zu gebrauchen. 

Sollte es übrigens in einem gegebenen Falle nötig worden, die Total- 
störung von ff' zu wissen, so hat man in der vorstehenden, allgemein giltigen 

dm 
Gleichung nur für -^j den Wert der totalen Greschwindigkeit einzusetzen. Es 

folgt dann: 

^ff {p-\-^)0' de a* i^;p dm a* ^p fdn a» t/p 

M ^~ «"1 «P ^- + — T«— -M + —^^J dt <'*-^—^- ■«»■■(^) 

u 

und, nach Substitution der früher gefundenen Werte von 

df dm 

dt ^^ It'- 
dq> (P'\'r)re (p — r^p C 

Denn die beiden letzten Glieder der Gleichung (Ä) reduzieren sich auf — , • 
In der Praxis wird indessen die Gleichung (B) wohl kaum zur Verwendung 
kommen Dächte man etwa daran, dieselbe zur direkten Integration der wahren 
Anomalie zu benutzen — um den, namentlich durch die Auflösung der trans- 
szendenten Kepler sehen Gleichung beschwerlichen Umweg über die mittlere 
und eiLzentrische Anomalie zu vermeiden — so müfste die Ausführung dieses 
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C 
Gedankens an dem Gliede ^ scheitern, da in diesem der gestörte Badius 

genommen werden muXs, der selbst bereits die Kenntnis der gestörten Ano- 
malien Yoraoflsetzt. 

XVI. Entwickelung der stSrenden Kräfte. 

Um den vorhergehenden Formeln eine zur Integration ge- 
eignete Gestalt zu geben , müssen die Kraftkomponenten 

PcosQS, Pcos^JB, P cos QW 
noch als Funktion der Koordinaten der Planetenörter dargestellt 
werden. Diese Seitenkräfte ändern mit der fortschreitenden 
Bewegung des gestörten und des störenden Planeten nicht blofs 
ihre Gröfse sondern auch ihre Biebtüng^ Um die hierdurdi 
hervorgerufenen Schwierigkeiten unserer Aufgabe zu trennen, 
werden wir zunächst die nach drei festen Kichtungen zerlegtem 
Komponenten der störenden Kraft bestimmen und n^hstdem 
jene ersten Komponenten mit veränderlicher Richtung auf diese 
zurückführen. 

Wir wollen also die Kräfte bestimmen, durch welche die 
Bahn eines Planeten um die (als ruhend angenommene) Sonne 
bedingt wird, wenn au&er der gegenseitigen Anziehung dieser 
zwei Körper auch noch die Attraktion eines zweiten — des s. g. 
störenden Planeten — in die Bewegung eingreift. 

Nimmt man die Sonnenmasse zur Einheit und nennt 
m die Masse des gestörten Planeten, 
m| „ „ „ störenden „ 
r den Radiusvektor des gestörten Planeten, 
r| „ „ „ störenden „ 

X, y, z die rechtwinkl. Koordinaten des gestörten Planeten 1 . 
^uVu^i » » » » störenden „ J 

Bezug auf ein System dessen Ursprung in der Sonne 
liegt, 
(i die Entfernung zwischen dem gestörten und störenden 

Planeten, 
f die Krafteinheit (d. h. die beschleunigende Kraft der 
Masse 1 in der Entfernung 1), 
so sind es folgende vier Kräfte, durch welche der Lauf des ge- 
störten Planeten um die fest gedachte Sonne bestimmt wird; 



54 B. Theorie der speziellen Störungen. 

f 

1. die Kraft -%-, welche vom Centralkörper äüf de» gestörten 

. Planeten ausgeübt wird, 

2. „ „ -^-, mit welcher der gestörte Planet die Central- 

^ körper angreift, 

3. „ „ ^—^9 welche den gestörte» Korper nach dem 

^ störenden hinzieht, 

4. „ ,. •^— ^, welcher in Wirklichkeit den Centralkörper 

nach dem störenden treibt, aber hier — wo nur die relative 
Bewegung des gestörten Körpers gegen den unbewegtem Central** 
körper untersucht wird •— als eine Kraft behand^ werden lürnfs, 
welche den gestörten Körper in einer dem Vektor r^ paratteleD 
aber entgegen gesetzten Richtung zu bewegen strebt. 

Die Kräfte 1 und 2 lassen sich, da sie beide den gestöi'ten 
Körper dem Centralkörper zu nähern suchen, summieren, also 
in den Ausdruck 

/(l4-m) 

vereinigen. Werden nun diese Kräfte nach Richtung der recht- 
winkligen Koordinatenachsen zerlegt, so erhält man — wenn 
denjenigen Komponenten, welche den Abstand des gestörten 
Planeten vom Centralkörper zu vermindern streben, das negative 
Vorzeichen beigelegt wird — : 



d^ /(l-fm) X ß^ ^ Xx — x 

^ /(t-fm) ^ ßh h—^ 



fnii , ar, ] 




n' ' n 




fmi . Vi 




r,2 r. 




M ^1 




r,^ r. 





(I), 



dt' r- r (ß^ Q 

in welchen Gleichungen — , —^ Kosinusse bezeichnen, 

r Q 

deten Bedeutung aus einer. leicht zu entwerfenden Figur sofort 

erhellt. 

• Die ersten Glieder der rechten Seiten stellen die Kompo- 
nenten dar, welche den Planeten in einer reinen Ellipse um die 
Sonne herumführen würden, wenn alle anderen Kräfte mit einem 
Male zu wirken aufiiörten. Diese; durch die übrigen Glieder 
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ausgedrückten, Kräfte sind also die momentanen perturbieren- 
den Kräfte. 

Werden demnach die Störungen der Koordinaten in irgend 
einem MomeBte mit 

6x, dy, da 
angedeutet, so hat man: 



d*3x fmi(xi—x) 




<P<Jy ßhitfi-y) 


fmiPi 


dt* q' 


■ r,8 


d^dg fmt(gi—e) 
dt' p* 


fmtZt 



(!«)• 



Es sind dies nichts anderes als die rechtwinkligen KompoQenten 
der Kraft P. Um nun zunächst die Kraft 

Pcos^iS 
zu bestimmen, bezeichne C die Flächengeschwindigkeit des ge- 
störten Planeten, sowie 

die Projektionen von C in der Ebene der Y—Z, der X—Z und 
der F— X, wenn man, der Einfachheit wegen, jetzt und später 
die Ebene der F— X mit der in irgend einer Epoche statt- 
findenden Ebene der Ekliptik zusammenfallen läfst (s. Fig. 9). 
Sind 

hy hy h 

die Neigungswinkel der gestörten Bahn gegen die Koordinaten- 
ebene, so hat man bekanntlich nach den ersten Lehren von den 
Projektionen: 

Ct = C . cos i^ j 

02 = Ccwiij > • • • • (»*) 

C3 = Ccosia ) 

C2«Ci2 + Cj^2-f-C3^ (n). 

Gehen wir nun von der Gleichung 

aus, so erhalten wir durch Variation zunächst 
dC «dV T> na 
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Es ist aber nach Gleichung (n) 

dC C|^ dCj C2 dc2 C3 d(^ 
dt ^ C dt ^ C dt ^ C dt " 

so dafs 

Bedenkt man nun, dafs 



Ci-y 



dy 
dt 
dz 
dt 
dx 



dz 
dt 
dx 

dy 



X 



mithin nach dem Yariationsprinzip und den Gleichungen (la): 
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dc2 «PSy t Pde 

dt " " dP 

dCi d^dy _ d^dx 

dt "" dt* ^ dt'' 



(11), 



-dt-y-dp--'-^- -*^»».(y^i-^«) [t;t - p) 

80 folgt: 

^'C^{xy^ — yxJ],.,,{(t), 
Dies ist der gesuehte Ausdruck der Komponente, senkrecht zum 
Vektor in der Bahnebene. 

Gehen wir nun zu der in der Richtung des Vektors wir- 
kenden Seitenkraft 

Pcos^JB 
über. Offenbar erhält man dieselbe, wenn inan 4ie ?W)h den 
Parallelkoordinaten gezählten Kräfte 

d'dx 
.^^u.s.f. 

auf den Radiusvektor projiziert und die Projektionen- summiert- 
Nun sind aber die Richtungskosinusse des Vektors mit den 
Achsen 







±, y., 


z 






r r 


r 


Also hat man: 








P cos QR = 


X 

r 


dHx y 
dt^ r 


d^dy z d'^dz 
dP- ^ r dP 



^ — K^m^ 



(>3J 



...w, 



xx^ +yy^ +zz^ I 1 1^\ 

und damit auch den Ausdruck fiir die Vektor -Komponente. 

" Schliefslich ergiebt sich die zur Bahnebene senkirechte 
Seitenkraft Pcos^TT 

aus der Gleichung 

p cos QW - yp2 — P2 cos 2(2Sf - P* cos2 ^B . . . {y)> 
nachdem zuvor die Mittelkraft 

gefunden worden ist. 
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Anmerkung: 1. Encke bezieht, wie schoQ an einem ffdberen Ojft^ 
erwähnt, zum Zwecke dftf Sräfttbestitnmung die Hcmetenorter auf «in var- 
änderliclieg K-Oordinatensystem, dessen T— X-Ebene in der gestörten Bahnebene 
und dessen X-Achse in ^ Bichtunir ^» Badioftvektoirs fällt, Üemgfimäls 
ergiebt sieh beispiehweise für die Komponente Pcoei^S (vergl. Äerl. Astrött. 
Jahrb. für 1S37, S. 329) die Gleichung: 



PcoßQS^-k*m\{jj - ^j • ^,. 



Dieselbe Form nimmt aber auch uDAer AuBdtuck in (r<) an, wenn manbeachted]^ 
dafs im Encke sehen Koordinatensystem: 

a?=«r, y«0, = 

wird. Allerdings erscheinen unter diesen Umständen die StÖrungskompönenten 
in überaus einfacher Gestalt. Allein der so erlangte Vofpteil düiftor duroh di« 
umständliche Koordinaten-Transformatkm doch wohl zu teuer erkauft sein. 

Amnerkang 2. Die oben entwickelten Komponenten der störenden 
Kraft erscheinen — wenn man erwägt, dfef« 



dz dy . 



sämtlich als Funktionen der gestörten Koordinaten 

und der Geschwindigkeiten 

äx dy dz . 
If* 'W ~dt* 

ein Besultat, welches sich natfirlich leicht voraa«eehen liels. Die genaue 
Berechnung dieser Kräfte setzt demnach die Kenntnis der gestörten Koordi- 
naten bereits voraus, und es wärde ammö^lkh iseift) atf ^eseM Weg« sa irgend 
einem Ziele zu gelangen, wenn nicht der sehr kleine Störungsfaktor 

es gestattete, in erster Annäherung statt der gesfiötteb difi» eltipüsdlMi 
Koordinaten zu gehraucheik. Dasaelbe gilt von don in den obigjon Störungs- 
ausdrücken, z. B. in 

«^ k vp sm t 

auftretenden Gröfsen 

h P, i^Si, r u. s. f., 
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für welche ja strengenommen ebenfalls die gestörten Worte genommen werden 
mtifsten. Auch hier genügt es, zunächst die elliptischen — mit den oskulie- 
renden Elementen berechneten — Werte einzuführen. 

Man kann denmach für eine Beihe (gewöhnlich gleichweit abstehend 
gewählter) Zeitpunkte die Geschwindigkeit der Störungen numerisch berechnen 
und hieraus dann die Störungen selbst durch das Verfahren der mechanischen 
Quadratur bestinmien. Von Zeit zu Zeit wird man die Eechnung unterbrechen, 
um aus den gestörten Koordinaten sich neue Elemente zu verschaffen, welche 
dann fGr ein weiteres Integrationsintei^all die Basis der Bechnung bilden. 
Wann und wie oft eine solche Erneuerung der Elemente Platz zu greifen hat, 
dav«l«r wSa siftn ^^ so- bfilcl :al& uiögtioti, «nt^- i«*g]KiJi4g#t A'bwÄi^p^ aller 
im gegebenen Falle entscheidenden Verhältnisse, schlüssig machen. 



c. 

Anderer Weg zur Berechnung der speziellen 

Störungen. 



Die oben gegebene Entwickelung der Komponente 
P cos QS 
legt den Gedanken nahe, die Bestimmung der perturbierten 
Flächengeschwindigkeit 

C 
und ihrer Komponenten 

C|, Co, C3 

überhaupt zum Ausgangangspunkt der Störungsrechnung zu 
machen. 

Denn einerseits steht G in einer selir einfachen Verbindung 
mit dem Semiparameter jp, indem 

andererseits gehen die Neigungen 

iu i^y h (s. Fig. 9) 
leicht aus jenen Gröfsen hervor, zufolge der Gleichungen 

C4 .Co . C3 

cos tt *= 7T; <50S ^2 ** 7^ ; cos tj = yr • 

Sind aber i^ und ij gefunden, so läfst sich auch die veränderte 
Länge des Knotens K^, den man per excellentiam den Knoten der 
Bahn nennt, nämlich der Winkel 

ohne sonderliche Mühe feststellen. Denn da man die Lage der 
Achse OY gegen den Frühlingspunkt V kennt — es wäre ja ge- 
stattet, sie von von vorneherein mit der Äquinoktiallinie OV 
zusammenfallen zu lassen — da man ferner die Seite K^ Y nach 
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Ermittelung der Winkel i\ und i^ aus dem AKiK^Y finden 
kann, so ergiebt sich die gesuchte Knotenlänge durch die Glei- 
chung 

Der grofse Vorzug dieses Rechnungsgangs Kegt in der ge- 
radezu klassischen Einfachheit der Störungsausdrücke für die 
Flächengeschi^indigkeiten, insofern dieselben namentlich von allen 
Koeffizienten —• wie sin y u. ä. — befreit sind, ein Vorzug, der 
auf jedem anderen Wege verloren geht. Nach Bestimmung von 

folgt der Berechnung der Anomalien M, u und q>, wodurch 
gleichzeitig die gestörte Halbachse 

a 
bekannt wird, da die mittlere Anomalie die Kenntnis der ge- 
störten mittleren Bewegung 



-/; 



verlangt, diese aber bekanntlich gleich 

k 

ist. Aus p, a und ^ liefert sodann die Polargleichung den 
Vektor r. Um endlich auch die gestörte Breite ß, die Länge X, 
die Parallelkoordinaten x, y, z nebst ihren Geschwindigkeiten 
finden zu können, bedarf es nur noch des Winkels 

w, 
dessen Bestimmung — wenn diejenige der Knotenlänge voran- 
gegangen ist — in jedem Integrationsintervalle auf die Berech- 
nung der Perihelverschiebung hinauskommt. 

StSrimg der FlSchengeschwlndlgkelt. 

Bezeichnet (s. Fig. 9), wie bisher, 
G die Flächengeschwindigkeit, 
Ci ihre senkrechte Projektion in der F— ^-Ebene, 
^ » » » » » -^ ^~ » 

so hat man nach dem Obigen die Gleichungen: 
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-^ -= -^Mm-y^i)[y - ^s) ■ 

Sucht man also für verschiedene äquidistante Zeitpunkte mit 
Hilfe dar oskulierenden Elemente: 

die elliptischen Koordinaten x, y, des gestörten Planeten, 
„ „ „ i^uVu^t „ störenden 

den Radiusvektor r^ des störenden Planeten, 

die Entfernung q beider Planeten 
und berechnet damit die Geschwindigkeiten 

^«.8. f., 

dann liefert die mechanische Quadratur zunächst die gestörten 
Komponenten C|, Cj» ^3 ^^^ hierauf die* Gleichung 

die resultierende 'Flächengeschwindigkeit selbst. 

Störung der Neigung. 

Die vorstehende Bestimmung der Flächengeschwindigkeit C, 
oder da 

des Sempiparameters p erfordert drei selbständige Integrationen. 
Schon dieser Umstand lä&t erwarten, dals aus den gewonnenen 
noch zwei weitere Elemente gefunden werden können. In der 
That erhält man zunächst die Neigunswinkel der gestörten Bahn 
gegen die drei Koordinatenebenen vermöge der Gleichungen: 

cos ij « Cj : C, 

cos 15 — «2 « C> 
cos i^^c^: C. 
Wenn man, wie in Fig. 9, die Ebene X — Y als die Fun- 
damentalebene der Polarkoordinaten annimmt, so ist es inbeson- 
dere der Winkel 

welchen man die Neigung der Bahn nennt. 
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StSning des Knotens. 

Das dritte Element, welches aus unsem drei Integralen 
bestimmt werden kami, ist die Länge des Knotens Ki (siehe 
Fig. 9). Von der Lage der beiden Knoten K^ und -BTj können 
wir, weil fBr die weitere Rechnung ohne Bedeutung, hier absehen. 

Aus dem rechtwinkligen Dreiecke K^ K^T {olgi 

cos Kl Y« cos ii : sinij. 

Erwägt man nun, dafs auch die Länge yY des unserer Wahl 

überlassenen Punkts Y als bekannt zu betrachten ist, so ergiebt 

sich die gestörte (meist kurz durch ^ bezeichnete) Knotenlänge 

VZ| -= yY — arc cos (oosi, : sin tj). 



StSrung der mittleren Anomalie. 

Zur Bestimmung der gestörten Anomalie bedient man sich 
der früher entwickelten Gleichung (x): 



M^ 



•«•+/w*l+f"''+/'»/5'« 



in welcher 

«I« die mittlere Anomalie der Epoche (des Zeitpunktes d.08kulation) 

w, „ „ Bewegung „ „ und 

t die seit der Oskulation verlaufene Zeit bezeichnet, während 

am 1_ \U i> cos « 

dt ~ hV~a 11 



^] Pcos ^E+^siny Pcos QS 



{y) 



dn „ 
M — "* 



W 



lc-fä{\ e / ' e 

'>^ya yap 

zu setzen ist. 

Man wird also zunächst die den verschiedenen Zeitpunkten t ent- 
sprechenden Geschwindigkeiten der mittleren Anomalie 

_^ - . * 



dM dm 



fdn 
J dt 



dt 



= -3r + 



dt ai 



zu b^echnen haben, was berdts die Herstellung eines Integrals, 

k 
nämlich der gestörten mittleren Bewegung -j, erfordert. Zudem 

Ende kann man entweder die Gleichung (is) integrieren und so 
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die Störung der mittleren Bewegung unioittelbar feststellen, oder 
man kann vermöge der uns ebenfalls bekannten Gleichung: 

dt kr ^ k^p 

zuerst die gestörte mittlere Entfernung und daraus die mittlere 
Bewegung bestimmen. Nehmen wir das letztere an, so mufs 
nach dem Verfahren der mechanischen Quadratur aus den ver- 
schiedenen, in den gewählten Zeitpunkten stattfindenden Werten 

von 

da 

dt 
letztere Gröfse ak empirische Funktion der Zeit^ (vergl. u»A. 
Astromech. d. Verf, S. 28) dargestellt werden, ehe die Integra- 
tionen für die einzelnen Intervalle ausgeführt werden können. 
Erst nachdem so die verschiedenen Spezialwerte von 

dM 

dt 

gefunden und letzteres wiederum auf eine empirische Funktion 

der Zeit zurückgeführt, liefert eine nochmalige Integration den 

Wert der mittleren Anomalie 

M. 

Mit diesen Andeutungen müssen wir uns hier begnügen, indem 

wir wegen des Details der Rechnung auf Enckes Abhandlung 

„Über mechanische Quadratur", Berl. Jahrb. 1837, verweisen. 

SiSrung der Exzentrizität. 

Die Berechnung der mittleren Anomalie involviert bereits 
diejenige der mittleren Entfernung 

a. 
Da wir aufserdem auch den Semiparameter jp kennen, so folgt 
die Exzentrizität aus der Gleichung 

p « a(l — 6^). 

Störung der exzentrischen und wahren Anomalie* 

Die exzentrische Anomalie u ergiebt sich aus der mittleren 
M und der Exzentrizität e durch die Kepl ersehe Gleichung 
M — u — esinw. 
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während die wahre Anomalie 9> aus u und e mit Hilfe der 
Relation 

i/r=T.tg| - v/m-tgg 

gefunden wird. 

StSrHüg des Badlusrektors. 

Aus a, €, 9) oder a, «, m erhält man den Vektor mittelst 
der Polargleichungen: 

r(l+€cos q>) ^ p 
r = a(l — €cosu). 

Störung des Winkels w. 

Die unter No. 7 mitgeteilte Gleichung 

dco __ d^ d{l — Sl) 

dt dt dt 

d(p . dSi) 

oder in Integralform: 

(da man bei jeder einzelnen Integration cos i als beständig be- 
trachtet) lehrt, dafs die Berechnung der Störung des Winkels 
zwischen dem Perihele und aufsteigenden Knoten nur noch die 
Ermittelung des Integrals 

f^^dt - r(-(l±P>^^Pcos QS + ^^Pcos Qb] dt 
^ dt «/ \ ekyp ek\p I 

erfordert. Denn die Störung der Enotenlänge, d. h. das Integral 

ist bereits bekannt. 

Störung des Arguments der Breite« 

Das Argument der Breite K^B (s. Fig. 9) ist 
=* </) — (jt) « l — S^, 
wenn, wie bisher, unter l die Länge ^K^B des Planeten P in 
der Bahn verstanden wird. Da die Störungen von 

Israel-Holtzvart. ^ 
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qp, CO und Sl 
bereits gefunden sind, so hat man sowohl die Störung 

des Arguments 9>— w 
als auch der Länge l in der Bahn. 

Störung der Breite. 

Aus dem Dreiecke K^BM (Fig. 9) erhält man die gestörte 
Breite ß durch die Gleichung 

sin /? = sin 4 sin (y — w). 

StSrung der LSnge. 

Die gestörte Länge in der Fundamentalebene ist 

Die Knotenlänge fi ist bereits bekannt, der Winkel K^M aber 
durch die Gleichung gegeben: 

tg ^i-3f = cos ^3 • tg (y — (X)). 

StSmiig der Parallelkoordinaten x, y, z. 

Zuerst findet man die Koordinate 
z = rsin/?. 
Zur Bestimmung der beiden Koordinaten 

y und z 
wird man die Gleichungen in Anwendung bringen können: 
r2 c=: a;2 + t/* + ^2 

CjiC + Cay + Ca^ei^O. 
Die letztere Belation — welche bei konstantem Ci, c^ und C3 
ausspricht, dafe die Koordinaten Xy y, z einer Ebene genügen, 
die Planetenbewegung also in einer Ebene vor sich geht — er- 
giebt sich sehr einfach, wenn man die erste der Gleichungen 
dz dy 

dy dx 

mit Xy die zweite mit y und die dritte mit z multipliziert und 
die Summe der drei Gleichungen bildet. 
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Einfacher zum Ziele kommt man mit den Gleichungen 
IS ^ rsinß \ 

y «= rcosjSsin(X — ©) ^ (m), 

X ^ r cos ß sin {X — ©) j 
wenn mit der Winkel VX (s. Fig. 9) angedeutet wird, den 
man gleich Null setzt, im Falle die Achse der x durch den Frühlings* 
punkt V gelegt wird. 

StSrung der Koordinatengesohwlndigkeiten. 

1. Folarkoordinaten. 
Die Geschwindigkeiten der wahren Anomalie und des 
Radiusvektors sind wie in der elliptischen Bewegung, durch 
die Gleichungen ^ 

d(p h^p 
dt ^ r^ 
dr he , .,. 

äT-y^«^'^^ <^) 

gegeben; nur sind für p, «, r und tp die gestörten Werte zu 
nehmen. Ersetzt man in den Gleichungen sub II und III der 
Hilfsätze 

|aurchg(pi,.6), 

SO ergiebt sich für die Geschwindigkeit der Länge 

dl __ cosi d(p 

dt ~" cos ^ß dt ' 
also dl k }/p • cos i .^^. 

dt r^cös^ß~ ^^ 

und für die Geschwindigkeit der Breite 

oder dß _ ky^ptgßcotg (l—S l) 

dt r^ " 



(3). 



2. Parallelkoordinaten. 

Aus den Geschwindigkeiten der Polarkoordinaten folgen 
die der Parallelkoordinaten: 
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dz =* sin ßdr'\'r cos ß dß 

dy «= cos ßmiXdr + r cos ß cos l dl — r sin A sin /!/ dß 
dx = cos ß cosk dr — r sin |ö? cos X dß — rcosßsiukdl, 
wenn man in den Gleichungen (m) der Einfachheit wegen 

annimmt. Übrigens könnte man zu dem yorliegenden Zwecke 
auch die Gleichungen der Flächengeschwindigkeiten benutzen: 

dz dy 

dx dz 

dy dx 

in Verbindung mit der Relation: 

dr ke . x dx y dy z dz 

dt |/j) ^ r dt r dt r dt 

Denn die drei ersten Gleichungen sind für sich allein zur Be- 
stimmung der Geschwindigkeiten nicht ausreichend — infolge 
der zwischen ihnen stattfindenden Beziehung 
CiX -{- c^y -\- c^z ^ 0. 



D. 

I 

Die Störung der Flächengeschwindigkeit und 
der mittleren Entfernung als Basis der ge- 
samten Störungsrechnung. 



Für die eben abgeschlossenen Entwickelungen ist besonders 
bezeichnend, dafs die zur Bahnebene senkrechte Störungskom- 
ponente 

Pcos QW 
in dieser expliziten Form überhaupt nicht zur Anwendung kommt. 
Denn die einzigen Elemente, deren Störung durch sie bedingt ist, 

die Neigung und Knotenlänge, 
werden im vorliegenden Falle aus den gestörten Flächengeschwindig- 
keiten und diese selbst unmittelbar aus den rechtwinkligen Kom- 
ponenten der störenden Kraft hergeleitet. 

Aber auch alle anderen Störungen lassen sich leicht, mit 
Umgehung der Kräfte 

PcosQB und PcosQS, 
durch die orthogonalen Komponenten 

d^dx d^dy dHz 
dt^ ' dt^ ' dt^ 
ausdrücken, oder, richtiger gesagt: Man kann ebenso wie die 
Störung der Flächengeschwindigkeit und ihrer Komponenten 
auch diejenige der mittleren Entfernung in sehr einfacher Weise 
durch die rechtwinkligen Komponenten darstellen und dann alle 
weiteren Störungsgeschwindinkeiten auf diese beiden 



zurückführen. 



dC , da 

■dt ""^ dt 
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Das Charakteristische dieses Verfahrens besteht also darin, 
dalfi an Stelle der nach Richtung der Polarkoordinaten zerlegten 
Kräfte 

Pcos QW, Pcos QS, Pcos QB 

gewisse Funktionen derselben, insbesondere die Störun^sch windig- 
keiten dC .da 

eingeführt werden. 

Um eine vollständige Übersidit der gesamten Rechnung 
zu gewinnen, sei hier nochmals daran erinnert, dafs aus den 
Störungen der Flächengeschwindigkeiten 

in der oben gezeigten Weise zunächst die drei Elemente 

Parameter, Neigung und Knotenlänge 
entwickelt werden. 

Da ferner aus der Gleichung 

sich ergiebt: 

dp _ 2 ^dC ^ 

dt k^ dt' 
aber nach dem Früheren: 
^dC dci . dco dc^ 

^It^'^-df-^'^^-^'^^ 



(«i) 



"^'^^if^ä ~ T3]ki(2^^i-~^2/i)+C2(^^i~-a;^i)+C3(^yi— y^i)] • • -(«2) 



so ist auch 



dp , ^dO 

dt ""^'^ ^ dt 



durch die vorangegangene Rechnung ohne Weiteres als gegeben 
zu betrachten. 

Die weitere Rechnung gestaltet sich dann folgender- 
mafs^u: 



Störung der mittlerou Entfernung«, der Exzentrizität s. 

1. Störung der nüttlereu Entfernung a. 

Durch Variation der bekannten Gleichung: 
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a 



wird erhalten: 
da 
dt ° 



2^ 
r 

r 






m-&h(^)' 



dx dHx dy ^ d^dy de 
di"W^di' ~W^di 



dHz 



df 



oder nach (la), Seite 55: 



da 
dt 



=» 2»!«« 



/ dx 



^y^dt^'^dt 



dg\ n _ n 
dt)\Q' n^j 

, 6?y , d0\ 1 



/ dic dy ^ . 



(««). 



Vermöge der beiden Gleichungen («2) und (ag) ist es nunmehr 
gestattet, ^e Geschwindigkeiten 

dp , da 

als bekannte Gröfsen zu behandeln. 

2. Störung der Exzentrizität e. 

Obwohl an dieser Stelle der Semiparameter p und die mitt- 
lere Entfernung a bereits gefunden sind, demnach die gestörte 
Exzentrizität einfach durch die Relation 

p-^a(l — €^) 
bestimmbar ist, so findet doch auch die Geschwindigkeit 
dieser letzteren Störung in Zukunft als blofse Hilfsgröfse Ver- 
wendung, nämlich als Funktion der Fundamentalgeschwindigkeiten 
dp . da 

dt "'^^ dt- 



Die vorstehende Gleichung liefert: 
de 1 l da 

dt 



2a'^e 



/ da dp\ 

[^dt-''di) 



(«*). 



3. Störung des Perlhels. 

Nach dem Früheren hat man: 
d(p 1 

dt E sin (f 



cos 



d€_ 1^ dp\ 

"^'dt r dt) 



(«5)- 
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4. StSrung des Winkels lo zwischen Perlhel und Knoten. 

Es ist, wie oben bereits gezeigt: 

/w*-/^*--/f« «-.)■ 

Da wir die gestörte Knotenlänge, also das letzte Glied der vor- 
liegenden Gleichung, schon kennen, so bedarf es nur noch der 
Herstellung des Integrals 



/ 



dt 



auf Grund der Relation («5). 

5. Störung der mittleren Anomalie mo der Epoche. 

Dieselbe ergiebt sich aus der im Vorhergehenden ebenfalls 
begründeten Gleichung: 

dm r / r d(p p + r . de\ , . 

-di--Y^[Tlt--Y-'''"^dt) <'''>• 

6. StSrang der mittleren Bewegung Wq der Epoche. 

Im Falle die mittlere Entfernung a bereits festgestellt ist, 
was wir hier voraussetzen, so erhält man die gestörte mittlere 
Bewegung n durch die Gleichung 

k 

Eine direkte Bestimmung der Störung yon n^ würde mit der 
Gleichung: 

dn ^ ^ ^ 

dl '^ ~2~J'm 
auszuiühren sein. 

7. Störung der mittleren Anomalie M zur Zelt t 
nach der Epoche, 

Die Wirkung der beiden vorhergehenden Störungen zeigt 
sich in der mittleren Anomalie Jf, welche zur Zeit t nach der 
Epoche durch die Gleichung gegeben ist: 



E. 

Die Störungen 

der Parallelkoordinaten und Elemente, 

dargestellt durch allgemeine Reihen. 



Die EntwickeluDgen der Abschnitte J5, und D betrafen 
ausschliefslich die zur Berechnung spezieller Störungen geeig- 
neten Formeln. Dieselben liefern die Abweichungen von der 
elliptischen Bahn nicht — wie dies in der Theorie der sogen, 
absoluten Störungen der Fall ist — in der Form allgemeiner 
Funktionen der Zeit, sondern geben nur die Geschwindig- 
keiten der Störungen als Funktionen von Gröfsen, welche mit 
der Zeit sich ändern und für jede Zeit mit einem beliebig hohen, 
wenn auch erst allmählich erreichbaren Grade von Annäherung 
bestimmt werden können, so dafs die Möglichkeit vorliegt, aus 
diesen für eine Reihe von Zeitpunkten berechneten Geschwindig- 
keiten durch mechanisches Quadrieren die Beträge der Störungen 
selbst zu ermitteln — zunächst für einen beschränkten, dann 
aber für jeden noch so grofsen Zeitraum, der ja stets in Inter- 
valle von angemessener Kleinheit zerlegt werden kann. 

Es liegt nun zwar unserem Zwecke fern, auch die absoluten, 
d. h. allgemeinen Störungen in den Bereich unserer Betrachtungen 
zu ziehen. Doch glauben wir dem Leser einen Dienst zu er- 
weisen, wenn wir zum Schlüsse seine Aufmerksamkeit noch auf 
einige Reihen lenken, welche die Parallelkoordinaten eines ge- 
störten Körpers mit völliger Allgemeinheit, selbst mit Berück- 
sichtigung der Störungen des störenden Himmelskörpers, als 
Funktionen der Zeit und gewisser Anfangskonstanten darstellen, 
welche deshalb als die allgemeinen Integrale der Differential- 
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gleichungen des s. g Problems der drei Körper zu betrachten 
sind. Schon aus diesem Grunde und weil ein anderes allgemeines 
Integral bis jetzt nicht hergestellt ist, wohl auch kaum jemals 
gefunden werden wird, dürften jene Reihen ein erhöhtes Inter- 
esse für sich in Anspruch nehmen. Daran ändert nichts der 
Umstand, dafs jene Reihen wegen der im Verlaufe der Zeit 
mehr und mehr abnehmenden Konvergenz zu allgemeinen Rech- 
nungen unbrauchbar sind — zumal sie vielleicht das einfachste 
und natürlichste Mittel darbieten, zunächst für kleinere Zeit- 
abschnitte und nächstdem für willkürlich grofse, aus solchen 
Intervallen zusammengesetzte Zeiträume • die numerischen Werte 
der Störungen aufzufinden. 

Nehmen wir aü, für den Anfangspunkt der Zeit t (den 
Oskulationszeitpunkt) seien gegeben: 

1. die Koordinaten und Geschwindigkeiten des gestörten 
Körpers 

dx„ dyo d0o 
dt' dt ' dt 

2. dieselben Gröfsen für den störenden Körper 



^0 > 



dxQ^ dyo^ d0o^ 
dt ' dt ' dt 
so sind nach dem Gravitationsgesetze auch die zweiten Diflferential- 
quotienten als bekannt zu betrachten, da 



df- 

df" 

d%^_ 
dt' 

d^Xj^ 

dt'' ■" 
df' 



f(i-\-m) Xq 



+ 



' 

/(l+m) 



Co' 

f^iiyo^-yn) 



I ßhi^o'- 

MyjLzijill 



f(\+nn\yo^ 
' 




ai). 



dargestellt durch allgMueine Beihen. 
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Alle höheren Diffeireatialqaotieateii aber lassen sich gleichMs 
anf die Koordinaten und ihre Geschwindigkeiten zurückfuhren. 
So ist beispielsweise der Quotient dritter Ordnung 






/•(i 



+m) 



fm 






dt 



drg _ 
dt 

_ dx^ 
dt 



dx^ 
dt 



+ 



/Wj 



3x„^ 



dt 



dxo' 



dt 



dt 



dt 



_ 3/mi(a;o'— a;„) dp„ 
(?o* dt 

dr 
wobei man sich — in Bezug auf die Bildung von -~9 

— nur zu erinnern hat, dafs 

Die in den Quotienten vierter und noch höherer Ordnung immer 
wieder erscheinenden Quotienten 

cPx„ cPxq « 

df- ' dt^ '^•^•^• 
sind selbstverständlich stets durch ihre in (I) und (II) gegebenen 
Ausdrücke zu beseitigen. 

Betrachten wir nun die Koordinaten x, y, ^ als Funktionen 
der Zeit t, indem wir setzen 

x^m 

g-^xpit) 
und entwickeln, yrie zuerst Euler gethan, nach dem Lehrsatze 
des Maklaurin, so ergiebt sich: 



X'=Xi 



■*" dt 



t . d^x^ t^ d%,. t^ 



y-Vo-^ 



dVo 



■^n + 



dt 
+ 

deo 



dt 



1 ■ 

d^ 
dt'' 

■\- 

dt'' 
^4- 



dt 
t" 

dt^ 
nl ^ 



+ 



+ 



1-2 • dt* 1-2-3 
d"+^ f(art ) ft+t 

d^»+i ■ (n+1)! 
t^ . d'y„ t* 



+ 



1-2 ' dt* 1-2-3 
dt"+^ ' ' (n+1)! 



ä!»«„ . ** . d'«„ t* 



1 
dP* 



dt^ 
nl 



+ 



1-2 dt' 1.2-3 



(?^»+i 



(n+1)! 



(III). 
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Die letzten Glieder dieser Gleichungen, in welchen a^, a^, «3 
drei positive Zahlen zwischen Null und der Einheit bedeuten, 
bilden die s. g. Reste der Reihen. Dieselben sind konvergent 
unter den für die Maklaurinsche Reihe giltigen Bedingungen. 
Vergleicht man sie mit der Exponentialreihe 
t t^ t^ 

^=n._+_+__+ 

welche für jeden endlichen Wert von t konvergiert, so erkennt 
man leicht, dafs auch die Gleichungen (III) eine endliche Grenze 
haben, wenn sämtliche Differentialquotienten 

. ^«' "dT' ~W' dt^' 

,, äyo_ d^ dVo. 

^«' dt' dt^ ' dt* ' 

d^^ dH^ dH^ 

^«' dt' dt''' dt*' 

endliche Werte haben. 

Da im vorliegenden Falle die Koordinaten 

^, y, ^ 
aus den zweiten, durch das Gravitationsgesetz gegebenen 
Differentialquotienten 

d^ d^ d^ 
df" ' dt^ ' dt^ 
erzeugt werden, so sind die als bekannt vorausgesetzten Wert- 
paare dx^ 

^^' ~dt 

^<>' dt 

^«' dt 
als die Integrationskonstanten zu betrachten. 

Man kann aber der Entwickelung auch eine andere, in 
manchen Fällen vorzuziehende Form geben. Bezeichnet man 
nämlich die Integrationskonstanten vorläufig mit 

A2, X>2 
■^3 » -^3» 

SO darf man setzen: 
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x = A^-\-Bit +fätf^ dt 
y = ^2 + B4 +fdtf^ dt 
0^A^-hB^t-\-fdtf 



d^0 



dt 



dt^ 
Zerlegt man nun mit Hilfe der Formel der teilweisen Integration 

fu dv « UV — yV du 
die Doppelintegrale in zwei einfache Integrale, indem 

dv'^'dt 

U.S. f. setzte so erhält man weiter: 

,^A, + B,t + tf^dt-ft^dt 

Wendet man hierauf zur Entwickelung der einfachen Integrale 
entweder viederum die partielle Integration oder auch die Reihe 
des Joh. Bernoulli an, so ergiebt sich:- 
/•(T'a; r< d^x t^ ä^x t' d*x 

J dt* 1 dt* 1-2 dt* "^ 1-2-3 dt* "*" 

^ d'+^x ft^ d '+^-x ," 
=^«I d^+» ^J»i d^+2 *** 

_ r.d^ i^d^ t* <Pa; t* d*x 

J dt* ~ 1-2 dt*'^h2'i dt* 4:1 dt* '^"' 



J (n-{ 



Damit wird: 

. _^Ti._^i-t^d*x 2t* d'x 3t*x 



dt 






^ (»+1) 1 dt'>+* 

wo die oberen oder unteren Vorzeichen zu nehmen sind, je nach 
dem n ungerad oder gerad ist. 
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y 



■B4- 



31 



d^x 



3t* 

41 



d*x 
dt* 



Ebenso findet man: 

21 dt 
g^A, + B,t + ^-^- 

Man überzeugt sich leicht, dafs die Integrationskonstanten 
-Ä| , A2 , A^ 
mit den zur Zeit ^«=0 stattfindenden Werten 

der Koordinaten, hingegen die Integrationskonstanten 

mit den Anfangswerten 

^, ^, ^ 
dt dt dt 

der Geschwindigkeiten zusammenfallen. 
Hiernach folgt schliefslich: 



y = yo 



eteo 


. . <^'^ 


1-f* 


d'x 


2*» 


_J iL 


•^+d<^ 


1-2 
. d*x 


dt* 
St* 


31 






+ dt* 


4! 


"■ 


^ dt 




12 


dt^ 


2t* 
3! 






"*" dt* 


St* 
4! 


— 


. d^o 


_ d^e 


l-f 


d^g 


2t* 


+ d^ 


•^^ dt^ 


1-2 


dt* 
St* 


31 






' dt* 


41 


■~ 



(IV). 



Die Bedingungen der Konvergenz sind dieselben wie die der 
Reihen (III). Beide Reihensysteme unterscheiden sich nur da- 
durch, dafs in dem einen sich sämtliche Diflferentialquotienten 
auf den Anfangspunkt^ in dem anderen auf den Endpunkt des 
Zeitintervalls beziehen, mit Ausnahme der als Integrationskon- 
stanten zu betrachtenden Gröfsen 

dxn 



X 



e » 



dt' ^«' ' 

welche stets dem Anfangspunkte zugehören. 
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Obgleich den Reihen (III) und (IV) der Charakter eines 
allgemeinen Integrals zukommt, so ist ihre Anwendung zu all- 
gemeinen numerischen Rechnungen doch, schon durch die höchst 
umständliche Entwickelung der höheren Differentialquotienten 
ausgeschlossen. Hingegen zu speziellen Rechnungen scheinen 
die Reihen um so besser geeignet. Wir wissen, dafs irgend eine 
momentane Ellipse völlig bestimmt ist, wenn die ihrer Zeit ent- 
sprechenden Koordinaten nebst Geschwindigkeiten, also: 

dx dy dz 

dt dt dt 
bekannt sind. Die Rechnung wird sich deshalb für die einzelnen 
zum voraus bestimmten Zeitpunkte lediglich auf diese Gröfsen 
zu richten haben. Wählt man diese Zeitpunkte einander nahe 
genug (in dem unten folgenden Beispiele um je 20 Tage aus- 
einanderliegend), so wird man die Reihen kaum weiter als bis 
zum dritten Differentialquotient einschl. zu berücksichtigen haben 
— sofern es sich um die Störungen handelt. Dennoch würde 
selbst die gewöhnliche Berechnung dieses letzteren Quotienten 
zu beträchtlichen W^eiterungen führen. Wir woUen deshalb — 
ehe wir zur Berechnung eines konkreten Falles übergehen — 
zunächst einen sehr einfachen Weg andeuten, auf dem man zu 
den Werten der dritten und vierten Quotienten gelangen kann, 
ein Verfahren, welches übrigens — wie wir in dem Abschnitte 
„Störung der Elemente" sehen werden — bis zu jedem ver- 
langten Grade astronomischer Genaunigkeit erweitert werden 
kann. Bezeichnet man die zu drei verschiedenen, durch die 
Intervalle t^ und t^ 

getrennten Zeiten stattfindenden Differentialquotienten zweiter 
Ordnung durch 

dP^ d^ d^ 
dt'^ ' df' ' dt^ ' 
so hat man nach dem Maclaurinschen Satze: 

d^2 df^ dt' ' l '^ dt'' ' 2 d^5 * 6 "*■ 

d^ ^ d^ d^ t^ d^x^ t^ d^x^ ^ 

woraus bis auf Gröfsen fünfter Ordnung sich ergiebt; 
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d*X2 
~1F 

d*X2 



_ _ _2 

dt* ~ ij^ih + t^) 



** IF 



-<,' 



t.^ + t, 



dt^ (*i «2 ) 



ih + h) 






- dp ' ' dP ^'' • "'' dp 
Wählt man, wie gewöhnlich, gleiche ZeitintervaUe, so verein- 
fachen sich die Gleichungen in: 



d^Xi 
~dP 
d*Xi 
1F 



1 

2ti 

1 


d'^Xi 
dp 
d^Xi 

[dP 



d'^x^ 



dP 
d^Xj 
dp 



— 2 



ä^Xi 
~dP 



(V). 



Dieselben bieten dann ein sehr bequemes Mittel, um aus drei 
aufeinanderfolgenden Differentialquotienten zweiter Ordnung 
den dritten und vierten Quotienten des mittleren Orts zu 
finden. Wird nur eine Genauigkeit bis zu den Gröfsen vierter 
Ordnung verlangt, was der Kegel nach völlig befriedigende 
Resultate liefern dürfte, so kann man nach den obigen Reihen 
auch setzen: 



d%2 



1 

ti 

1 


d'^Xz 
dP 
d^Xi 
dP 



d^Xi 



dP 



dp 



oder auch: 



(Vax 



Berechnung der EoordinatenstSrungen. 

Subtrahiert man von der gestörten Koordinate 



CvOCq r^ Cb ccq 



^•=^-+if-i + 



dp 



t 2 



d^Xg 
dP 



f 3 

6 ^ 



ihren. ungestörten elliptischen Wert 

d^ ti . dH^ 
dt 



Sl — fco + jj. 






— + 
2 ^ 



d'lo 



6 



dp 2 • dt* 

80 ergiebt sich, wenn man beachtet, dais der Index Null sich 
auf den Zeitpunkt der Oskulation bezieht, demnach die Glei- 



chungen bestehen: 



die Beihe: 



yo-l. 

dx, ^ d^o 
dt dt 



d'(a;o— So) 
dt* 



2 



+ 



dK^o-io) <,' 



dt* 



T-»- 



Berechnnng der kuordinatenstörungen. 
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Wird die Störung 
allgemein durch 



X- 



Ax 



bezeichnet, so folgt für die Störung am Ende des 1. Zeitinter- 
valls nach der Oskulation: 



AaTi — 



dt^ 2 ^ dt* ■ 6 "*" 



Durch die Euler sehe Reihe erhält man für die Werte der 
Störungen am Ende der folgenden Intervalle: 



AXi = A«! + 
Aa;,— AXi+ 



dAXj 

dt 
d AXj 



ti_ ^Ax^ ^ d^AXj 
2 ■*" dt^" ' ^ ^ dt' 



t. 



d^Ax. 



dt^ 



2 ti^ d*AXi 



dt> 



6 ^ 
6 ^ 



z^Xn ^^,-1 ^ — 2 + ^^, 2 ^ dt^ 6^ ^ 

Werden auf die weiter unten, bei Betrachtung der Störungen 
der Elemente, angegebene Art die dritten, vierten und fünften 
Differentialquotienten durch solche der 2. Ordnung ausgedrückt, 
so verwandeln sich die Gleichungen in: 



Aa;i = 



'360 



d*Ax -i 
~~di* 



36 



d*Ax-i 
dt' 



. ,^, d*Axo , ao d'Axt 



_ A t^Aa:, ti ti* 



^ •d'Aa;-! _ „g d'Zia;^ 



d«' 



+ 171 -^5'- + 38 



d*Aa;2 



Aa; -Ax'^-i+^^^.k+Jll. 
Aa;„-Aa; + ^^ ^ +3^^ 



df 



,d*Aar„_3_^ 



dt' 

d'AXn-i 



+ 171 



df' 

d»Aa:,_i 
dt' 



+ 38 



d<' 
d'Zia;, 



di' 



(I), 



WO die Stellenzeiger — 2 und — 1 auf die Anfänge des 2. und 
1. Zeitabschnitts vor der Oskulation hinweisen. 



Israel-Eoltzwarfc. 
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Entwickelt man nun ferner auch clie ersten Differential 
quotienten nach Potenzen der Zeit: 
dAXj d'Axo ^ d»Aa;o <|^ _ 

dt '^ dt' ' 1'^ dt^ ''2 '^ 

ti I d'Ax-i d*Ax-i d'Z\a;„ dV, \ 

24 \ dt^ dt^ "^'*^ dt* "■*" df I 

dAcCi _ dAxi , d^Axi <, , d'Axi f^* , _ dAxt 



dt 



d AXj d*Z\Xi tj d^AXj tj* 
dt dt* ' l'^ dt* " ' 2 



dt 



tt (d*Ax-i ^ d*AXa . ,^ d'Axt . ^ d*A^i \ 
24\ dt* ^ dt* ^^^ dt* '^^~dt* f 



dA^n-\ dAXn-1 d*AX n-i tj d*AXn-i V , 

dt ~ dt "^ dt* ' 1 '*' dt* ■ 2 "^ 



dAx, 



■«-2 



dt 



ty (d^Ax„.-i . d* Ax„-3 • , „ d*Axn'-i „d'Zia;«-i\ 
'^2i\ dt* ~^ dt* '^^^ dt* ^^~W~) 



(la). 



80 bekommt man durch Addition: 



dAx, 



n-l 



dt 



24 



( 



d*Ax-2 d*Az-i 
dt* "^ dt* 

d*Ax-i d*A x^ 
dt* "^ dt*' 



, d*AXn-i\ 

^—dt*—) 

d*AXn-S 



+ 19 
+ 9 



d*Ax„ 



dt* 



+ 



d*Axi 



I d*AXi 
\ dt*~ 



dt* 

d*Ax^ 



+ 



dt* 



dt 



5-- + -4 



d*AXn-i \ 

dt* j 

d*AXn-l 



dt* 



<ißy 



Die Substitution in die letzte der Gleichungen (I) liefert endlich: 



AX„— AiCa-l = + -—^ 



'd*Ax-i 



df. 
.(d*Ax-i 
- ^[—df-- 
d*Ax„ 
dt^ 



+ 
+ • 



d 



'Aa;.-, \ 
dt^ I 

d*Ax^% 



+ 19 



+ 9 



3150 



'Aa;n-a \ 
dP j 

( d*Axy d*Ax,-i \] 

\ dt^ -»-•••+ ^fi jj 

\jd* AxH-3 _ 3g d*Axn--i .„. 



( 



+ 38 



dt' 

d*AXn 



d*AXn-i 

dt^ 

■ ■ ■ m 
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Hierdurch ist die im Laufe des n. Zeitintervalls nach der Osku- 
lation eintretende Störung 

zurückgeführt auf eine Summe von Diflferentialquotienten der 
2. Ordnung, ißt also bestimmbar, sobald die letzteren ge- 
geben sind. 

Beispielsweise hat man hierdurch für die Störung während 
des 3. Zeitintervalls: 

/\x, - A^2 - ^ [\dF~ "*■ ~~di^ "^ ~~dt~) 

\ dt' "^ dt' "^ dt' ) \ dt' "^ dt' "*" dt' ) 

_^ c,(d'Ax, _^ d'Ax^ ^ d'Ax^ \ 
'^^\~dr' '^ dt' ^ dt'~~) 
- d'Ax^ ^ d'Ax^ , .„. d'Ax^ ^^ d'Ax^ 



^360 

Ähnliche Formeln ergeben sich für die beiden anderen Koor- 
dinaten. 

Auch unterliegt es keinen Schwierigkeiten die Störungs- 
formeln für andere Zeitpunkte als die Grenzen der Intervalle 
aufzustellen. So findet sich für die Mitte des zwischen und 1 
liegenden Intervalls (worunter man keineswegs gerade das der 
Oskulation folgende Intervall zu verstehen braucht), wenn man 
rückwärts interpoliert: 

Die Geschwindigkeit der Störung: 
dAx Q^i ^ dAXj ___ d 'AXj j, d^A Xj t^' _ 

dt "" dt "~dt* * 2 ■*■ dt^ '2.2» 

__ dAx^ ^ tt L d'Ax^i -o d'^Ax. , f;^ ^^Z^, ^ d'Ax^] 

^^W '^m[^~dt' ^^"d"^~ ^^~dt~^^^^li^y 

Die Störung selbst: 

dAXi t^ . d'Axi t^' d^Ax^ t^^ . 
Ax^^^Ax,^-^ .-^^——^•^^^-.-^^^^— + ... 

"^^* dt 2^11520 

-87^-i + 261^^+1269^^-53^^^ 



dt' "^ dt' ^ '^^-^ dt' dt' 
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Entwickelang der Differentialquotienten 3, Ordnung. 

Wir müssen nun die Quotienten 2. Ordnung, aus denen 
sämtliche Störungen herzuleiten sind, näher betrachten. Wie 
bekannt, ist die den gestörten Planeten angreifende Gesamtkraft 
d^x Ä'u+m)^ Tc^m^ x^ lc^m\x^ — x) 

Durch Subtraktion des elliptischen Betrags . 



ergiebt sich die störende Kraft 

d^Zso) ,. (x |\ Ic^m^x^ 



dt^ 



Wird 






Iz^m^x^ — x) 



rl3 



ß3 



xr=^-\'Ax 

(r,+Ar)» ^(§+Aa;)M-(i? + Ay)^ + (^+A^)^ 
gesetzt und läfst man die höheren Potenzen der Störungen vor- 
läufig aufser acht (was wegen des Faktors h^ auf eine Vernach- 
lässigung von Gröfsen 6. Ordnung hinausläuft), so folgt 

S"- A-=(^*'~3g^^^a;-3?^.Ay-.3Sr.A;2; 
r r^ 

und damit der Ausdruck der störenden Kraft in der Richtung 

der a:- Achse 

d^Ax i»' 



(l+m) 



dt^ 



[(r.'-3g")A^-3J^;Ay-3gC-A^] 



ri^ 



und analog in der Richtung der beiden anderen 
Achsen: 

-^~ - **^"^ l(n'-3i/*)Ay-3i?^Aa;-3v^A;^J 



dt* 


r,' 


d*As 


*'(l+m) 


dt* ~ 


r/ 



L.13 



Q' 






rl3 



Q* 



(III). 
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Man ersieht hieraus, dafs die störenden Kräfte der Koordinaten 
bereits von den Störungen der Koordinaten selbst abhängig sind. 
Doch können letztere wegen des kleinen Faktors V = 0,0003 in 
erster Annäherung unberücksichtigt bleiben. 

Bei Aufstellung der Gleichungen (I) und (la) wurde vor- 
ausgesetzt, dafs 

Ao;, und~^ 

sich unmittelbar an die Oskulation anschlösse. Ist dies nicht der 
Fall, besteht bereits die Störung Aa;© und die Geschwindigkeit 

—77^, so ist der Wert von A x^ durch A x^ -I — -rr^ • -r- , der 
von —Tf^ durch — -, ® und folgeweise der Wert von Axn — Axn-\. 
in Gleichung (II) durch — ^r^ • v- zu ergänzen. 



Berechnung der Störungen. 

Erste Annäherung. 

Bezeichnet man Kürze halber die in den Gleichungen (III) 
vorkommenden Gröfsen 

^[(r,»-3|").Aa; ] durch S^ 

i^ - ^ durch Q' 
i_[^^._3,2).Ay ] durch 8^ 

U. 8. f., 

80 hat man: 






(in«). 
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Durch 



Sj^,. S_\, S/, 8- 

Ql^. «A, «0* 

Si,, 81, 

u. s. f. 
werden alsdann die zu den verschiedenen Zeiten stattfindenden 
Werte dieser Gröfsen angedeutet, und zwar für gewöhnlich 
durch 8-2 der um 2 Intervalle vor der Oskulation bestehende 
Wert u. s. f. 

Mit diesen Abkürzungen erscheint die Gleichung (II) in 
der Form: 

fPmV 



n— 4 »— S n— 2 n— 1 

-8—1 1 



k'^Htl 



^QnU- 36Ö,1,+17 1 Q,U+^SQ^ 



-f- 



360 

n--4 M— 3 n — ^2 n— 1 



}' 



-2-1 I 



}"■ 



24 

^^ 
■*■ 360 

wo das Summenzeichen 2 die in der Mathematik allgemein 

übliche Bedeutung hat. 

Betrachtet man die Störungen 

Ax, Ay, AZy 

wegen des Faktors i* als Gröfsen 2. Ordnung, so gehören auch 

die 8 derselben Ordnung an und der Klammerausdruck 



{ r 



ist von der 4. Ordnung. Daraus folgt aber, dafs: 

1. in erster Annäherung dieser Klammerausdruck un- 
berücksichtigt bleiben kann, 

2. man nur einen Fehler der 6. Ordnung begeht, wenn 
man in zweiter Annäherung die Gröfsen 8 mit Hilfe 
der ersten Näherungswerte von Ax, Ay, Az be- 
stimmt werden. 

Hiemach gestaltet sich das erste Näherungsverfahren 
folgendermafsen; 
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Man bestimme die Gröfsen 

und berechne damit für die aufeinanderfolgenden Intervalle den 
Klammerausdruck 

{■■•■•■}'• ... 

dessen Werte alsdann die während der einzelnen Zeitabschnitte 
eintretenden Störungen der Koordinaten darstellen. 

Zweite Annäherung, 

Um nun zur zweiten Annäherung zu gelangen, sind zunächst 
mit den oskulierendeu Elementen die elliptischen Gröfsen 

3?-?7:r/ U.8. f. 
(deren Bestandteile schon aus der vorigen Rechnung fertig 
vorliegen) für die einzelnen Zeitpunkte zu bestimmen» hierauf 
mit Anwendung der gefundenen Näherungswerte von 

Ax, Ay, A0 
die Ausdrücke 

und aus diesen die Werte des Klammerausdrucks 

{ ■r 

für die verschiedenen Zeiten zu berechnen. Die Störung während 
des n- Zeitabschnitts ist alsdann 

Air«— Aa?„_i=-{ V— | j" 

und zwar, eutsprechend der Genauigkeit der Integration, bis zu 
den Grö&en 5. Ordnung einschl. genau. 

Wird — etwa infolge zu umfangreicher Zeitintervalle — 
eine noch gröfsere Genauigkeit begehrt, so mufs man sowohl 
das allgemeine, in der Gleichung (II) ausgesprochenen Inte- 
grationsverfahren durch Zuziehung der Gröfsen 6. Ordnung 
(vergl. nachher Störungsrechnung der Elemente) ändern, als 
auch die in den Ausdrücken der S bisher vernachlässigten 
Quadrate der Störungen berücksichtigen, überhaupt aber in 
einer wiederholten Rechnung alle bi^ dabip unbeachtet ^ebliebenefl 
Gröfseu 6, Ordnung heranzieheut 
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Rechnungsbeispiel. 

Störungen der Asia durch Jupiter. 
Man hat nach Frischauf (Grundr. d. theor. Astr., S. 93): 
Masse der Asia: m = 

„ Jupiters: m* = 1:1047,88 
Oskulationszeit: 1868, Januar 22 
Länge der einzelnen Zeitintervalle: t^ = 20 Tagen. 

Ungestörte Koordinaten der Asia. 



1868 


X 


y 


i 


logr 


löge 


Jan. 12 

Febr. 1 

„ 21 

März 12 


— 2,5620 
2,6077 
2,6387 
2,6543 


+ 1,0545 
0,8693 
0,6791 
0,4850 


—0,2063 
0,1902 
0,1731 
0,1549 


0,4438 
0,4402 
0,4362 
0,4318 


0,8883 
0,8873 
0,8860 
0,8845 





Für dieselben Zeiten Koordinaten Jupiters: 




x 


y 


j' 


-20».*»«»V 


— 20'*»«y 


-20»*>«V 


+ 4,8100 
4,8448 

4,8748 
4,9001 


-1,2759 
1,1222 
0,9673 
0,8116 


-0,1041 
0,1054 
0,1067 
0,1078 


-44,06 
44,44 
44,83 
45,14 


+ 11,69 

10,30 

8,89 

7,48 


+ 0,95 
0,97 
0,98 
0,99 



Die Zahlenwerte der drei letzten Kolumnen sind in Einheiten 
der 7. Dezimale ausgedrückt, also ist beispielsweise 
-44,06 statt —0,000004406 gesetzt. 



Aus 
0?* — X 



Erste Näherung. 
Berechnung der Werte von Q. 



^3 r** ' ^ Q^ y.13 ' V ^3 



^13 



folgt, wenn man gleich mit ^j-i^w* multipliziert, in Einheiten 
der 7. Dezimale 



Ente Nbherang. Berechnung der Werte von Q. 
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1868 


Q' 


0» 


Q' 


Jan. 12 


— 26,05 


+ 6,00 


+ 1,20 


Febr. 1 


-26.10 


5,40 


1,18 


„ 21 


—26,18 


4,80 


1,14 


März 12 


— 26,19 


4,23 


1,11 



Da die Koordinaten der beiden Planeten nickt für das Osku- 
lationsmoment (22. Januar) und die von da ab gerechneten (je 
um 20 Tage auseinanderliegenden) Zeitpunkte, sondern für die 
Mitten dieser Intervalle gegeben sind, so müssen wir die Rech- 
nung mit einem kleinen Umwege beginnen — der also wegfallt, 
wenn die Koordinaten, entsprechend den oben entwickelten Vor- 
schriften, von vornherein mit der Oskulation anfangen. 

Auch ist die Rechnung, in Ermangelung ausreichenden 
Zahlenmaterials, nur mit einer geringern, jedoch die Gröisen 
4. Ordnung noch mit einschliefsenden Genauigkeit durchfuhrbar 
und demgemäfs auf folgende Formeln zu stützen: 



Axj — Axt ' 


dAx 
•" dt 


1 


A/Sx^ ä/Sxx 


8 


cPAXü 


dt dt 


dt^ 



df^ 



-10 



dt^ 



— 3 



d»A«, 



df^ 



'«1 
24 



■6 



dt^ dt^ 



(9). 



Um zunächst die Störungen und ihre Geschwindigkeiten für 

den 1. Februar 

zu bestimmen, was aus dem angegebenen Grunde auf besonderem 
Wege zu geschehen hat, bedienen wir uns der Gleichungen: 



dt 



dAxi^ _ ti 
dt' 24 

"" dt Y 



2^!^ + 11 ^'^^* 



d^Axo 



dt^ 



dt^ 



dt^ 



\ d^Axo . d^Ax^ d'^Ax^ 



384 



Bezieht sich in diesen (allgemeinen) Gleichungen der Index 1 
auf den 1. Februar, so geht der Index 0,5 auf den Oskulations- 
zeitpunkt zurück, i. Da nun 
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dAüOQji 



dt 
Axüfi - 0, 

so ergiebt sich, weinii man für ,,., hier seinen Näherungswert 

Q^ setzt; wob^i man sich jedoch zu erinnern hat, dafs in der 
obigen Tabelle die Q bereits mit ^1^«= 400 multipliziert sind: 

~dr ^ 4§ö^-2.26,05-11.36,10-f.26,l8)^-0,65 

Aa:^-— 10.0,65 — -g|^(—5-26,05-4ß-?6,10-h3. 26,18)— 3,2 

(Störung am l. Febraar). 
Damit folgt aus den Gleichungen (9): 

—^ 0,65^^(^?6,06+6.26,l0+3.26,18)--2 

Aä?3-~3,2-13 — ^(-36,05 + 10.26,10+3.26,18)"^ — 29,3 

(Störung am 21. Februar), 
ferner aus denselben Gleichungen: 

Z^:c3«-29,3-40 — 2^(— 26,10+10.26,18+3.26,19) 82,4 

(Störung afn 12. März). 
Mit Hilfe der korrespondierenden Gleichungen für die 2. und 
3. Koordinate findet man ebenso aus den Qv und Q^: 

dt "'^* 

Ayi = + 0,7 (Störung am 1. Februar) 

Aj(2 = + 6,1 (Störung am 21, Februar) 

A;8?i = +0,15 (Stprung am 1. Februar) 
£^02 ^ + 1,3 (Störung am 21. Februar). 
Eine Vergleichung dieiser Werte mit den Ergebpissen von 
Frischauf 

Ax,^- 3,3(-- 3,2); Ay, -+0,8(0,7); A^i -+0,2(0,15) 
Aa;i--29,4(-29,3); Ay^'^« +6,1(6,1); Aier,- + 1,8(1,3) 
Aa:3--82,2(-82,4) 
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zeigt in sofern völlige Übereinstimmung, als die kleinen Ver- 
schiedenheiten in den Dezimalen auf die verschiedene Axt der 
Abrundung dieser Dezimalen zurückzuführen siad« 
Um nunmehr zur 

zweiten Armakerung 
überzugehen, müfsten die Werte von 

bestimmt werden. Dieselben sind jedoch an dieser Stelle noch 
so geringfügig, dafs man von ihnen absehen kann und bei Be- 
rücksichtigung nur einer Dezimale in den Störungswerten auch 
absehen mufs. So ist beispielsweise der gröfste Wert der Si 

nämlich 

fifg-t^^cc— 0,002, 

was man mit Anwendung der dem 12. März entsprechenden 
Werte: 

r2«7,3047; r^« 144,21; a?2«7,0454; y2_o,2352; j92 = o,0240 
t/x 1,2873; yis 0,0751; Ä:^«-h0,4lll 

leicht findet, wo zur Bequemlichkeit die lateinischen statt der 
griechischen, auf die oskulierende Ellipse sich beziehende Buch- 
staben genommen sind. 

Das Vorstehende wird hinreichen, von diesen Rechnungen 
eine klare Vorstellung zu geben. Wendet man, wie es hier ge- 
schehen, die Reihen nur zur Berechnung der Störungßn (nicht 
der gestörten Koordinaten selbst) ^n, so dürften die Resultate, 
selbst bei noch gröfseren Zwischenzeiten, stets befriedigend aus- 
fallen, insbesondere wenn auch die höheren Dififerentialquotienten 
der Störungen zu Hilfe genomipen werden. 

StSrungen der Elemente. 

Erheblich einfacher gestaltet sich das Verfahren in seiner 
Anwendung auf die Störung der Elemente, da hier das all- 
mälige Anwachsen der Störungen einflufslos bleibt. 

Es sei q irgend eins dieser gestörten Elemente, so kann 
man entwickeln: 
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WO der Index Null sich wiederum auf den Anfang der Zeit t 
bezieht. In diesem Falle stellt ^o den vollen elliptischen Anteil 
von q dar, während alle folgenden Glieder ausschliefslich der 
störenden Kraft ihr Dasein verdanken. 

Führt man nach der früher gegebenen Vorschrift den zweiten 
und dritten Diflferentialquotienten von ^o si^if ^^^ ersten zurück 
und nennt, wie bisher, das Zeitintervall ^i, so ergiebt sich: 



?i - ?o+ 



dq^ 



dt 



4.. 



wenn 



dq- 



ti' 1 


dqt dq- 


i 






2 2^1 


dt dt 








^ 6 tt^ 


dqi dq-i 
[dt ^ dt 


9 dqo 
-^ dt \ 



dt 



dem Intervalle vor der Oskulation entspricht. Durch 



Zusammenziehung kommt: 



Qi-Qo-^ 



oder allgemein: 



dqp __ dq. 
dt 



dt ^"^ dt 



qn = qn-i + 



12 



dq^ 



n-l 



dq 



n-l 



<-5% 



dt , dt dt 

Diese Gleichung liefert das gestörte Element g„ am Schlüsse 
des w-ten Intervalls nach der Oskulation aus dem gestörten 
Werte desselben Elements am Schlüsse des (w— l)-ten Inter- 
valls. Um 

qn direkt aus ^o 

zu erhalten, hat man nur qn-i durch qn-29 Ö'n-g durch ^»-3 u. s. f 
auszudrücken und in die letzte Gleichung zu substituieren. Es 
wird dann gefunden: 



Ö'o- 



äK 



dqo , dqi 



dt dt 



+ 



dqn- 



dt 



•-) 



also beispielsweise 
/ dq-i , dqo 



\ dt 
l-5(^ 



dq^i j_ äqo ^ dqi 



dt 
dq^ 



dt ^ dt 



dt 
•• + 



dqn\ 
dt I 



dqn-z \ 
dt ) 



dqi dq^ dggX 
dt dt dt ) 



äqx _L dq. 



dt ' dt dt dt 



')'^''\df^~¥^'W;'^'W) 



l)ie St(>rung irgend eines Eiemeniefl q. 
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Iq Gleichung (A) wird uns ein ganz allgemeines Verfahren 
geboten, um 

die StSrung irgend eines Elementes q 

aus den verschiedenzeitigen Geschwindigkeiten dieser Störung zu 
entwickeln, selbstredend in Verbindung mit den Werten q^ des 
Elements zur Zeit der Oskulation. 

Man würde hiernach, wenn es sich etwa um die Störung 
der Flächengeschwindigkeit in der X — Y-Ebene handelte, nur 

ZU setzen, sodann für die Zeiten 

t—i > »0 > *t > ^2 • • • • 
diesen Ausdruck zu berechnen und die erhaltenen Werte in die 
Gleichung (Ä) einzuführen haben. An Einüachheit läist das 
durch (Ä) angedeutete Verfahren geviTs nichts zu wünschen 
übrig; aber auch an Genauigkeit dürfte es wohl den zu steUen- 
den Anforderungen meist genügen. 

Übrigens hindert Nichts, dem Verfahren jeden beliebigen 
Grad von Strenge zu Terschaffen. Man mufs in diesem Falle 
nur mehr als drei aufeinanderfolgende Geschwindigkeiten heran- 
ziehen. Will man z. B. auch noch die Gröisen fünfter Ordnung 
in Betracht ziehen, so setze man: 



dt 
dt 
dt 



dt 
dq^ 
dt 
dqi 
dt 



£q^ 2tt . d*q, 
dt^ ' 1 
d^Qi t^ 



dt* 
d*q. 



dt^ 
d*q. 
dt^ 



1 



dt* 
d% 



2%^ 

2! 


d*q^ 
dt'' 


21 


dt'' 


h' 


, d*ff. 



1 



dt* 



2! 



2%« 

TT"*' 

4;+ 



Durch Elimination von je zwei der höheren Differentialquotienten 
erhält man sodann: 



dt' 

d% 

dt* 

d*q^ 

dt* 



i_ 
1 



dt dt 



dq. 



d«. 



dqj 

dt 



— 2 



dgj 
dt 



dt 



dt 






dqi 



dqt 



dt '^ dt ^^ 



dt 



dqt 
dt 



(b). 
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Werden diese "Werte in die Reihe 



«3 = «2 + 



dt 



1 ■*■ df^ 



eingeführt, so ergiebt sich: 



2s - 3» + ü 



h 



24 
oder allgemein: 

q» = qtt-x + 24 












dg-j 



4! ^ 



^-^19^ + 9 



dt dt dt 



9 



dt 



(So)- 



Bezeichnet Qq wiederum den Zustand des gestörten Elements im 
Momente der Oskulation, so findet man hiernach dessen Wert 
am Ende des w-ten Zeitintervalls: 



3o + 



-( 



k 

24 



fdq- 



dt 



^ dq-\ ^ dqo ^ 



dqn-^s 



dq-\ dqo 
dt dt 



dt dt 
dqn- 



dt 



dt 



')-"( 



<^go 



i-^+ 



dt ' dt 



dq, 



W-i 



dt 



+ 9 



(dqt_ 
\dt 



4-* + 



dt 



oder, bei Anwendung des bekannten Summenzeichens 

(r^n — 9 



' r-— 2 r- — 1 



r-i) 



r=l 



eine Gleichung, welche den Wert des Integrals um den Betrag 
einer Gröfsenordnung genauer liefert als (Ä), 

Bei Berücksichtigung der Gröfsen sechster Ordnung be- 
darf es noch eines weiteren Differentialquotienten -J^ . Das näm- 
liche Verfahren führt, dann auf die Gleichung: 



ri~2 



r^—l 



rrn+l 



Es ist einleuchtend, wie man auf diese Art weiter gehen 
kann, was jedoch nur in solchen Fällen sich als notwendig er- 
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weisen dürfte, in welchen das gewählte Zeitintervall ^, eine un- 
gewöhnliche Gröfse besitzt. - Auch in den mustergiltigen Kech- 
nungea Enokes beschränkt sich die mechanische Quadratur 
auf eine Berücksichtigung der vierten Differenzen derStörungs- 
geschwindigkriten , also auf eine Genauigkeit, welche unserer 

Gleichung {B) entspricht, in welcher gleichfalls noch -r^ in 
Rechnung gebracht ist. 

Was die Koefizienten der Gleichungen {A), (B) und (C) 
betrifft, so finden, wie leicht erklärlich ist, die Beziehungen statt: 

8 — 1+5-12 ] 

1—5+19+9 — 24 ) — den Nennern zu t^, 

11 — 74+456+346 — 19-720 ) . 
Sollen die Störungen für andere Zeiten als gerade die Endpuqkte 
der Intervalle bestimmt werden^ so lassen sich die erforderlichen 
Rechnungsvorschriften ganz auf demselben Wege herleiten. 

Für die Mitten der Zeitabschnitte erhält man z. B., wenn 
man vom Schlufspunkte des betreffenden Abschnitts rückwärts 
geht (man könnte ebenso gut auch vom Anfangspunkte aus sich 
vorwärts bewegen) und die der Gleichung (C) zukommende 
Genauigkeit anwenden will: 



9i^-22-^7jm 



83%-832i^-3442# 
at dt dt 

+ ^ dt ^^ dt 



■(D) 

Für einen beliebigen Zeitpunkt innerhalb dnes Intervalls 
folgt, wenn man diesmal vom Anfangspunkte ausgeht, ferner mit 

die seit letzterem verflossene Zeit bezeichnet, aulserdem die 
Gleichung {B) zu Grande gelegt Mrird: 

wo -3^9 U.S. f. vermöge der Gleichungen (ft) zu ersetzen sind. 

Die Ausführung mehrfftcher Integrale erfolgt nur durch 
eine wiederholte Anwendung der vorßtehenden Regeln. So wird 
das in der Störung der mittleren Bewegung (vergl. S. 48) auf- 
tretende Integral 
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Mm 



dt 



gefunden, indem man aus den gegebenen Geschwindigkeiten 

dn 

dt 
mit Hilfe der Gleichung (J5) zunächst für die einzelnen Zeit- 
punkte das Integral 

herstellt, sodann diese Integrale als neue Geschwindigkeiten -^ 
auffafst, also 

setzt und mit diesen -^ ebenso verfährt wie zuvor mit den -jj • 
dt dt 



Kechnungsbeispiel. 

Die Störungen des Perikels der Vesta durch Jupiter. 
Nach Encke (vergl. Berl. Jahrb. f. 1837, S. 270) hat man 
(Zeit der Oskulation: Januar 0,1810)» 



Ol" Par. Zt. 


42^ 


Ol» Par. Zt. 


dq 


1809, IX, 17 


30"550 


1810, Vm, 19 


93"259 


X.29 


56,829 


IX, 30 


85,839 


XII, 10 


76,602 


XI, 11 


77,461 


1810. 1,21 


90,348 


xn,23 


68,876 


III, 4 


98,589 


1811, II, 3 


60,822 


IV, 15 


102,308 


III, 17 


54.004 


V,27 


102,193 


IV, 28 


48,982 


VII, 8 


98,947 







Die Intervalle umfassen hier 42 Tage, also 
^1-42 . 
und mit dieser Zahl sind bereits sämtliche Werte der Geschwindig- 

keit -^ der Perihelstörung multipliziert. 
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Da die Oskulation in die Mitte eines Intervalls fällt — wie 
es von dem Enck eschen Verfahren verlangt wird — nämlich 

zwischen XII, 10 und 1,21 
liegt, so wollen wir zunächst mit Hilfe unserer Gleichung (D) 
die Störung 

von 1810 Jan. bis Jan. 21 

berechnen, indem wir von letzterem Zeitpunkte an rückwärts 
gehen. Man erhält: 



g'-%=ffo + 



7680 



83 %1^ 832%- -8442% 
dt dt dt 



+ 408.%-57# 
dt dt 



wo der Index dem 21. Januar entspricht 

Durch Einführung der tabellarischen Werte folgt: 

g_^= q^ + J^ [83 . 56,829-832 • 76,602-3442 • 90,348+408 • 98,589 

-57 . 102,308] = ^0 — 43",698. 
Weil q-v^ in die Oskulation fällt, so stellt offenbar 

g,~g_i,, = 43",698 
die gesuchte bis zum 21. Januar eingetretene Störung dar. 

Wir mufsten diese Rechnung vorausschicken, um vor Allem 
einen festen Punkt zu erhalten, an welchen unsere Gleichungen 
(J.), {B) oder (C) unmittelbar angeschlossen werden können. 
Denn diese Gleichungen beziehen sich auf die Geschwindigkeits- 
worte, welche am Anfange oder Ende eines Intervalls stattfinden 
während in der Enck eschen Tabelle die Werte für die Mitten 
der Intervalle gegeben sind. 

Um nun zweitens die Störung 

vom 21. Januar 1810 bis 3. Februar 1811 
zu bestimmen, können wir uns einer der drei Gleichungen {A\ 
(J5) oder (C) bedienen. Wir wählen die letztere und erhalten 
zuerst, wenn man berücksichtigt, dafs jetzt 

-^ dem 21. Januar 1810 
dt 

^ dem 3. Februar 1811, mithin 
at 

Israel-Holtzwart. 7 
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dt 



dem 29. Oktober 1809 



dq- 



i dem 10. Dezember 1809 



dt 

U. 8. f. 

zukommt, aus der Tabelle: 

r=n— 3 



422^ = 42 

r-— 2 



dq^2 ^ dq^i ^ dqp ^ dq^ ^ ^ dqn-si 
dt dt dt dt dt 



dt 

-= 56,829^+ 76,602 + 90,348 + • • •+ 85^839 « 804,914 



42 2 '-t- «25.546 

r =n— 1 

422-^-817,820 

r-ü 

42 2-^-788,294 

und dann weiter nach Gleichung (C): 

- 19 • 743,709] ]qn-qo-+ 804",595 
als Störungsbetrag vom 21, Januar 1810 bis 3. Februar 1811. 
Oben war die Störung bis zum 21. Januar 1810 
43",698 
gefunden worden, womit sich als Gesamtwert der Störung von 
der Oskulation bis zum 3. Februar 1811 

43",698+804",595 - 848",293 
ergiebt. Vergleichen wir dies mit dem Resultate der Encke- 
schen Rechnung 

848",289, 

so zeigt sich eine vollständige Übereinstimmung. Denn der 

kleinen Differenz von 

0",004 
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kann bei einem Zeitraum von mehr als 400 Tagen eine reale 
Bedeutung nicht beigelegt werden. 

Um schliefslich auch eine Anwendung der Gleichung {B^ 
zu geben — welche, wie schon hervorgehoben, der Genauig- 
keit des Encke sehen Verfahrens wohl am nächsten steht — 
mag noch die Störung vom 3. Februar 1811 bis 17. März 1811 
berechnet werden. Es findet sich (wenn man, um sämtliche 
Multiplikationen sofort an der Tabelle ausführen zu können, den 
Koeffizienten 19 = 20-1 setzt) : 
2„ - g„_i « ^jV [(77,461+20. 60,8224-9.54,004;~(5. 68,876+60,822)] 

= ix [(1719,115 ~ 344,380) = 57",281. 
Nach Encke (S. 271) beläuft sich dieselbe auf 
905",573 - 848",289 - 57",284. 
Auch diese Resultate können als übereinstimmend angesehen 
werden. 

Znsatz. Das eben entwickelte Verfahren der empirischen 
Quadratur ist lediglich aus der Eulerschen (Maclaurin sehen) 
Reihe entnommen, während die sonst üblichen Methoden auf die 
Interpolationsrechnung gegründet zu werden pflegen. Im ersten 
Falle schliefst sich die Rechnung unmittelbar an die gegebenen 
DiflFerentialquotienten zweiter Ordnung an, und man ist dabei 
der im zweiten Falle geforderten, immerhin lästigen Bildung der 
höheren Differenzen überhoben. Es liegt hierin eine Bequemlich- 
keit, die um so weniger abzuweisen ist, als die Summenbildung 
der Differentialquotienten in beiden Fällen verlangt wird. 

Dafs aufserdem das erste Verfahren keineswegs an gleiche 
Intervalle gebunden ist, vielmehr ebenso leicht auch für un- 
gleiche Zwischenzeiten entwickelt werden kann, ist ohne 
Weiteres einleuchtend, findet aber auch schon im Anfange unserer 
Darstellung (Seite 79) seine Bestätigung, insofern dort ausdrücklich 
die Ungleichheit der Zeiträume angenommen wurde. 

Übersehen wir nun — am Schlüsse unserer Betrachtungen 
über Berechnung spezieller Störungen angelangt — nochmals 
das ganze Gebiet, so teilt sich dasselbe, wie folgt: 

1. Berechnung der Koordinaten und Geschwindigkeiten 

dx dy dz 
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des gestörten Planeten, sowie der Koordinaten 

^1^ yo ^1^ 
des störenden Planeten für die Zeitpunkte 

und zwar unter Anwendung der zur Zeit der Oskulation 
(^j ^ 0) stattfindenden Elemente, 

2. HerstelluDg der Ausdrücke für die Geschwindigkeiten der 
Störungen aus den sub 1 gefundenen Gröfsen nach Mafs- 
gabe des Abschnitts D, 

3. Bestimmung der gestörten Elemente selbst mit Zuziehung 
der Gleichung (Ä). 



F. 
Mathematische Form der Yariationsmethode. 



Wir haben im Verlaufe unserer bisherigen Untersuchung 
lediglich von der dynamischen Form des Variationsprinzips Ge- 
brauch gemacht. Zur Erreichung einer gewissen Vollständigkeit 
ist es jedoch notwendig, auch die mathematische Form, wenig- 
stens ihrem Wesen nach, hier noch kurz zu erörtern. Die Be- 
handlung eines bestimmten Falls dürfte zu diesem Zwecke am 
geeignetsten sein. 

Stellen wir uns also die sehr allgemeine Aufgabe: von der 
geradlinigen, durch keine äufseren Kräfte beeinflufsten Bewegung 
fortzuschreiten zu der von den beliebigen Komponenten 

PcosQX und Pcos^r 
erzeugten Bahn. 

Im Sinne der mathematischen Variationsmethode müfste 
dieselbe folgendermafsen formuliert werden: 

Es sind gegeben die Differentialgleichungen 

dt^ ( 



sowie ihre Integrale: 

dx . 
dt-'' 

?-« 
dt 






(11) 
•••(III). 
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WO a, h, a, ß Konstante bedeuten, nämlich die Werte der Koor- 
dinaten und ihrer Geschwindigkeiten bei Beginn der Zeit t 
Man sucht die Integrale der Differentialgleichungen 



dp 
d^y 
dt^ 



= Pco8QX 
^ PcosQY 



.(!«). 



Auflösung 

Mit dem Eingreifen der Kräfte P cos QX und P cos Q T 
können die Gröfsen a, b, a, ß nicht länger mehr als beständig 
angesehen werden. Man mufs vielmehr setzen: 



dx . / 



dt 
d^x 
df 



da db 
di'^di 



dß\ 
dtj] 

^'dt'^dtylt'^'i 
d^^da_ d( da dß\ 
dt^ dt'^dt\ dt'^di) 



(II«) 



)!.... 



(Illa). 



Die beiden letzten Ausdrücke sind nun jedenfalls, den Bedingungen 
der Aufgabe gemäfs, gleich 

Pcos QX, bezw. Pcos QY 
anzunehmen. 

Allein damit haben wir erst zwei Bedingungsgleichungen 
zur Bestimmung der vier Gröfsen 

da db da dß 
dt dt dt dt 
so dafs wir noch zwei weitere Bedingungen den Gleichungen (IIa) 
zu entnehmen haben. 

Man könnte über diese zunächst ganz willkürlich verfugen 
und sich in diesem Falle lediglich von dem Gesichtspunkte der 
gröfstmöglichen Einfachheit leiten lassen. Diese würde aber 
offenbar erreicht, wenn man die eingeschalteten Glieder der 
(IIa) gleich Null setzte, da diese auch den Fortfall der Klammer- 
ausdrücke in den (Illa) nach sich zöge. Es ergäben sich dann 
die Bestimmungsgleichungen: 
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da 
^ dt 


db 
^dt' 


= 


^ dt 


^dt 


= 


da 
dt 


= Pcos 


QX 


da 

dt 


= PC08 


QY 



(IV). 



Erwägt man nun aber, dafs in den Aufgaben der Astroaiechanik 
die Gleichungen (III) stets den Anfangszustand der nachfolgenden 
Bewegung darstellen, dafs, mit anderen Worten, die Form 
dieser Gleichungen auch für die künftige Zeit fortbestehen soll, 
so darf man nicht nur, sondern mufs in jedem Zeitmomente 

dt dt 

Ut^dt ^ 
wählen. Denn dadurch reduzieren sich die Gleichungen (IIa) 
auch bei Beginn der Zeit t, wie es sein mufs, auf: 

dx 

dy 

Wir gelangen also auch durch diese Überlegung zu denselben 
Bestimmungsgleichungen (IV). Substituiert man nun die hieraus 
folgenden Werte von 

da db da dß 

dt dt dt dt 
in das System (III), so erhält man für die durch die Kräfte 
Pcos QX und Pcos QY gestörten Koordinaten die Gleichungen: 

^ -^ [«0 +/Pcos QX.dt] t-h[ho—ft. Poos QXdt] 1 
y - [cto-hfPcos QY,dt]t-h [ßo-ftPooB QYdt] i 
Eine allgemeine Ausführung der Integrale 

/^Pcos QXdt und ftPcos QYdt 
würde die Kenntnis der Kräfte als Funktionen der Zeit, also 
die Lösung der Aufgabe, bereits voraussetzen; sie ist deshalb 



(V). 
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nicht einmal in dem einfachen Falle der elliptischen Bewegung 
möglich, wo 

Pcos(2X=--^f^ 

Will man die Gleichungen (V) dessenungeachtet benutzen, so 
bleibt nichts übrig als ein empirisches Verfahren — nach Ana- 
logie des seither in der Störungsrechnung gebrauchten -^ zu 
Hilfe zu nehmen, durch welches aus einzelnen numiörischen 
Werten der Kräfte diese selbst als empirische Funktionen 
der Zeit erhalten werden. 

Es verdient besondere Erwähnung, dafs die beiden Integrale 
in jeder der Gleichungen (V) sich in ein einziges zusammenziehen 
lassen. Denn nach der Regel von der partiellen Integration 
hat man: 

t .fP cos QX dt—ft P cos QX dt SdtfP cos QX dt 
tfPcos QYdt-ftPcos QY dt -^fdtfP cos QYdt 
Hiernach sind die Gleichungen (V) auch durch die folgenden 
ersetzbar: 

x^\-\-aQt-\'fdtfPco^QXdt\ ,„ . 

y^ßo + ao^-i-fdffPcosQYdti'"^ 
Wir kommen so wieder zu denselben Gleichungen, welche 
im Abschnitte E den Ausgangspunkt unserer Betrachtungen 
bildeten) und diese Thatsache läfst die Natur der Variations- 
methode — als einer blofsen Reduktionsmethode — klar erkennen. 
Man kann aber den Gleichungen (Va) aber auch noch die 
folgende Form geben: 

^ -^/[«o +-/PCOS QX dt] dti-bo^ . . 

y ==/k+/Pcos QYdt'ldt^ß^ y-'-y. Ph 

welche wohl als der natürlichste Ausdruck des Problems betrachtet 
werden dürfen, insofern sie die Störung durch 



lediglich auf 
übertragen. 



Pcos Q^-^ und Pcos QY^ jß- 
«,«^unda,=.^ 
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Überhaupt aber lehren die verschiedenen Formen 
(V), Va) und (Yß\ 
dafs derselbe Zweck auf mannigfache Weise erreicht werden kann. 
In den Gleichungen (V) unterliegen sämtliche Oskulations- 
konstanten a h a ß 

der Variation. Die Gleichungen (Vor) gestatten die Auffassung, 
dafs nur 6^ und /J^, also die Anfangswerte der Koordinaten, 
Variationen erleiden, während in den Gleichungen (V/5) um- 
gekehrt nur^o undofo^ also die Oskulations werte der Geschwindig- 
keiten, gestört erscheinen. 

Im Vorstehenden haben wir die mathematische Form der 
Variationsmethode nur auf einen bestimmten Fall angewendet, 
nämlich zur Ableitung der gestörten aus der oskulierenden 
geradlinigen Bewegung. Dennoch wird der Leser den all- 
gemeinen Grundgedanken unschwer herausfinden. Wir dürfen 
uns daher auf folgende Andeutungen beschränken: 

Angenommen, man kenne die Koordinaten der oskulierenden 
(etwa der elliptischen) Bewegung als Funktion gewisser Konstanten 

a, €, fj, (Oy Sl^ i 
und der freien Zeit t, so können wir diesen Zusammenhang durch 
folgende Gleichungen 

x=^f(a, €, jii, CO, Sl, i, t) 
y = f/)(a, e, f,i, w, ft, i, t) 
g'^ipia, €y fj, (0, ü, i, t) 
zum Ausdrucke bringen. 

Differenziert man nun diese Gleichungen zweimal nachein- 
ander, wobei auch die bisherigen Konstanten als veränderlich 
betrachtet werden, so ergeben sich 6 Gleichungen, aus denen 
die 6 Bestimmungsgleichungen für die Variation der 6 Konstanten 
dadurch gewonnen werden, dafs man zunächst die Aggregate der 
ersten Differentialquotienten (nach Abzug der auf die freie 
Zeit t sich beziehenden) 

= Null, 

sodann, unter Berücksichtigung dieser 3 ersten Bedingungs- 
gleichungen, die zweiten Differentialquotienten 

*= den störenden Kräften 
setzt. 
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Das Prinzip ist klar und einlach, die Durohfillirung der 
Rechnung aber in den meisten Fällen — z, B. bei dem Übergange 
von den elliptischen zu den gestörten Elementen einer Planeten- 
bahn — überaus mühevoll. Dieselbe wird besonders erschwert 
durch die hier stets notwendig werdenden Eliminationen, obwohl 
die Unbekannten 

da d€ 

dt' dt'"' 
nur linear in den Bestimmungsgleichungen auftreten. Grerade in 
der Vermeidung dieser Ehminationen liegt der wesentlichste Vor- 
zug der dynamischen Form der Variationsmethode. 



G. 

Berechnung der Planetenbahnen 
aus Beobachtungen. 



Die allgemeinen Bedingungen des Problems. 

Der eben abgeschlossenen Untersuchung über die Bewe- 
gungs-Störungen der Himmelskörper lag die Kenntnis der unge- 
störten elliptischen Bahn als wesentliche Voraussetzung zu 
Grunde. Sowohl dieser Umstand, als die, wie wir bald sehen 
werden, mit den vorhergehenden Entwickelungen obwaltende 
methodische Verwandtschaft veränlafst uns, auch dem 

Probleme der elliptischen Bahnberechnung 
eine etwas eingehendere Besprechung zu widmen. Ohne an 
dieser Stelle schon eine strenge Definition dieses Problems zu 
geben, sei nur im allgemeinen bemerkt, dafs dasselbe zum 
Gegenstande hat, die als elliptisch angenommene Bahn eines 
meist bis dahin noch gänzlich unbekannten Himmelskörpers aus 
geeigneten, stets nur einen kurzen Zeitraum umfassenden Beob- 
achtungen zu bestimmen. 

Bedeuten 

^z, y%y H 

die zu erforschenden rechtwinkligen heliozentrischen Koordinaten 
dreier Planetenörter in Bezug auf die Ekliptik, sowie 

tx und ^2 
die Zeitabschnitte, welche beim Übergange des Plsi^neten vom 
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ersten zum zweiten, bezw. vom zweiten zum dritten Orte ver- 
fliefsen, so hat man nach der Eul er sehen Reihe: 



Xi = x^ 



Vi - Vi 



^1 « ^2 — 



X^ =^ Xo ^ 



2/3 == ^3 + 



^3 + 



dt 

d02 

dt 

dx2 

'df 

dt 

d02 
dt 



1 



+ 






'dF 

d^ 
dt^ 
d^ 
df' 
d^ 
df^ 

dt^ 
d% 
dt^ 



2! 

A!_ 

2! 

A!_ 

2! 

21 ^ 
21 ^ 
2! ^ 



Qf Xn ^1 



(I)- 



Da die Differentialqaotienten der zweiten Ordnung nach dem 
Gravitationsgesetze durch die Koordinaten ausgedrückt werden 
können, indem 

d^Xi k^Xi 



dp 
d% 



k% 



U. 8. f., 



dt* r,' 

alle höhere Quotienten aber auf die Koordinaten und die beiden 
ersten Quotienten reduzierbar sind, da beispielsweise 



d^Xj 
dt^ 



= k^ 



oder*) 
d^Xj 



3x 



*» 



0«2/2 / "^ 



dx^ 
~df 



dr^ _ 
dt 



dx2 

~dt 



02 ^\ L. ^ 

r. dt 1 ^2^ dt 



^2 dt 

u. s. f., 
so enthält das Gleichungssystem (I) nur zwölf Unbekannte, 
nämlich die neun Koordinaten und die drei Dififerentialquotienten 
der ersten Ordnung. 



*) mit Berücksichtigung der Gleichung 



dr^^x^dx^ ^ y^dy^ ^ z, dz^ 



dt 



dt 



dt 



dt 



folgt. 
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Liefert also die Beobachtung eine Reihe Ton Gröfsen, 
mit deren Hilfe sechs weitere Beziehungen zwischen den Unbe- 
kannten aufgestellt werden können, so hat man für die zwölf 
Unbekannten eine genügende Anzahl von Bestimmungsgleichungen 
und unser Problem ist, wenigstens sachlich, gelöst. 
Der Ort der Erde, d. h. 

ihre heliozentrische Länge L und 
ihre Entfernung von der Sonne B 
ist aus der bekannten Theorie der Erdbewegung fiir die drei, 
durch die Intervalle 

t^ und ^2 
getrennten Beobachtungsmon^nte leicht und genau bestimmbar. 
Wir dürfen hiernach 

als gegeben betrachten. 

Nehmen wir nun ferner an, es seien die den drei Erdörtern 
entsprechenden geozentrischen Längen und Breiten des 
Planeten, nämlich 

^u ßi 

A2, ß'2 
^3» ßz 

durch Beobachtung gefunden worden, so haben wir eine ge- 
nügende Anzahl Data, um die fehlenden sechs Bestimmungs- 
gleichungen bilden zu können. 
Bezeichnet man nämlich mit 

Q 

den (unbekannten) Abstand der Erde von dem Planeten, so hat 

man, wie wir aus der analytischen Geometrie wissen, aufserdem 

durch eine leicht zu entwerfende Figur bestätigt wird (vergl. 

u. A. des Verf. Elem. der Astron., IL. S. 100): 

x^-^Q^ sinA j 

Vi "= ^1 ^os ßi sin A| + iJ| sinZj > • • • • (") 

Zx = Qi COS ßi cos A^ -\' ü| cosXi j 

und ähnliche Gleichungen für die 6 anderen Koordinaten. Dabei 

ist also die X-Achse senkrecht zur Ekliptikebene angenommen. 



HO 
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(11) 



Eliminiert man die (>, so folgt: 

yi = iC| . cot ß^ sin A| -{- i?i sin L^ 

^1 = Xj . cot ßi cos A| ■+■ JB| cosZri 

2/2 = % . cot ßo sin Ao + Ro sin Zo 

jSo^^oCo' cot /^2 cos A2 -I- Ro cos Zo 

^3 = ^3 • cot /?3 sin A3 4- R^ sin X3 

Z3 «= 0^3 . cot ß^ cos A3 -h JB3 cos ig 
und damit sind wir im Besitze der gesuchten sechs weiteren 
Relationen, welche verbunden mit dem Systeme (I) das Problem 
lösen. 

Wir erkennen hieraus — was mit Hilfe allgemeiner end- 
licher Gleichungen erst durch eine sehr umständliche Unter- 
suchung möglich ist — dafs auf Grund folgender Data: 

1. der drei Erdörter, 

2. der drei Paare geozentrischer Längen und Breiten des 
Planeten, 

3. der beiden Zwischenzeiten, 

das Bahnproblem ySlUg bestimmt ist. 



Erste näherungsweise Bestimmung der Bahn. 

Beseitigt man aus dem Systeme (I) die ersten Differential- 
quotienten, so ergiebt sich: 



1.2 



dt' 
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"*" 41 dt^ 

■*" 5! W^ ^"^^ 

1.2 dt' '^ 3 ! dt^ 

^1^3 + ^2^1 — (^1 + ^2>2 

tMh+k) d% . tMk^-h'^) d% _ 

1.2 "di^"*" 31 d^"*" ^ ^'' 



Xni) 
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. Wir entnehmen nun diesen Gleichungen — deren Entwicke- 
lungsgesetz ebenso klar wie einfach ist — weiter, dafs man die 
Zwischenzweiten t^ und ^2 zwar vom theoretischen Standpunkte 
ganz beliebig wählen kann, dafs es sich für die Rechnung hin- 
gegen empfiehlt dieselben gleich, oder wenigstens nahezu gleich 
zu wählen, da alsdann die dritten (überhaupt die ungeraden) 
Differentialquotienten yerschwinden^ dafs mau ferner dieselben 
nur von mäfsiger Gröfse (etwa 8 bis 14 Tagen) wählen darf, 
da sonst die höheren Differentialquotienten zu stark ins Gewicht 
fallen und ihre Vernachlässigung bei einer ersten Nälierungs- 
rechnung allzu schädlich wirken kann. Die Bemühungen der 
Astronomen sind deshalb auch stets darauf gerichtet, die Beob- 
achtungen diesen beiden Bedingungen zu unterwerfen. Sind 
dieselben erfüllt, so lehrt die Erfahrung, dafs man stets zu 
brauchbaren Näherungswerten gelangt, wenn man in erster 
Rechnung die rechten Seiten der Gleichungen (III) 

annimmt, mithin bei Beginn der Rechnung die Gleichungen an- 
wendet: 

hx, + t^x, - (t, + t,)[l- -^j x,^0 

t,yB + t,y,-{t,+ t,)l^l-^'^y, = 0\ . . . (in«). 

<i^3 + <2^i - (f. + *2)(l -■ ^) ^2 = 

Der fernere Verlauf der ersten Näherungsrechnung ist nun leicht 
zu überblicken. Setzt man nämlich 

^2'-(^2'-Hy2'+^2'Ä 
so ist klar, dafs sich mit Zuziehung der 3. und 4. Gleichung des 

Systems (II) dieser Ausdruck sich in eine blofse Funktion von 

a?2 verwandelt und dafs dieselbe Reduktion mit dem t/z und 02 

der (Illa) vorgenommen werden kann, so dafs überhaupt also 

die drei letzten Glieder der Gleichung (Illa) blofs noch die 

Unbekannte X2 enthalten. Werden nun gleichfalls mit Hilfe der 

(II) auch die Koordinaten 
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aus den (Illa) entfernt, so nehmen letztere eine Form an, aus 
welcher auch die beiden Koordinaten 

x^ und Xi 
ohne Weiteres fortgeschafft werden können. Die Finalgleichung 
enthält dann lediglich noch die Unbekannte 

deren Bestimmung schliefslich mit irgend einer empirischen 
Methode — etwa der regula falsi — zu erfolgen hat. Der Eü- 
minationsprozefs gestaltet sich zwar etwas umständlich, aber, 
da er durchängig linearer Natur ist, ohne jede Schwierigkeit. 
Wollte man die Schlufsgleichung rational machen, so würde sie 
in Bezug auf X2 bis zum 8. Grade ansteigen. Übrigens kann 
man — wenn dies vorgezogen wird — sogar mit blofsen linearen 
Gleichungen operieren, indem man zunächst auch die (das Quadrat 
der G aufs sehen Konstanten k enthaltenden) Differentialquotienten 
der zweiten Ordnung aus den Gleichungen (Illa) tbrtläfst, so- 
dann mit den so erhaltenen Werten von x^y ^2» ^2 ^i® Gröfse 
von ^2 berechnet, letzteren Wert in die Gleichungen (III«) ein- 
führt, die Koordinaten X2, y^y H ^^^ neuem berechnet und dies 
solange fortsetzt, bis eine Verbesserung von r^, nicht mehr 
eintritt. 

Selbstverständlich lassen sich statt der rektangulären Koor- 
dinaten 

x^y y^, ^2 

^3> yzy ^3 

vermöge der Gleichungen (a) auch die geozentrischen Entfer- 
nungen (J„ Q2, ps 

des Himmelskörpers einführen. Geschieht dies, so ergiebt die 
Auflösung des Systems (Illa) folgende Bestinmiungsgleichungen: 

1 a^itj + ^2) — a^tj—aih h\{ti + ti) 

1- Qi- ^ ^^"2/2^^ 



£f fl« = . — : — 1 ; — : — ; • —x, r- 






\. .(IV), 
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wenn zur Abkürzung die aus den Beobachtungen folgenden 
Grölsen 

tg ßi sin (A3 — -^2) — tg ßz sin (A3 — A^ ) + tg/?3 sin (Ao — A| ) « a 
[tg ß^ sin (ii — A3) — tg /J3 sin (i, — 2j)] ^^^^^ — a, 

TD 

[tg /?3 sin (Li — A J - tg/i^i sin (A — h)'] "^"^^^ '== *i 
[tg ß, sin (L, - AO - tg Ä sin (L, - A,)] ■ ^' - -- 



acosi6?3 
gesetzt werden und man unter 

die analogen L^ und 1^3 statt L^ in sich aufnehmenden Aus- 
drücke versteht. 

Fügt man zu den Gleichungen (IV) die aus dem Dreiecke 

Sonne - Planet - Erde 
folgende Relation 

ifi- ^ (i2^ + 2iJ2 cos/?2 cos {L^-'h)Qi + -^2^ {ß\ 

so sind jene Gleichungen auf sehr verschiedene Weisen lösbar, 
am einfachsten dadurch, dafs man in {ß) zunächst 

__ &3^i 4- fei ^2 ~ ^'*(^l + ^2) 

setzt und mit dem daraus hergeleiteten Werte von r^ in die 
zweite der Gleichungen (IV) eingeht, sodann den solchergestalt 
verbesserten Wert von q^ von neuem in {ß) substituiert, mit 
dem korrigierten r^ einen abermals genauem Wert von q^ sucht 
und diese Korrektionen zwei- oder dreimal, überhaupt so lange 
wiederholt, bis etwa die drei ersten Dezimalen von ^2 ^ß" 
verändert bleiben. Den gefundenen Wert von q^^ bezw. rj 
betrachtet man als den definitiven 

ersten Näherungswert, 
der dann auch zur Berechnung von q^ und ^3 aus der ersten 
und dritten Gleichung des Systems (IV) verwendet wird. 

Es würde Zeitvergeudung sein, diese erste Näherungs- 
rechnung weiter, d.h. bis zur Erlangung streng genügender 

Israel-Kolt/.wart. 8 
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Zahlenwerte auszuführen , da den Gleichungen (IV) selbst noch 
ein Fehler der dritten Ordnung (d. h. ein Fehler von der Ord- 
nung TcH^) anhaftet, welcher der Regel nach bereits auf der 
3. oder 4. Dezimale zur Geltung kommt. Daher können auch 
unbedenklich bei diesem ersten Teile der Rechnung vier- oder 
fünfstellige Logarithmen gebraucht werden. 

Verbesserung der Zelten wegen Aberration. 

Die Zwischenzeiten 

t^ und ^2 
der drei Beobachtungen bedürfen stets noch* einer kleinen Be- 
richtigung wegen der s. g. Aberration, die man — wenn der 
Himmelsköper, also auch die Entfernung q noch völlig unbe- 
kannt ist — nach der ersten Näherungsrechnung anzubringen 
pflegt. 

Das Licht braucht nämlich zur Zurücklegung des Erdbahn- 
halbmessers 

0,0057 Tage 

zum Durchlaufen der Strecke q (Entfernung des beobachteten 
Himmelskörpers von der Erde) -demnach 

Die Zeitpunkte der Beobachtungen sind mithin um diesen 
Betrag zu vermindern, damit sie den beobachteten Örtern ent- 
sprechen. 

Nennt man T^ und T^ die verbesserten, zwischen den 
drei Beobachtungen liegenden Zeiträume, so erhält man dieselben 
also, nach vorgängiger Berechnung von p^, ^2 '^^^ Qz^ aus den 
Gleichungen 

yi = ^i-((>2-ei). 0,0057^ 
y2 = ^2-(?3-?2).0,0057^ 
und diese korrigierten Zwischenzeiten sind nun den weiteren 
Näherungsrechnungen zu Grunde zu legen. 

Berücksichtigung der Grfffsen 3, und 4, Ordnung. 

Die Gleichungen (IV) sind im Verein mit {ß) zur Gewin- 
nung guter erster Näherungswerte sehr geeignet, vorausgesetzt, 
dafs die Zwischenzeiten eine geringe und möglichst gleiche Gröfse 
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besitzen. In dem Mafse jedoch als diese (immer recht lästige) 
Bedingung unerfüllt bleibt, vermindert sich auch der durch jene 
Gleichungen erreichte Grad von Annäherung. 

In solchen Fällen verdient das folgende Verfahren besondere 
Beachtung, welches von vornherein auch die 

dritten und vierten Differentialquotienten 

demnach die Gröfsen 3. und 4. Ordnung mit völliger Strenge 
berücksichtigt und die Rechnung somit von der postulierten 
Bedingung der Gleichheit und Kleinheit der Zwischenzeiten 
möglichst frei macht. 

Aus unseren früheren Entwickelungen (S. 80) kennen wir 
bereits die bis auf Gröfsen der 5. Ordnung genauen Gleichungen : 
d^o[)2 I L ^d^ ^ d^Xj , . d^] ] 



dt'' tMh + k) 
d'^x^ 2 



dt* tit^it^+t^) 

d^V 
und die ähnlichen Relationen für --^ u. s. f. 

Substituiert man diese Werte in das Hauptsystem (III) 
und setzt abkürzend 

Ti to" ^2^1 rp 

j^ ^ 

^2 H " ^1 ^2 rp 

\^ °^' 

12 ^* 

so erhält man die Gleichungen: 

U. S. f. 

Werden mit Zuziehung der (a) wieder die Entfernungen q 
eingeführt und dann die Gleichungen gelöst, so ergeben sich die 
bis auf Gröfsen der 5. Ordnung strengen Relationen: 

8* 
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«2(^1+^) — a3^i~%^2 . ij+k 



T.-:-T.^^ 






rp k rf "^" 

/ ^1 TT" ^3 "TT 

7- «3 • 
^2 1_ 



2; 



3. 



iL 



Qi - 



^1 + ^2 






77 ^ m ^ 



1 



22 



Ä* 



\. 



T T 

-^2 r -'^1 — 5" 



l-^i^TT 



3. 









-'»TT ~ -'»TT 



1-T, 



*» 






m ^ /T7 ^'' 



(W), 



(«J) 



welche mit den Hilfsgleichungen 

^1^ = Pi^ + 2i?i cos ft cos (Li — A| )ßi + iJiM 

^2^ " C2* + 2JB2 cos ß^ cos (L2 — ^2)(>2 + ^' \ ' • 

rz^-Qz' + 2JB3 cos /?3 cos (ij — -13)^3 + ^3' I 
einer sehr einfachen Auflösung fähig sind. 

Zunächst setzt man in den (d) wieder (wenn noch keine 
Näherungsrechnung vorausgegangen ist): 

^' «2(^1+^2) — Ml — %^ _ ^3^1 -^hh — h(ti + ^2) 

^2 " t^-\- t.y 

cSi±k) — Cjt^ — C^tz 
g,^- -^ , 

berechnet mit den gefundenen Werten von 
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die Entfernungen Q2, Q29 Qz aus den (lUß) genauer, bestimmt 
aus diesen wieder durch Anwendung von (d) verbesserte Werte 
der r, sodann aus (Qlß) abermak genauere Werte der q und 
setzt dies Verfahren fort, bis man streng genügende Zahlen- 
werte gefunden hat — wenigstens, wenn man, nach sorgfältiger 
Prüfung der konkreten Verhäftnisse, glaubt mit einer Berücksich- 
tigung der Gröfsen 4. Ordnung die Rechnung abschliefsen zu 
können, was namentlich bei neuen Bahnberechnungen (nach Gaufs) 
die Regel bildet. Beabsichtigt man hingegen, noch eine weitere 
Näherungsrechnung folgen zu lassen, so wird es genügen, jene 
wechselseitigen Verbesserungen der q und r bis zur genauen 
Ermittelung etwa der 5. Dezimale vorzunehmen. 

Wird die auf Grund der (111^?) und (d) geführte Rechnung 
mit den von der Aberration noch affizierten Zeiten begonnen, 
so mufs man selbstverständlich im geeigneten Moment — nach- 
dem etwa die drei ersten Dezimalen der q bestimmt sind — die 
Intervalle von dem Einflüsse der Aberration befreien, und dann 
erst mit den verbesserten Zeiten die Rechnung vollenden. Aller- 
dings verlangt das vorstehende Verfahren die gleichzeitige Auf- 
lösung dreier Gleichungen. Aber dieser Umstand wird bei der 
einfachen Struktur der Gleichungen (HIß) kaum als eine Er- 
schwerung empfunden werden, zumal man sich am Schlüsse der 
Rechnung auch sofort im Besitze der drei Hauptgröfsen 

Qu Q2^ Qz 
befindet. 

Fortsetzung der Nitherungsrechnnngen durch Einschaltung 

neuer Örter. 

Das eben entwickelte Verfahren, dessen Wesen in einer 
Reduktion der 3. und 4. Differentialquotienten auf diejenigen der 
2. Ordnung besteht, läfst sich nun auch — allerdings mit einem 
erheblich gröfseren Aufwand von Rechnung — auf die noch 
höheren Differentialquotienten ausdehnen. 

Um zunächst dio Quotienten 5. und 6. Ordnung numerisch 
feststellen zu können, würden aufser den aut die drei beob- 
achteten Planetenörter sich beziehende!! Pifferentquotienten 
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2. Ordnung noch zwei weitere solcher Quotienten erforderlich 
sein — und zwar dürften diese sämtliche 5 Quotienten nur noch 
Fehler der ?• Ordnung enthalten, da sie sonst zur Bestimmung 
von Gröfsen 6. Ordnung untauglich wären. 

Angenommen, es seien auch die den Mitten der Zeit- 
intervalle 

tx und ^2 
entsprechenden Quotienten 2. Ordnung 

-dP' ^""^ -dP' 

mit der eben als erforderlich bezeichneten Genauigkeit bekannt, 
so könnten folgende Gleichungen angesetzt werden: 
cPXj d^X2 d^x^ ^ d^a?2 t^ 

'dP dt^ dt^ ' 1'^'"'^ dt^ ' 41 
^^•271,5 d^x2 dh:^ (ti\ , d^x^ ti^ 



ih 



dP dP- dp \2) dp 2*. 4! 

d^^2,5d^Xo d^X2 t^ d^ to* 

dP dP dP '2.r^'" *■ dp *2*.4! 
d^x^ d^Xn d^X2 ^2 , d^arg tj^ 

Hp dP'^ dP •"1+--+ dp '41 

Aus diesen 4 Gleichungen liefsen sich aber die 4 Gröfsen 

d^x^ d*X2 d^X2 d^X2 

'W' dp ' dp ' dp 
sofort numerisch bestimmen. Es bedarf dabei wohl keines be- 
sonderen Hinweises, dafs die früher gefundeneu und benutzten 
Werte von 

elf Xn 1 Cv Xtf 

im gegenwärtigen Falle nicht mehr ausreichen. 

Wir haben also nur noch zu überlegen, auf welchem Wege 
wir uns die beiden Quotienten 

d^Xi 5 d^X2 5 

-dt "'^•^ -dt^ 

verschaffen können. Denn die drei anderen Quotienten 
d-Xj d^X2 d^x^ 
dP ' ~dP' dP~ 
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sind wegen des inhärierenden Faktors Tc^ mit der verlangten 
Genauigkeit bekannt, da wir die Gröfsen 

und mit Hilfe von (of) auch 

^2> y2J H 
an dieser Stelle bereits bis zu den Gröfsen 4. Ordnung inclusive 
vermittelt haben. 

Zunächst ist einleuchtend, dafs wir die Quotienten 
1. Ordnung auf Grund der vorangegangenen Rechnung ohne 
Weiteres mit Hilfe der Vergleichungen 

^^~^^ dt ' l^ dt'' "2 dt' ''6\'^'W'Ä\ 
dvo ti 

dZ', t* . » 

hinschreiben können, nämlich: 

u. s. f. 
Damit haben wir aber Alles, um zwischen und neben den 
ursprünglichen Planeten örtern beliebige andere interpolieren 
zu können — mit der nämlichen Genauigkeit, welche den Koor- 
dinaten jener Örter zur Zeit innewohnt, im vorliegenden Falle 
demnach bis zu Fehlern der 5. Ordnung. Man findet nämlich: 

^1,5 *" *^2 

yifi = y2 — -^-^ 
^v^^ + ^-^-t--- 

u. s. f. 

Es erfordert diese Interpolationsrechnung kaum mehr als 

einfache Additionen bereits bekannter Gröfsen. Die so gefundenen 

Koordinaten der interpolierten Örter enthalten nur noch Fehler 

der fünften Ordnung. Ihre Diflferentialquotieuten der 2. Ordnung 



dXi tt d^Xi 


ti' 


d%. 


h' 


d*x.. t,* 


dt 2-1 ' dt' 


2--2 


dt"" 


2».3! 


' dt* 2*4! 


dy-i ti 
dt 2-1 ' '■' 
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folglich blofs Fehler der 7. Ordnung. Dieselben sind mithin zur 
Bestimmung von Gröfsen 5. und 6. Ordnung — welche hier 
verlangt wird — vollkommen geeignet. Nachdem nun diese 
letztere erfolgt, unterliegt auch die Verbesserung der Koordi- 
naten keinen weiteren Schwierigkeiten Es sei nur noch bemerkt, 
dafs die wirkliche Ausfuhrung dieser Rechnung zunächst eine 
allgemeine Auflösung der Gleichungen (III) bedingt, welche 
ergiebt: 

tA + ^2^ — (^1 + ^2)^2 . sin ( A3 — kj) 

tg ßi cos A3 — tg /? 3 cos kl tg ßi sin Ä^ — tg ßj sin A3 ^^ 

a cos ßz ' if'i cc cos ßo ^ ' 

Wie man auf diesem Wege — durch Einschaltung neuer 
Örter — weiter gehen könnte, ist leicht ersichtlich. Die Inter- 
polation selbst verursacht dabei die geringsten Schwierigkeiten; 
der Mifsstand liegt lediglich in der von Näherung zu Näherung 
sich mehrenden Anzahl der aufzulösenden Gleichungen. Be- 
schränkt man aber die eben von uns betrachtete Methode, 
höhere Differentialqnotienten auf niedere zurückzuführen, 
auf die Reduktion der 3. und 4. Differentialquotienten, wie die- 
selben in den Gleichungen (III^) zum Ausdrucke kommt, so 
giebt es keine andere Methode, welche bis zu dieser Stufe der 
Annäherung sich ihr an Einfachheit vergleichen liefse. 

Um aber eine nach allen Seiten befriedigende Lösung des 
Bahnproblems zu erhalten, wird man die letztgenannte Methode 
noch durch ein Verfahren ergänzen müssen, welche weitere 
Annäherungen gestattet, ohne dafs mit jedem Male die Schwierig- 
keiten sich wesentlich erhöhen. Ehe wir jedoch hierzu über- 
gehen, wollen wir eine kurze Znsammenstellung derjenigen 
Grleichungen vorausschicken, welche für die weitere Entwicklung 
des Bahnproblems von hervorragender Bedeutung sind. Die 
Ableitung dieser Gleichungen können wir dem Löser um so mehr 
überlassen, da dieselbe der Regel nach aus den Polargleichungen 
der Ellipse unter Zuhilfenahme fder Kepl ersehen Gleichung 
(S. 25) und elementarer trigonometrischer Relationen bewirkt 



Zusammenstellung wichtiger Gleichungen. 121 

werden kann. Nur bei denjenigen Relationen, für deren Her- 
leitung diese Mittel nicht ausreichen, werden wir den Weg ihrer 
Begründung kurz angeben. 

Bezeichnungen: 
a = grofse Halbachse der gesuchten Bahn, 
'P = Semiparameter, 
6 «= Excentrizität, 

i = Neigung der Bahn gegen die Grundebene, 
^1» ^2» ^3 ^ wahre Anomalien in der 1., 2. und 3. Beobachtung, 
%, ^2» % ^ excentrische „ „ „ „ „ „ „ „ 

fi^fi^F -^ den doppelten Dreiecksflächen, welche von dem 1. 
und 2., dem 2. und 3. sowie 1. und 3. Radius mit 
den zugehörigen Sehnen gebildet werden, so dafs: 
/i = ^t^2 8in(7>2 — ^i) 
h — ^2^3sin(y3 — <]P2) 
F = rir3sin(f/)3 — yi). 

Gleichungen: 
Aufser den aus dem Früheren (vergl. S. 24 und 25) be- 
kannten Gleichungen: 

r(l + € cos9))«=«p (1) 

r = a(l — ficosw) (2) 

(l+fi)tg*|-(l-e)tg=«| (3) 

u — esinu^nt (4) 

sind insbesondere noch folgende von allgemeiner Bedeutung: 

]/r^ sin ^^ '^ ^^ =- |/ap sin '^—^ .... (5) 

yriz-j sm ^2—^-*- « ^ap sm — -g-^ .... (6) 

rir2sin(y2— yi)==Äl/p ^1 — « i' [% — M| — sin(w2 — Mi)] • (7) 

^ 2 v^ cos 5^^ ^ cos^^-^' ^ 



2 2 



Mi 8-Mi-8in(% -Mi) 4 4.6 . j Mj-Mj 4.6.8 . , «^-Mi , ,„. 

=^ "^ "T TT^ 8in — TZ r ö~£~?7 sm — ;t |- . . . i j^ i 

31*2-— % 3 3.0 2 . 3.5.7 2 ^ ^ 



sm^ 2 



(durch die bekannten Mittel der Analysis zu entwickeln!) 



122 



G. Berechnung der Planetenbahnen aus Beobachtungen. 



fiyz-^hVi-Fy.^-^) (10) 

Es ergeben sich diese für die Lösung des Problems grund- 
legenden Gleichungen am leichtesten auf folgendem Wege. Er- 
wägt man, dafs nach der analytischen Geometrie 

y^Xx — x^yi 

yz^i—^zVi 
die Projektionen der Doppeldreiecke f2,fi,F auf die Grund- 
ebene der xy bedeuten, so hat man, — wenn i die Neigung der 
Bahn gegen diese Ebene — 

h 
Nun ist identisch: 

^1(^3^2 — ^3^2) + ^^3(^2^1 — ^2^0 — ^2(^3^1 — ^3^1) = 0. 

Wird mit cos i dividiert, so erhält man die erste der 
Gleichungen (10). Ebenso ergeben sich dann auch die zwei 
folgenden auf die beiden anderen Koordinatenebenen sich be- 
ziehenden Relationen. 

Aus den Polargleichungen der Ellipse für die 3 Örter 
folgt ferner: 

"-- — W--^ — —■■}"'' 

Durch eine etwas mühevolle Rechnung (indem man zunächst 
die Koordinaten x,y,z durch die excentrische Anomalie und 
die Bahnelemente ausdrückt, sodann differenziert und zwischen 
den Koordinaten und ihren Differentialen eliminiert) erhält man 
die merkwürdigen Ausdrücke: 



dx2 _ k\/p 

dt -^r 
dxj _ ^yp_ 

dt -"/;_ 

dy2 _ kj/p 
dt f,_ 

dr^ _^h^p 



2 

2 (fi— Vi 



^1 — ^1 ^-^ sin* -^ i 

p 

, 2r2Xi 

p 

Vi- Vi j,— sin 



2 



■r,+2r, 



P 



sin 



2 
2 £2— Vi 



2 



(12) 



u. s. f. 
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Hieraus ergiebt sich u. A, 

dXj dx2 d^X2 tx d^Xj t^^ 

Tt dF W'J'^'W^Y 

f^VP ^ 1.(13) 

dx^ dx^ _ d'^x^ ^ d^x^ ^ , 

Tt dt ~"*"d^'T"^ dt^ "T"^"' 

Es verdient bemerkt zu werden, dafs man vermöge dieser 
Gleichungen der Mühe enthoben ist, bereits zur Berechnung 
der Differentialquotienten der 5. und 6. Ordnung neue Orter 
einzuschalten. Doch muis ihrer Anwendung eine Berechnung 
des Parameters vorhergehen. 

Umgekehrt läfst sich aber auch der Parameter aus ihnen 
bestimmen, wenn und soweit man die Differentialquotienten kennt. 
Denn man hat: 

*(^i^2 + ^2^1) sin ^ ^ — 

Endlich mögen noch die Gleichungen 

dXj dt/i 

. Vi^i — x^Vi dt ^ dt ,^^. 

/ 1 k yp 

co8(r,-y.)=^^+^^+^- • • • (16) 

\yz Ti/ y^^^ ^^^^ ^^^^ \ / 

Erwähnung finden. Die letzte Gleichung folgt unmittelbar aus 
der bekannten Eegel für den Kosinus zweier Strahlen (hier r^ 
und ^2) im Räume, indem man bemerkt, dafs 
(f2 — (Ti den Zwischenwinkel von rj und ^2 

— den Kosinus des von ri mit der X-Achse gebildeten Winkels 

u. s. f. . 
darstellt. 
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Entwickelung eines allgemeinen Näherungsrerfalirens. 

Durch Auflösung der Gleichungen (10) folgt: 
Ff, 

/2 



(>2 62 + &34 + 6l4 



F 



/i /i 



(10«) 



Wird die erste dieser Gleichungen unter Mitwirkung der 
Relation (11) in nachstehender Weise umgestaltet: 



?i = «^ - 



0^3 > % »2 " /• 



f 
72 



-=«2 — % — (^ — ^z) 



ti 



2a, n sin ^'^^^^ sin 5^^^ 



^cos 



^3— ^2 



-(a.-a.)(S-|)..(V) 

SO ist die Förtn dieser Gleichung zur Bestimmung von q, sehr 
geeignet. 

Ohschon wir von derselben nur zur weiteren Korrektion 
des aus (HIß) folgenden, bereits bis zu den Gröfsen 4. Ordnung 
genauen Werts Gebrauch machen werden, so wollen wir doch 
mit einigen Worten andeuten, wie dieselbe auch zur Eröffnung 
der gesamten Rechnung dienen kann. 

Setzt man nämlich für die erste Näherung 

n ^ 






" dt ' 2 



, = ^^^...4-...(cf.p.31) 



sin 



(p3—V i „ WpJk+J^ 



cos 



2 
Va— ^2 



2r\ 



= 1, 



so verwandelt sich die Gleichung (V) in: 
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Qi 



__ (oa — ^i) h — (% — %) ^1 ^^^1(^1 + ^2) 



also genau in die erste, auf ganz anderem Wege gewonnene 
Relation des Systems (IV). 

Um nun die Gleichung (V) weiter zu entwickeln und zur 
Ausführung höherer Annährungen brauchbar zu machen, be- 
dienen wir uns der Gleichungen (7), (6) und (9), indem wir 
zunächst folgende Abkürzungen einführen 



sm-' 



2 



«^2, sin2- 



2 




sm 



sm 









'^2; 



Es ergiebt sich dann: 

A h ^ ^/i hfl 

/ 2 ^2 hfl 

o 3 

ti<y-\/p [t«3 — Mit— sin (% ~«2)]~^<**Vi' [w,— %— sin (Mj— «i)] 



<il/VVsin»5^^ 



^2/2 

3^3. 5^^ + 3.5. 7*^ ^•• 



i. y»— yi 



•t,v'r,»Vsm»^ 



4 4.6^ 
3+33 ^»-*-" 



^fi-i» 



2^1 ]/r2r3 sin' 



^ ys— »2 



3^ • cos 



^3— qPa 






(«) 



und nach Substitution dieses Ausdrucks nimmt die Gleichung (Y) 
die Form an: 
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ei = 



(Oj — ffij )^ — («3 — fla) ti 
k 



2nr,r, cos "T^^ cos ^^^I^ cos ^^^f^ • y.y^ 






6 (r/rjra cos 



y»- 



-( 



42^4(^(03 --02) 



bWr^Ts cos 



)V. 



35 



(v^";o7^ip^ 



•[?2--N'?.] 



[I,»-J^l.*]- 



0. Ordnung 



2. Ordnung 



4. Ordnung 



6. Ordnung 



^2^ 



u. s. f. 



• • • • (Va), 

In erster Näherung würde man, wenn die Rechnung 
auch mit vorstehender Gleichung eingeleitet werden soll, 



_ Äyi? ^1 _ Sektor 



= 1 



'^ /i Dreieck 

^2 •= 1 und ya = 1 
^2 = ra = rj 

cos?^-»cos2^==cos^!^ 

setzen und zunächst lediglich die Gröfsen 2. Ordnung mit be- 
rücksichtigen. Da der Faktor 

in diesem Stadium der Rechnung sich auf 



1 — 



^i^ 



reduziert, so erscheint es nur dann unbedingt geboten, das mit 
ihm behaftete Glied 2. Ordnung in die erste Annäherung auf- 
zunehmen, wenn die Zeiträume t^ und ^ eine erheblich ver- 
schiedene Gröfse besitzen. 

In zweiter Annäherung, welche die Gröfsen 4. Ordnung in 
Rechnung zu ziehen hat, müssen insbesondere die Verhältnisse 

^2» ^2» 72 

genauer bestimmt, sowie die Gröfsen (s. Gleichung 8) 
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1 — 



S,= l- 



2v'r,r, 



Ti + r^ 



2Mi' 



cos 



92 — ^1 



-2 



Ti + Ts 






2i/ror, cos 



Vi — 92 



{ 



2-i/r^r. cos 



9a— Vi 



M 



00 



(^0 



2 yy.-^-3--- 2 

berechnet werden. Zu dem Zwecke bedient man sich am ein- 
fachsten der Gleichung (7). Hiernach ist: 



yi 



TiTi sin (qpj - SPi) 



G\/rir2 cos 
und nächstdem: 



2 j 



■yi^ 



^ 6 V 



6.8 



.».»- 



5.7 



W 



•y» 






t^F 



n 



(0 



Anfangs war 
gesetzt worden. 



In zweiter Näherung wird man also 
^2— ^iV 



6(y^i^2 cos 



zu rechnen haben, um die Gröfsen 4. Ordnung vollständig heran- 
zuziehen,' während in den Gleichungen (^) und (l) sowie in der 
rechten Seite der Gleichung (x) die Werte von y^, ^2 vorläufig 
= 1 gesetzt werden, selbstverständlich aber auch in der bereits 
verbesserten Gestalt zur Anwendung kommen dürfen. 

Nachdem so ein genauerer Wert von ^j gefunden, wieder- 
holt sich die Rechnung, indem man die Gleichung (Va) in der 
bis zu den Gröfsen 6. Ordnung erweiterten Form zu Grunde legt. 

Da wir annehmen, dafs die Koordinaten nach der Methode 
der unendlichen Reihen (Gleichung Ulß und 16) schon bis zu 
den Gröfsen 4. Ordnung genau bestimmt sind, so wird man — 
um nun mit Hilfe der (Va) die Annäherung fortzusetzen — 
auch in den rechten Seiten der Gleichungen (/i) und (x) 
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yi = i- 



6(y^.co8^^)' 



zu nehmen und überdies in (x) das mit dem Faktor ^^ versehene 
Glied heranzuziehen haben, um y^ mit dem hier verlangten 
Genauigkeitsgrade zu erhalten. 

Es bedarf jetzt nur noch ein Punkt einer kurzen Erörte- 
rung. Aus dem Systeme (HIß) mit seinen Hilfsgleichungen 
ergeben sich gleichzeitig die Gröfsen 

Qiy ?2> QZJ ^1^ ^2^ ^3 

und dann vermöge der Gleichung (16) und der ihr ähnlichen 
Relationen die Winkel 

92— yij 93 — 9^2» 93— 9l- 

Operiert man nun aber mit der Gleichung (Va) — sei es, 
dafe damit die Rechnung überhaupt begonnen oder die Ver- 
besserung der aus (HIß) gewonnenen Ergebnisse vorgenommen 
wird — so ist noch die Frage zu beantworten, auf welchem 
Wege in diesem Falle die Werte von 

^2 und ?3 
bestimmt werden. 

Transformiert man zunächst die 3. Relation des Systems 

(10a) in ähnlicher Weise wie früher die erste Relation (vergl. 

f 
V und Va), so ergiebt sich — da yj» Yz t~ ^^® ^®^ vorausge- 

/2 

gangenen Rechnung mit der erforderlichen Genauigkeit bekannt 
sind — der korrespondierende Wert von 

?3. 
Weniger bequem gestaltet sich, durch Anwendung und 
gleiche Behandlung der 2. Relation des Systems (10a) die Be- 
stimmung von (>2. Vielmehr ist es vorzuziehen, zuerst ^ mit 
Hilfe einer der beiden anderen Gleichungen zu eliminieren und 
erst dann die resultierende Gleichung in der dargelegten Weise 
weiter zu entwickeln. — Man findet aber auch aus ()| und ^3: 

A _ (h0i—((^i + Qi) (^2 -F asq— (at + gi)(c3+03) ^ 

A 03^2— (63 + ^3)02 fz »2(^3+0 — ^3^2 

und dann durch Substitution dieser Werte in die zweite Relation 
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leicht den Strahl gj* 

Berechnung der Elemente. 

Nachdem die geo- und heliozentrischen Leitstrahlen sowie 
die Zwischenwinkel der letzteren (<f2 — (fi u. s. f.) ermittelt sind, 
berechnet man ki/^t^ 

y. = — ^ 

nach Gleichung (x) voil neuem und erhält daraus sofort den 

Semiparameter p, 
sowie aus (/<) die Gröfse ^i, ferner aus (5) die 
mittlere Entfernung a, 

während die Gleichung « = t/ 1 — ^ die Exzentrizität liefert. 

Diesen Bestimmungen läfst sich aber die Berechnung der übrigen 
auf die Bahn und die Örter des Planeten bezüglichen Gröfsen 
durch Vermittelung der Polargleichungen, der Keplerschen 
Gleichung u. s. f. ohne Schwierigkeit anschliefsen. Es werde nur 
noch daran erinnert, dafs (nach Seite 122) die Bahn-Neigung i 
am einfachsten aus 

QQ^i^M2:?lIZllML (Ebene der Ekliptik « Ebene der y—0) 
U 
gefunden werden kann. 

Auch erscheint die Bemerkung nicht überflüssig, dafs die 
Gleichung (11) 

4r,r,r, sin ^^ siu ^^ sin ^-^ 

trotz aller Eleganz ihrer Form zur Berechnung des Parameters 
im allgemeinen ungeeignet ist. Denn der Nenner, welcher das 
von den drei Planetensehnen gebildete Dreieck darstellt, ist eine 
durch Subtraktion entstandene Gröfse der dritten Ordnung, so 
dafs, wenn die Argumente noch selbst mit Fehlern 3. Ordnung 
behaftet sind, auch nicht einmal ein näherungsweise richtiger 
Wert von p zu erwarten ist. Zwar kann dem Ausdruck auf 
mannigfache Weise eine zur Rechnung geeignetere Form gegeben 



130 



G. Berechnung der Planetenbahnen aus Beobachtungen. 



werden — indem man z. B. Zähler und Nenner durch /i divi- 
diert, wodurch letzterer sich in 

n 

also in eine Gröfse 2. Ordnung verwandelt — doch wird die 
Genauigkeit des Resultats wegen des nicht zu vermeidenden 
Verlusts einiger Dezimalen niemals befriedigend ausfallen. 

Zahlenbeispiel, 

Um dem Leser Gelegenheit zu geben, das im vorstehenden 
entwickelte Verfahren in einem besonderen, der Natur entnom- 
menen Falle zu prüfen, mögen zum Schlüsse die Hauptdata der 
von Gaufs in der Theorie der Bewegung der Himmelskörper 
(§ 151 u. ff.) vollständig durchgeführten Berechnung der 

Juno -Bahn 
hier noch übersichtlich zusammengestellt werden: 
Beobachtungen: 



Beobachtungszeiten 



ß 



legi? 



1804 Okt 5,458644 
» 17,421885 
=» 27,393077 



354« 44' 3r',60 
352° 34' 22",12 
351« 34' 30",01 



— 4ö59'31",06 

— 6o2r55",07 

— 7^17' 50 ",95 



12° 28' 27",76 
240 19' 49"^05 

34^16' 9",65 



9,9980979 
9,9969678 



Hieraus (s.e.c.): ««0,0002051838. 



log aj« 1,2686397 

log -&^«= 1,5847192 

log ci « 1,3187119 



log aj = 1,5353883 

log ~&2= 1,8919579 

log Cj = 1,6489324 



log aa = 1,6657662 

log-&3= 2,0362753 

log Cj- 1,7988222 



ri" = Qx^ + 0,2779226 p^ -h 9,9993652 
ri" - Q2^ 4- 0,2260202 ^2 + 9,9961958 
^2 -u 0,1607421 (»3 + 9,9939356 



Die überstrichenen Zahlen 
bedeuten Logarithmen. 



^3" ^ ^3' 

Verbesserte Zeiten: log ^^ « 1,0778409; log ^2 = 0,9987339. 

Endergebnisse (mit Berücksichtigung der Gröfsen 6. Ordnung) : 

log Ti « 0,3307640 ; log r^ - 0,3259878 ; log r^ = 0,3222239. 

(f^- tp, « 4« 5' 53",34 ; ^3 - ^2 =- 3^ 29' 0",39 

logp« 0,3954836; loga- 0,4224389; log e« log sin 14^2' r',87. 

Neigung der Bahn « 13<»6' 44",10. 

Länge des aufsteigenden Knotens = 171^ 7'48",73. 

Mittlere tägliche Bewegung == 824",7989. 
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